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Robert F. Stengel and Paul W. Berry 
The Analytic Sciences Corporation 


1 . 


INTRODUCTION 


1 . 1 BACKGROUND 

As aircraft become capable of flying higher, faster, 
and with more maneuverability, prevention of inadvertent 
departure from controlled flight takes on added significance. 
To some extent, the airframe can be designed to provide in- 
herent protection against loss of control, as in the addition 
of nose strakes to regulate high angle-of -attack (a) vortices; 
however, performance objectives are likely to dominate the 
choice of such features as wing planform and chord section, 
nose shape, aircraft density ratio, and tail area. It is 
likely, therefore, that the freedom to configure the aircraft 
for intrinsic departure prevention will be restricted and 
that the flight control system will be called upon to provide 
additional protection. 

Flight at high a invariably complicates the control 
problem. Dynamic coupling between longitudinal and lateral- 
directional motions becomes apparent, aerodynamic trends 
vary considerably, and control surface effects diminish or 
become adverse. Coupling and nonlinearities can cause a self- 
sustained oscillation ("wing rock") at high a, degrading 
precision tracking tasks without necessarily causing loss of 
control. Abrupt maneuvering, external disturbances, control 
system failure, or pilot error can produce a "departure" 
(pitch, yaw, or roll divergence), possibly leading to high 
acceleration and to a fully developed spin. The recovery 
from spin or gyration is, at best, an emergency procedure 
which is not always successful. Clearly, it is preferable to 



prevent the departure before it occurs rather than to be 
forced to take emergency measures. 

Although departure and spin are related topics, 
there are at least as many dissimilarities between the 
phenomena as there are similarities. Departure is a trans- 
ient event, while spin is a quasi-steady condition. Depart- 
ure connotes instability with respect to the initial flight 
condition, while spin can be thought of as a bounded, 
periodic (and therefore stable) motion about a near- 
equilibrium flight path. Departure can occur in level 
flight, but spinning equilibrium ultimately results in ver- 
tical motion of the aircraft's center of gravity. 

Departure and spin both are beset by the difficul- 
ties inherent in describing fully coupled dynamic systems 
of high order and in describing static and rotary aero- 
dynamics at complex flight conditions. However, as sug- 
gested by the preceding comparison, the approximations and 
assumptions which hold for one are not necessarily appro- 
priate for the other. In particular, it appears that line- 
arized dynamic models may have a practical utility in pre- 
venting departure which does not readily carry over to spin 
recovery. The reason is that a closed-loop control law 
which continuously acts to prevent departure restricts 
fc .igular excursions to small values; hence, their dynamic 
effects can be described by linear models. A spin recovery 
strategy necessarily must operate with large angular changes 
which result in significant nonlinear effects. 

The key to developing a linear model which is satis- 
factory for the study and control of departure is in the 
recognition that the nominal flight path, used as a refer- 
ence for the variational (linearized) motions, need not 
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t 

) 

represent a steady, unaccelerated flight condition. A flight . ] 

i 

pAth in dynamic equilibrium is equally satisfactory as long ] 

a? it is understood that variational motions are referenced 
to the continuously changing nominal flight* path. The inter- 
mediate step of linearizing the aircraft model with non-zero 
out constant sideslip angle (6q), angular rate (p^, q 0 , r^) , 
and load factor (n^) provides the coupling between longi- I 

tudinal and lateral-directional motions which is so import - 
fnt in the study of departure. 


It is essential to recognize that the combined 
effects of non-zero mean motions lead to significant coup- 
ling which otherwise might be missed in a linear dynamic 
model. It has been demonstrated that large mean values of 
roll rate and sideslip angle separately produce signifi- 
cant coupling of the short period and Dutch roll modes. It 
is less well known that the combined effects of these two 
variables produce coupling which is qualitatively different 
from that induced by a single variable. This variability 
in stability effect is similar to the seeming unpredict- 
ability ui. the departure modes of some aircraft, in which 
the aircraft is known to have more than one departure mode 
for supposedly similar flight conditions. This also sug- 
gests that departure modes are more predictable than might 
have been assumed. 


The complexity of the coupled dynamics and the 
possibility for misinterpreted control cues at high a indi- 

o 

ce*e a need for departure preventing control systems in 
highly maneuverable aircraft. New developments are required 
in characterizing the evolution of motions during extreme 
maneuvering and in the computation of control solutions. 
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1.2 


PURPOSE 


The purpose of this investigation is to identify 
general rules for the design of departure-preventing control 
systems, In achieving this objective, the analytic founda- 
tions for linear-time-invariant modeling of aircraft dynam- 
ics are extended to include extreme maneuvering conditions. 
Using tools of linear systems analysis, the stability and 
control characteristics of a high-performance aircraft are 
examined over a wide range of flight conditions, and spe- 
cific effects of configurational modification are developed. 
The study culminates in the development and evaluation of 
control laws for a Departure-Prevention Stability Augmen- 
tation System (DPSAS) using linear-optimal control theory. 


1.3 SUMMARY OF RESULTS 

The major tasks of this project were defined at 
the outset as: 

• Dynamic Model Development 

• Characterization of Departure Modes 

• Controllability Effects on Aircraft Departure 

• Control Laws for Departure Prevention 

These tasks can be summarized briefly as follows: Dynamic 

Model Development provided a range of nonlinear and linear 
dynamic models for use in the analysis of departure and the 
design of DPSAS control laws. Characterization of Departure 
Modes addressed the unaugmented stability of high-performance 
aircraft. Controllability Effects on Aircraft Departure con- 
sidered the direct (open-loop) effects of control forces on 
aircraft departure. Control Laws for Departure Prevention 
resulted in the design and simulation of linear -optimal 
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regulator control laws which stabilize the reference aircraft 
during extreme maneuvering and which adapt to changing flight 
condition . 

Although aircraft equations of motion are developed 
in Several texts, no derivation which retains all coupling 
terms in the linearized equations of motion was found. 
Therefore, the present investigation began with the develop- 
ment of nonlinear equations in four axis systems (inertial 
or earth-relative, velocity, wind, and body). This was fol- 
lowed by derivation of the associated linear equations of 
motion, as well as equilibrium equations which define a 
generalized trim condition. The validity of these equa- 
tions was established by direct comparisons of the time 
responses of the linear and nonlinear equations. 

A small, supersonic fighter aircraft was chosen as 
a baseline for study. A comprehensive model of subsonic 
nonlinear aerodynamic coefficients was available for this 
aircraft and was used to generate aircraft local stability 
derivatives as flight condition was varied. These sta- 
bility derivatives formed a large part of the linear-time- 
invariant dynamic model, which was analyzed by eigenvalue/ 
eigenvector, transfer function, and time response methods. 
Thus, the significant aspects of aerodynamic nonlinearity 
and inertial coupling, i.e., the local sensitivities to 
initial condition, disturbance, and control perturbations, 
were retained in the aircraft model. 

Using the reference aircraft as a starting point, 
the analysis proceeded along two separate lines. The first 
approach was to assess the effects of maneuvers and flight 
condition on the reference aircraft, recomputing the linear 
model for each variation in nominal angle, angular rate, 
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altitude, and velocity. In order to distinguish between 
aerodynamic and purely inertial effects, a limited number 
of cases were evaluated with varying flight condition and 
fixed aerodynamic derivatives. The second approach was 
to vary individual coefficients of the linear model so that 
specific configurational effects could be analyzed. The 
quantitative results presented here strictly apply only to 
the specific configurations studied; hence, care should be 
exercised in generalizing these results to other configurations. 

The linear-optimal regulator was applied to the 
DPSAS design problem, and the present results demonstrate 
the substantial benefits offered by the linear-optimal con- 
trollers. A design procedure which also identifies gain- 
scheduling relationships is presented; it has the following 
features : 

• Complete longitudinal/lateral-directional 
coupling is accounted for in the design 
process . 

• All significant feedback gains, cross- 
feeds, and control interconnects are 
identified . 

• The control structure is guaranteed to 
stabilize the aircraft, assuming that 
aircraft parameters are known and 
motions are measured precisely. 

• Tradeoffs between control authority, con- 
trol power, and aircraft .notions are 
incorporated in the design process. 

• The DPSAS adapts to varying flight con- 
dition . 

The extension of this design procedure to ?. full command 
augmentation system is direct, as the control design algo- 
rithms are easily restructured to consider handling quali • 
ties requirements, control-actuator rate limits, noisy feed- 
back measurements, and digital implementation. 
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In the process of conducting this investigation, a 
flexible computer program (ALPHA) for the analysis of high 
angle-of-attack stability and control was developed. Pro- 
gram ALPHA generates linear dynamic models and trim condi- 
tions from nonlinear aerodynamic and inertial data; presents 
results for body-axis, stability-axis, and reduced-order 
models; computes eigenvalues, eigenvectors, transfer func- 
tions, departure parameters, and linear-optimal control 
gains; calculates time histories for various initial condi- 
tions, control inputs, and disturbances; and incorporates a 
logical (executive) structure which facilitates parameter 
sweeps, initial condition variations, and model modifica- 
tions during a single computer run. 

1.4 ORGANIZATION OF THE REPORT 

This report presents dynamic equations, stability 
and control charactei istics of high-performance aircraft, 
and control laws for departure prevention. Prior results 
related to extreme maneuvering of aircraft are reviewed in 
C hapter 2 , which then presents a validation of the linear 
■' model and describes the effects of extreme maneuvering on 

the dynamics of the reference aircraft. Configurational 
effects on maneuvering dynamics are discussed in Chapter 3 . 
Control laws for a Departure-Prevention Stability Augmen- 
tation System (DPSAS'i are derived in Chapter 4 , and the 
report is concluded by Chapter 5 . Appendix A is directed 
to the development of nonlinear equations of motion, linear 
equations, generalized trim conditions, and tools for linear 
systems analysis. The model for the reference aircraft is 
summarized in Appendix B. 
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2. DYNAMIC CHARACTEflSTICS OF HIGH-PERFORMANCE AIRCRAFT 


2 . 1 OVERVIEW 

The problems associated with extreme maneuvering 
have two common characteristics : loss of control and large 
angles and/or angular rates , i.e., angles and rates gener- 
ally beyond the range of normal, "1-g" flight operations. 
Extreme maneuvering difficulties fall in the following 
categories, which contain some overlap: 

• Decreased inherent stability 

• Degraded handling qualities 

• Longitudinal/lateral-directional coupling 

• Stall 

• Wing rock 

• Departure 

• Post-stall gyrations 

• Incipient spin 

• Fully evolved spin 

The ordering of this list suggests that the severity of 
these phenomena increases with angle of attack and is 
aggravated by angular rates and sideslip angle. 

Appendix A presents a formal development of fully 
coupled linear-time-invariant models for aircraft motion. 
These models are suitable for investigating perturbation 
motions which are referenced to large angles and large 
angular rates. After reviewing prior investigations of 
aircraft dynamics in Section 2.2 , the remainder of the 
chapter concerns the application of linear systems analysis 
to the stability and control of a high-performance aircraft 
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(which is described in Appendix B) . Section 2.3 compares the 
time responses of linear and nonlinear dynamic models. 

Section 2.4 is directed at aircraft stability, and Sectio n 
2.5 treats aircraft control. Variations in the aircraft's 
dynamic characteristics during extreme maneuvering are 
addressed in Section 2.6 , which also introduces rudimentary 
effects of the pilot's control actions while executing a 
tracking task. The chapter is summarized in Section 2.7 . 


2.2 PRIOR STUDIES OF AIRCRAFT AT EXTREME FLIGHT CONDITIONS 

Although published studies of aircraft dynamics 
shortly followed the Wright brothers' flight (Ref. 1), and 
the concept of stability derivatives was published in 1913 
(Ref. 2), the dynamics of aircraft which are executing 
extreme maneuvers received little attention until the late 
1940' s. (Investigation of the related problem of aircraft 
spinning had begun a decade earlier.) There are several 
reasons for this, but the most significant reason is that 
extreme maneuvers had not presented sufficient problems to 
merit detailed engineering study. The advent of fighter 
aircraft with higher speeds, higher roll rates, higher den- 
sity, lower inherent damping, and higher cost accentuated 
the importance of understanding extreme maneuvering dynamics. 
Furthermore, the improved analytical tools and techniques 
spawned by World War II became available for application to 
flight dynamic problems. 

In addition to the extensive flight testing which 
high-performance aircraft received, three fundamentally dif- 
ferent avenues have been followed in the investigation of 
maneuvering flight. The first approach is the study of 
rigid-body dynamics of the aircraft , the second is the study 
of aerodynamics , and the third is the study of control . The 
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first two areas have a cause-and-ef feet relationship — aero- 
dynamic forces modify the momentum and energy of the air- 
plane — and there is "feedback," in the sense that the 
changing velocity and attitude of the vehicle contribute to 
changes in the aerodynamic forces. Although dynamic problems 
result from the interaction of dynamics and aerodynamics, one 
can distinguish between these two areas in reviewing past 
work. The third area considers methods of augmenting the 
natural aircraft stability, of limiting excursions from the 
normal flight regime, of providing adequate response, and of 
recovering from fully evolved spins. 

2-2.1 Dynamics of the Aircraft 

The objective of study is the solution of nonlinear 
and linear equations of motion, e.g., those derived in Appen- 
dix A. Options for analysis can be classified as explicit , in 
which a direct solution of motion equations is sought, or 
implicit , in which the evolution of motions is inferred from 
characteristics of the system. The solution of these equa- 
tions describes the aircraft's response to initial conditions 
and disturbances, and it provides a basis for identifying con - 
trol policies. The stability of the solution describes its 
tendency to return to a nominal value. Given an initial dis- 
turbance, the stable aircraft's solution returns to the nominal 
solution (or its error is, at least, bounded); the unstable 
aircraft's solution diverges. These analytical methods can be 
summarized as follows: 

Explicit Analysis - Stability, Response, and Control 

• Analog integration of differential equations 

• Numerical integration of differential equations 

• Closed-form solution of differential equations 

• Equilibrium solution of algebraic equations 
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Implicit Analysis - Stability 

• Generalized energy balance (Liapunov method) f 

• Absolute stability bounds (Popov cri- 
terion, circle criterion, etc.) 

• Stability bounds of "classical" dif- 
ferential equations (Mathieu’s equa- 
tion , etc . ) 

• Eigenvalue analysis (Routh-Hurwitz 
criterion, root locus, etc.) 

• Quasi -linear Eigenvalue analysis 

Implicit Analysis - Response and Control 

• Eigenvector analysis 

• Time-domain methods (Impulse or indicial 
response, auto- and cross-correlation 
functions , etc . ) 

• Transform methods (Frequency response, 
transfer functions, spectral density, etc. ) 

Applications of some of these techniques to the 
maneuvering flight problem are documented in the litera- 
ture. Much of the work related to handling qualities, sta- 
bility, coupling, and departure is based upon linear-time- 
invariant models and uses eigenvalue and transfer function 
analysis. Work of this type is reported in Refs. 3 to 10. 

In addition, parameters of linear-time-invariant models 

( C n n j . LCDP , etc.) have been correlated with flight test 
“p,dyn 

or piloted simulation data using little or no direct 
analysis of the equations of motion (Refs. 11 to 14). Quasi- 
linearization of a significant sideslip nonlinearity is 
applied to the wing rock problem in Ref. 15, and closed-form 
solutions for a class of large maneuvers are presented in 
Ref. 16. Since the early 1960's, a large number of investi- 
gations have used analog and numerical integration in the 
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study of departure, stall, post-stall gyrations and spin 
(Refs. 17 to 27). Equilibrium solutions of nonlinear equa- 
tions of motion have been used to determine spin conditions 
and are discussed in Refs. 24 and 28 to 30. 


A number of linear-tiroe-invariant departure param- 
eters have been identified, as reported in Refs. 5, 13, and 
14. These parameters relate to transfer function numerators 
and d nominators and are expressed in terms of stability and 
control derivatives (C ng , Ci g , Cq^, C].^, C n6i> , C^) , 
angle of attack (a Q ), moments of inertia (I x , I z ), direc- 
tional stability augmentation gain (k^ ) , and ailert o -rudder- 
interconnect (ARI) gain (kg): 

Directional Stability Parameter 


* C cos a n - •= — C, sin a n 

n 8,dyn n 6 0 h H 0 


( 2 . 2 - 1 ) 


Lateral Control Departure Parameter 


LCDP = C 

n 




( 2 . 2 - 2 ) 


Augmented Lateral Control Departure Parameter 


LCDP. 

A 





(2.2-3) 


ARI Lateral Control Departure Parameter 


LCDP. dT = C - C. * .. — ^ 

ARI n„ la + k 2 C x 



C +k„ C n 

n 6 2 n 6 
a r 


(2.2-4) 
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Acceleration Response Departure Parameter ("8 plus 6 
Stability Indicator"') 


ardp B = a o 


tan 


-1 


C I 
n 8 X 

C 1 *z 

X 8 Z 


ARDP 6 



(2.2-5) 


(2.2-6) 


The first four criteria indicate resistance tc 
departure when their magnitudes are greater than zero, while 
the last requires ARDP^ to be greater than zero and greater 
than ARDP^. The first criterion relates to the open-loop 
static stability of the Dutch roll mode; C n ^ ( j yn > 0 is an 
approximate requirement for stability. The LCDP's are approxi- 
mations to the closed-loop static stability of the Dutch roll 
mode when lateral control is used to maintain constant roll 
rate; when they ai 3 greater than zero, the Dutch roll mode 
is -vaticaliy stable, but when they are less than zero, the 
Dutch roll mode is statically unstable. The "8 plus 6" cri- 
terion is an attempt to combine stability and control con- 
siderations in a single departure indicator. 

There are a number of inadequacies in the above 
parameters, although they provide insight for future develop- 
ments. They are approximations to the exact transfer func- 
tion coefficients and do not indicate actual pole-zero loca- 
tions; they neglect damping terms entirely; and they do not 
account for longitudinal/lateral-directional coupling induced 
by large . sideslip angle (8) and angular rates ( p,q,r) — 
in. fact, the longitudinal dynamics are ignored completely. 
References 8, 9, 31, and 32 introduce coupling effects due 
to 8, illustrating the importance of some of the neglected 
terms, and Ref. 23 treats the dynamics of steady turning flight. 
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2.2.2 Aerodynamics 

i 

Equally important developments have been made in the 
area of aerodynamics . Measurements of forces and moments and 
visualization of flow phenomena have indicated the large vari- 
ations in aerodynamic conditions to be expected at high angles 
of attack (a) and sideslip (B), and there is an increasing body 
of data related to the effects of large angular rates (Refs. 

34 to 46). While these several references cover a variety of 
topics, they provide an introduction to the kinds of aero- 
dynamic problems which can be expected when aircraft fly at 
high angles and high rates. 

The two dominant phenomena which complicate the 
collection of valid data and the flight of actual aircraft 
are vortices and separated flow . The vortex is a by-product 
of aerodynamic lift, and each surface or body which generates 
lift has a corresponding vortex that trails downstream from 
the lifting source. This swirling airflow affects pressure 
distributions on the downstream surfaces of the aircraft, and 
it can combine with vortices generated on other parts of the 
aircraft to produce a very complex flow field. At low angles 
of attack or sideslip, the vortices from nose, wing, and 
tail usually form a unified flow field, which varies smoothly 
as the aircraft's attitude with respect to the wind changes. 

As the incidence increases, this smooth variation may break 
down, causing the flow to become less unified. 

The difficulties in predicting the actual forces 
and moments on full-scale aircraft from wind-tunnel data 
are greatest for large-angle flight conditions, not only 
because tunnel correction factors can become significant 
but because separation effects depend on the Reynolds num- 
ber of the flow (and, therefore, on the size of the air- 
craft). High-performance aircraft are most likely to perform 
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extreme maneuvers at subsonic velocity, in which case Mach 
number effects may not be significant; however, scaling of 
the flow to provide representative Reynolds number is re- 
quired, if model test data is to be applied to the full-scale 
aircraft . 


Stability problems associated with large aerodynamic 
angles may arise fiorn either the nose, wing, or tail, depend- 
ing on aircraft configuration. Consequently, it is impossible 
to identify a single aerodynamic solution to problems of 
departure (other than to make all aircraft use the same con- 
figuration). Aerodynamic solutions include wing-root leading- 
edge extensions, nose strakes, redesign of the nose cross- 
section and profile, maneuvering (leading-edge) flaps, and 
adjustment of horizontal tail anhedral . 

The aerodynamic forces and moments discussed above 
are static , in that they arise from fixed values of a and 6. 
These terms establish the static stability and trim points 
of the aircraft. Forces and moments which result from 
angular rates (p,q,r) and accelerations (u,v,w,p,q,r) are 
dynamic and thus contribute to damping and transient response. 
There is indication that assumptions which conventionally are 
made for low-angle flight conditions, e.g., that the 6 and 
yaw-rate effects are simply additive, break down at high 
angles. Unfortunately, dynamic forces and moments are dif- 
ficult to measure in practice, and relatively few facilities 
are equipped to measure dynamic forces, much less to separate 
B and r effects. 

2.2.3 Control 

The third subject for study is control of flight 


motions during rapid maneuvering, and it is clear that the 
emphasis of recent studies has shifted away from spin 
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recovery to departure and spin prevention . At best, spin 
recovery is an emergency procedure, and it is not always 
successful. Safety is important, but it is not the only 
issue: an aircraft which is prone to spin is less likely 

to complete its mission successfully. It is preferable, 
therefore, to prevent the spin before it occurs. 

Nevertheless, if a spin occurs, it is important to 
understand what control actions can be used to recover. The 
most favored technique for recovery is to command constant, 
anti-spin control settings (Ref. 47). The proper control 
settings depend on the type of spin (flat, steep, oscillatory, 
or erratic) and on the aircraft configuration — particularly 
the tail damping, aircraft density, and mass distribution. 

In many cases, the available anti-spin control moment is 
less than the restoring moments which maintain spin equi- 
librium, i.e., the spin cannot be broken with constant con- 
trol settings. The idea of "resonating" the aircraft out 
of the spin by applying oscillatory controls was proposed as 
early as 1931 (Ref. 48) and as recently as 1974 (Ref. 49). 
While this task may be difficult for the pilot to execute, 
simple logic for pulsing the controls automatically can be 
designed. 


The concept of automatic control systems which pre- 
vent stall, departure, and spin has gained momentum, and 
it is now recognized that departure prevention can be built 
into the stability augmentation system (SAS), which vir- 
tually all modern high-performance aircraft contain. The 
basic approaches to departure prevention taken to date can 
be classified as limiters (or inhibitors), stability aug- 
menters , control interconnects , or some combination of 
these three. A dual-mode spin-prevention system is developed 
in Ref. 50. This system applies constant antispin controls 
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when a and r exceed separate threshold values, then switches 
to a rate-damping mode once the spin is neutralized. Ref- 
erence 51 presents a departure-prevention system which 
inhibits a, increases directional stiffness (by stability 
augmentation), and restricts the aircraft to roll about its 
flight path. A stall-inhibitor system for a variable-sweep 
aircraft is described in Ref. 5?. This system incorporates 
an a limiter, a- depenaenc command- and stability-augmentation 
gains, increased directional stiffness and damping, and 
aileron-rudder interconnect. A similar philosoDhy is adopted 
in Ref. 53, where speed stability also is augmented to account 
for a-limiting effects in the landing approach. Departure pre- 
vention considerations are evident in the designs for two 
additional high-performance aircraft (Refs. 54 and 55), and 
the effects of stability augmentation and roll/yaw interconnect 
are demonstrated in Ref. 56. 

While a common thread runs through the designs im- 
ported in Refs. 50 to 56, these reports suggest the need for a 
unifying control theory to aid the design of future departure 
prevention systems. These studies have made extensive use of 
experience, nonlinear simulation, and flight testing to arrive 
at successful designs, but the underlying concepts of stability, 
response, and control remain to be identified. 

Summary - This section has presented a brief survey 
of prior developments related to maneuvering flight, dis- 
tinguishing between investigations of dynamics, aerodynamics, 
and control of the aircraft. It is shown that the range of 
problems, from degraded handling qualities to fully evolved 
spin, can not be completely solved by focusing on only one 
area. New developments are required in characterizing the 
evolution of motions; in the measurement and understanding 
of forces and moments at extreme flight conditions; and in 
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the computation of control solutions. In following sections 
of this report, the problems of dynamics and control are 
addressed in detail. 


2 3 COMPARISON OF RESULTS FROM LINEAR AND NONLINEAR 
SIMULATIONS 

The use of linear models in highly dynamic situa- 
tions has been restricted, in the past, by a lack of linear 
models which include complete dynamic effects and by the 
lack of a general method of finding the proper nominal flight 
condition. The linear models developed in Appendix A include 
all the effects of a dynamic nominal flight condition. To 
verify these models and to develop methods of using them, 
this section presents a comparison of linear and nonlinear 
results. The nonlinear results are in .he form of test tra- 
jectories generated by a nonlinear aircraft simulation using 
aerodynamic and mass data for the reference aircraft. 

During the early part of this investigation, large 
differences between the linear and nonlinear results appeared 
along highly dynamic flight trajectories. These were traced 
to the use of an incorrect nominal state vector. From these 
observations, the concept of generalized trim (Section A. 3. 2) 
was developed, and a method of finding generalized trim 
points was derived. (Section A. 3. 2 describes the generalized 
trim calculation computer program.) 

The generalized trim condition is one in which the 
derivatives of the velocity and angular rate states are as 
close to zero as possible. Dimensionality considerations, 
as discussed in Section A. 3. 2, lead to the conclusion that 
the generalized trim problem involves six of the aircraft 
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state equations, the corresponding six states, and six con- 
trol parameters (in this case, four control settings and two 
Euler angles). The problem becomes a search for those values 
of nominal body-axis velocities and angular rates that null 
the selected nominal state rates. 

The following subsections examine specific results 
of the comparison of linear and nonlinear trajectories to 
support these points. 

2.3.1 Elevator Control Input 

Elevator deflection produces a change in pitch 
moment, causing an immediate change in the aircraft angle 
of attack. This causes the aircraft to climb or dive. In 
combination with the throttle, elevator position establishes 
the aircraft flight speed, angle of attack, and flight path 
angle. The tests presented here involve small-amplitude 
elevator inputs when the aircraft is in straight-and-level 
flight at slow speed and high angle of attack. Figure 2.3-1 
illustrates the time history of the most important longi- 
tudinal motion variables for eight seconds following the 
control application. All lateral variables are approxi- 
mately zero for the nonlinear model and exactly zero for the 
linear model. 

Comparison of the linear and nonlinear curves indi- 
cates excellent agreement. It is important to note that the 
nonlinear aircraft response verifies that the lateral and 
longitudinal modes are truly uncoupled in this flight con- 
dition. The nominal flight condj.cion is a steady-trim flight 
condition and satisfies the generalized trim condition. 
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Figure 2.3-1 Small Amplitude Elevator Input 


2.3.2 Aileron Control Input 

The ailerons primarily provide roll moment, and 
the trajectories shown in Fig. 2.3-2 illustrate the air- 
craft response to a small amplitude aileron doublet. The 
linearized trajectory, whose nominal flight condition 
is again straight-and-level flight, differs only slightly 
from the true nonlinear response, and the linear and non- 
linear trajectories exhibit lateral-longitudinal separation. 


2.3.3 Rudder Control Input 

Large-input, large-response trajectories resulting 
from rudder deflection are examined in this subsectiou, with 
the goal of testing the trajectory matching capabilities of 
a linear simulation for a highly dynamic flight condition. 
The nonlinear test trajectory lasts eight seconds after the 
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Figure 2.3-2 Small Amplitude Aileron Input 


control is applied; linear trajectories starting at the 
initial time and at four seconds into the trajectory are 
tested . 


Figure 2.3-3 compares the nonlinear trajectory to 
a linear trajectory starting at the time of control appli- 
cation. The nominal trajectory for linearization is the 
original static trim flight condition of straight-and-level 
flight. The trajectory match is acceptable for about two 
seconds, and the angle-of-attack plot illustrates the cause 
of the deviation. Because it exhibits lateral-longitudinal 
separation, the linear trajectory does not capture the change 
in angle of attack that the nonlinear trajectory contains. 

This change in angle of attack has a large effect on the 
subsequent dynamics which the linear model fails to duplicate. 
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By examining the trajectory beginning four seconds 
after the control is applied, methods of linearization for 
highly dynamic trajectories can be derived. Figure 2.3-4 
illustrates an early attempt. Here, the point of lineariza- 
tion is approximate, i.e., it does not satisfy the generalized 
trim condition discussed below and in Appendix A. The re- 
sulting linear trajectory diverges from the nonlinear tra- 
jectory fairly quickly, and the slopes do not match at the 
initial point for some states. Furthermore, the frequency 
of the resulting motion is considerably different from that 
of the nonlinear motion. Due to its dependence on ad hoc 
estimation cf the nominal flight condition, the results of 
this approach are highly variable in quality. 

One of the most striking errors in the linear tra- 
jectories shown in Fig. 2.3-4 is that the slopes of the 
states do not match at the beginning of the linear trajec- 
tory. This observation, which implies that the nominal 
state r ates are not zero, led to the development of the 
generalized trim concept. In this context, this concept 
indicates that to provide an accurate representation of a 
nonlinear system by a linearized one, it is necessary to 
choose a point of linearization that exhibits zero nominal 
state rates . 

Applying this generalized trim procedure to the 
point four seconds after control application produces the 
results shown in Fig. 2.3-5. Compared to the previous fig- 
ure, the generalized trim procedure produces clearly superior 
results. There are no initial slope errors evident, the 
match is excellent for two seconds, and it is reasonably 
close for much longer. Additionally, the frequency of the 
linearized motions is close to that of the nonlinear 
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trajectory. This, along with the amplitude match, supports 
the use of a properly linearized model for the analysis of 
a nonlinear vehicle along a highly dynamic trajectory. The 
significance of this result is put in proper perspective 
when it is realized that the vehicle has perforim d a 360-deg 
roll between t * 0 and t * 7.2 sec, and the pitch angle goes 
from 15 deg to -45 deg from t * 0 to t = 5 sec. 

It should be noted that the nominal flight condi- 
tion for linearization was found by an analytic method that 
does not require the solution of a nonlinear trajectory 
from which to estimate nominal values. The generalized trim 
procedure is a useful method for calculating nonunal flight 
conditions even along highly dynamic trajectories. 

Summary - These comparisons presented here establish 
that nominal flight conditions which satisfy the generalized 
trim condition produce good trajectory matches and that the 
corresponding linear models should provide accurate infor- 
mation about the nonlinear system dynamics. 

2.4 EFFECTS OF ANGULAR MOTION AND FLIGHT CONDITION ON 
AIRCRAFT STABILITY 

The effects of altitude and velocity variations, 
angle-of-attack and sideslip angle variations, and steady 
angular rates on aircraft stability are examined in this 
section using the linearized dynamic models and eigenvalue/ 
eigenvector analysis technique presented in Appendix A. The 
purpose of this analysis is to show the effects of indivi- 
dual flight variables, as well as the combined effects of 
flight variables which normally are zero in "1-g" straight- 
and-level flight. For this study, the aircraft is trimmed 
initially for "1-g" flight at an angle of attack of 15 deg 
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and at an altitude of 6100 m. As flight variables 
change, the lead factor may change accordingly; however, 
the primary objective of this chapter is to isolate the 
individual effects of each specific flight variable being 
examined, so all other variables are held at their initial 
values. 

2.4.1 Altitude and Velocity Effects 

Altitude affects the air density and, therefore, the 
dynamic pressure. This causes the aerodynamic forces and 
moments to be reduced, relative to the inertial effects, 
as altitude increases, as shown in Fig. 2.4-1. Higher 
altitude causes both the natural frequencies and damping ratios 
of the Dutch roll and short period modes to decrease. The roll 
mode also slows down as altitude increases. 

Changes in velocity affect the dynamic pressure, as 
well as the angular rate normalization terms (b/2V and c/2V) 
and the velocity-angular rate cross-product terms. These 
changes cause significant increases in Dutch roll and short 
period frequencies as velocity increases (Fig. 2.4-la). The 
damping ratio of the short period mode is affected only 
slightly by velocity changes over the range shown in Fig. 
2.4-lb. The Dutch roll damping decreases as velocity 
decreases, so that the Dutch roll is unstable at the lowest 
velocities presented here. The reference aircraft's roll mode 
(Fig. 2.4-lc) is changed only slightly as velocity varies, 
contrary to the result obtained from the approximate lateral- 
longitudinal equations discussed below. 

Figure 2.4-1 indicates only small increase in 
spiral mode speed and phugoid frequency and damping at 
lower altitudes. Low velocity results in low phugoid 
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damping, so much so that the mode is unstable over a signifi- 
cant portion of the velocity range examined in Fig. 2.4-1. 


As noted above, the lack of roll mode variation with 

0 

velocity is contrary to results obtained with accepted 
approximations. The expected linear change in roll mode with 
velocity is deduced from approximate lateral-directional 
equations, which can be derived by neglecting the roll angle 
equation (and, therefore, the spira’ mode), by assuming that 
the Dutch roll mode consists of wind-axis yawing motion, and 
by assuming that the roll mode consists of wind-axis roll. 

The approximations that result from these approximate lateral- 
directional equations are 
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where A ^ is ^e r °H m °de eigenvalue, 
as in Eq. (2.2-1) and 
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Equations (2.4-1) and (2.4-2) predict values for 
A ro j^ and £o) n DR as shown by the dotted lines in Fig. 2.4-ld 
for the 6100-m case. The actual values, taken from Figs. 
2.4-lb and c, are shown by solid lines on the same figure. 
The approximate equations do a poor job of predicting mode 
speed because the subject aircraft is fuselage-heavy (high 
I /I v ratio); hence, the Dutch roll contains more rolling 
response than is assumed when deriving the approximate 
equations. As can be seen from Fig. 2.4-ld, this leads to 

o 

a damping interchange such that the roll mode is faster than 
expected and the Dutch roll mode is more poorly damped than 
predicted. 

This examination leads to the following conclusions 
for the reference aircraft: 


• Higher altitudes result in lower damping 
and frequency of the Dutch roll and short 
period modes, as well as increased roll 
mode time constant. 

• Lower velocities result in a decrease in 
short period frequency at constant damping 
ratio, as well as decreased Dutch roll 
frequency and damping. 

• The approximate lateral-directional equa- 
tions do not predict roll mode or Dutch 
roll damping accurately for the subject 
aircraft. The complete equations should 
be used for an accurate determination of 
these parameters. 


2.4.2 Aerodynamic Angle Effects 


The aerodynamic angles, a and g, specify the orienta- 
tion of the vehicle relative to the velocity vector, and, to a 
large extent, they define the flow field around the vehicle. 
For this reason, the aerodynamic angles are prime determi- 
nants of- the aerodynamic forces and moments. Consequently, 
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significant differences in the speeds and shapes of the 
normal modes occur as a Q and Bq are varied. 

Figure 2.4-2 illustrates the boundaries between 
stability and instability which result from these variations. 
These boundaries define the and Bq for which the real 
part of one or more eigenvalues migrates from negative 
(stable) to positive (unstable) sign (see Section A. 4.1). 

The phugoid mode is a slow mode and is unstable at low ciq. 

The Dutch roll mode, a fast mode, becomes unstable at high 
Gq. The dashed line in Fig. 2.4-2 is an important boundary, 
indicating the transition of a relatively slowly divergent 
phugoid oscillation into two real roots, one of which is 
highly unstable. This transition line occurs at high Bq — 
about 10 to 15 deg. 



Figure 2.4-2 Effects of Aerodynamic Angles on 

Aircraft Stability 

The shape of the Dutch roll stability boundary 
indicates that moderate values of nominal sideslip angle 
(two to five deg) stabilize the mode. This is due to 
lateral-longitudinal coupling; a close examination of the 
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Dutch roll/short period eigenvalues indicates that, up to 
about five deg of sideslip, Dutch roll damping increases as 
short period damping decreases. 

The eigenvectors of the linearized model provide 
information about the normal mode shapes which indicate the 
involvement of each state in each mode. Figure 2.4-3 illus- 
trates some specific eigenvalue/eigenvector variations with 
angle of attack. Real eigenvectors, such as those associated 
with the roll mode, are characterized only by the relative 
magnitudes of each state, as the phase angles are either 0 
or 180 deg. A time history of this mode would show a con- 
stant ratio between the various state amplitudes. These 
amplitudes would evidence exponential decays with equal time 
constants, given by the negative inverse of the eigenvalue. 
Complex eigenvectors, such as those of the Dutch roll, are 
characterized by the relative magnitudes of the involved 
states and by the phase angle between them. A time history of 
this oscillatory mode is generated by the projections of the 
eigenvectors on the real axis as the entire eigenvector set 
rotates with angular rate given by the imaginary part of 
the eigenvalue. The magnitudes decay exponentially with the 
time constant given by the negative inverse of the eigen- 
value's real part. 

Despite the large changes in eigenvalues with angle 
of attack, Fig. 2.4-3 shows little corresponding change in 
eigenvector shape. The only major changes involve the pro- 
portions of angular rates in the fast modes, and these changes 
are due to mode speed variations, as described above for the 
roll mode. The short period mode contains increased pitch 
rate at large ciq for this reason; in the Dutch roll 
eigenvector, the roll rate-to-sideslip ratio increases and 
decreases with the Dutch roll frequency. Overall, the Dutch 
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Figure 2.4-3 Angle-of -Attack Effects on Aircraft Eigenvectors 


roll mode of this aircraft involves a great deal of rolling 
motion, underlining the low rolling inertia typical of modern 
fighters. 


The short period eigenvector shows that this oscil- 
lation typically 'involves angle-of-attack perturbations at 
constant velocity, as axial and normal velocity perturbations 
are approximately 180 deg out of phase with each other and 
are related in magnitude by tan c*q. The short period mode 
is faster at high angle of attack, and it includes more pitch 
rate than at low angle of attack. 

The changes in specific eigenvalues and eigenvectors 
with sideslip angle are illustrated in Fig. 2.4-4. Lateral- 
longitudinal coupling is quite prominent for the asymmetric 
flight conditions portrayed in this figure. Modes of com- 
parable speed couple most readily. Roll angle response is 
found in the phugoid mode, and pitch angle becomes a component 
of the spiral mode, so that both modes involve slow roll-pitch 
motion. Angle of attack appears in the Dutch roll eigenvector, 
and a roll-sideslip combination becomes important in the short 
period mode, so that both modes involve an angle of attack- 
sideslip oscillation. In both cases, Aw and Av (or, equivalen- 
tly, Act and A0) are almost 180 deg out of phase. Note that 
the changes in the speeds of these modes are small and gradual 
as the sideslip angle is varied. 

To demonstrate some of the causes of aerodynamic 
angle effects observed above, the aerodynamic coefficients 
are held constant (at the values for a q * 15 deg and 6 q * 

0 deg), and the body orientation with respect to the velocity 
vector is varied over the same range of aerodynamic angles 
used in Fig. 2.4-2. The results, shown in Fig. 2.4-5, differ 
significantly from those shown in Fig. 2.4-2. There is only 
a slight sideslip effect. The Dutch roll mode, rather than 
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ANGLE OF ATTACK, a 0 (DEG) 

Figure 2.4-5 Effects of Body Orientation on 
Aircraft Stability 

becoming unstable at high olq , is destabilized by lower a Q . 

The phugoid stability boundary near olq =* 12 deg is roughly 
similar to that found in Fig. 2.4-2, indicating that the lack 
of phugoid stability in this area is not due to aerodynamic 
variations . 

Figure 2.4-6 assists in the evaluation of the high 
angle-of-attack Dutch roll instability. This figure com- 
pares the Dutch roll eigenvalue to the departure parameter 
C ne dyn (see Section 2.2) and to Cn r and C nj . dyn , the last 
of which is defined in Eq. (2.4-5). 

The results of Fig. 2.4-6 indicate that, at least 
in this case, C n g dyQ is a good indicator of the Dutch roll 
mode's imaginary part. Neither C nr nor C nr dyn provide a 
particularly useful indication of Dutch roll stability. 

This example indicates that Cng ( j yn has only limited value 
as a departure parameter. For this aircraft, Dutch roll 
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Figure 2.4-6 Variations of Directional Aerodynamic 
Coefficients with Angle of Attack 

instability is due to negative damping, and C n g ^ is in- 
adequate as a predictor of departure. 


The following conclusions concerning aerodynamic 
angle effects on stability of the reference aircraft can be 
made: 


• The Dutch roll mode becomes unstable due 
to negative damping at high otQ. This is 
caused by changes in the aerodynamics as 
ao increases. 

• Mean angle of attack variations have sig- 
nificant effect on the eigenvalues, but 
mode shape (eigenvector) changes are 
small relative to other effects. 

• Mean sideslip angle introduces lateral- 
longitudinal coupling; therefore it 
has a large effect on the mode shapes 
(eigenvectors), without causing large 
changes in the eigenvalues. This lateral- 
longitudinal coupling primarily occurs 
between modes of similar speed and can 
lead to a transfer of damping, as in the 
situation where small sideslip angles 
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stabilize the Dutch roll mode at the 
expense of short period damping. 

• The parameter C' n g ^ gives a good indi- 
cation of Dutch roll frequency, but it 
is not useful as a departure parameter 
for the subject aircraft. 

2.4.3 Angular Rate Effects 

Non-zero nominal angular rates have two effects on 
the linearized aircraft uynamics. The first, an aerodynamic 
effect, results in a change in the nominal forces and moments 
due to the steady angular rates. The second is dynamic, and 
it is due to the cross product of angular rate with velocity 
(in the force equations) and with angular momentum (in :he 
moment equations). The specific terms involved (for I 3 0) 

XZ 

are given in Table 2.4-1. A close examination reveals that 
mean pitch angular rate, q Q , enters both the lateral and 
longitudinal equations but does not affect lateral-longitu- 
dinal coupling terms. Mean roll and yaw rates, p Q and r Q , 
enter as lateral-longitudinal coupling terms. Steady roll- 
rate capability of most high-performance aircraft is much 
higher than pitch- or yaw-rate capability, so roll-rate 
effects are especially important. 

Stability boundaries as functions of pitch rate and 
yaw rate are illustrated in Fig. 2.4-7. The destabilizing 
influence of q^ is the major effect, and it has an especially 
severe effect on the Dutch roll mode. Yaw rate has a mild 
stabilizing effect on the Dutch roll and spiral modes. This 
is due partially to lateral-longitudinal coupling, because 
short period and phugoid damping decrease as Dutch roll damping 
increases . 


38 


mm 


v 


TABLE 2.4-1 

DYNAMIC EFFECTS OF STEADY ANGULAR RATE 


Angular Rate 

Multiplied By 

Enters Term 

p o 


3q/3r 

p o 

-l 

3w/3v 

p o 

l 

3v/3w 

p o 

(I x -Iy)/Iz 

3r/3q 

q 0 

-1 

3u/3w 

q 0 

1 

3w/3u 

q 0 

dx-Iy)/Iz 

3r/3p 

q o 

(Iy-Iz)/Ix 

3p/3r 

r 0 

1 

3u/3v 

r 0 

(*z~*x)/*y 

3q/3p 

r 0 

-1 

3v/3u 

r ° 

( ly-^z^/^X 

3p/3q 


P7T3 DUTCH ROLL STABLE 
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Yaw-Rate/Pitch-Rate Effects (a Q *15 deg) 


Figure 2.4-7 






The eigenvector changes that accompany increases in 
are shown in Fig. 2.4-8. Steady pitch rate does not 
introduce lateral-longitudinal coupling, but some changes in 
mode shapes appear in the roll mode and in the separation 
of the complex phugoid mode into two real roots. Both fre- 
quencies and damping ratios of the Dutch roll and short 
period modes change, but changes in the mode shape are minor. 

Steady roll rate is important because fighter air- 
craft are capable of higl. Pg, and air combat maneuvers often 
include such motions. For the aircraft to roll with con- 
stant aerodynamic angles, the roll rate must occur about the 
wind x-axis (which is the same as the stability x-axis for 
constant nominal aerodynamic angles). Sideslip variations 
also are considered, since piloting error can easily result 
in non-zero Bq during a rolling maneuver. Both positive and 
negative Pq are considered, to account for roll "into" or 
"out of" the sideslip. 

The stability boundaries that result from combined 
ro? .I rate and sideslip are shown in Fig. 2.4-9. These 
boundaries indicate that Pg has only a small effect on the 
fast modes, primarily the Dutch roll mode. The combination 
of Pg and small values of Bq of opposite sign serves to 
destabilize the Dutch roll mode. Roll rate destabilizes 
the phugoid mode in general, but there is a combination of 
Bq and Pg that maintains phugoid stability. High Bq results 
in a fast divergence for all values of Pg tested. 

Eigenvector variations due to steady rolling are 
illustrated in Fig. 2.4-10. Eigenvalue changes are signi- 
ficant, considering the angular rates involved. The mode 
shapes also change, so that lateral-longitudinal coupling 
is important. Large roll-rate/sideslip perturbations in 
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Figure 2.4-9 Stability Boundaries for Sideslip/ 

Roll-Rate Variations (<*q = 15 deg) 


the short period mode and large angle of attack perturba 
tions in Dutch roll mode are examples of this coupling. 


Conclusions about steady angular rate effects are 
as follows: 


• # Mean yaw rate and roll rate cause 
lateral-longitudinal coupling and 
therefore change the mode shapes 
significantly. Roll rate is by far 
the more significant because of the 
large values it can exhibit. 

•* Mean pitch rate changes the speeds of 
the normal modes without affecting 
their shapes significantly. Even low 
values of qQ (about 5 deg/sec) can cause 
the Dutch roll mode to be unstable. 
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Figure 2.4-10 Roll-Rate Eigenvector Effects 


2.5 


EFFECTS OF ANGULAR MOTION AND FLIGHT CONDITION 
ON AIRCRAFT CONTROL 


The transfer function provides a primary measure 
of the quality of aircraft control, as it is the Laplace 
transform of the ratio between a specific output and a spe- 
cific input (Appendix A). The transfer function gain, Kp, 
is the stady-state value of the transfer function after all 
transients damp out, asstuning that all transients are stable. 
The transfer function gain, Kj, is (for most aircraft), the 
initial state rate response to a transfer function's control 
step. The poles of the transfer function are the eigenvalues 
of the unforced system, as described in Appendix A. 

The zeros affect the magnitudes of excitation of 
the normal modes, which are related to the distance between 
the zeros and the appropriate eigenvalues in the s plane. 

In the limiting case, a zero and pole in the same location 
cancel, and the. corresponding mode does not appear in that 
response. Zeros located in the right half-plane are called 
nonminimum-phase zeros (due to their effects on the phase- 
shift of sinusoidal inputs), and they have major impact on 
the aircraft's transient response and on controller design. 
For example, an undesirable reversal in the initial res- 
ponse is caused by such zeros, as illustrated in Fig. 2.5-1. 

The nonminimum-phase type of response is undesirable 
because it makes closed-loop control difficult. The pilot 
can be misled by this type of response, as the magnitude and 
sign of the motion are uncertain. Addition of a high-gain 
feedback loop around a transfer function that exhibits non- 
minimum-phase properties can result in instability of the 
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Figure 2.5-1 Typical Step Response Forms 


closed-loop system (Ref. 9). Finally, this type of response 
can make it impossible to implement some simple forms of 
adaptive control, as they can suffer from instability for 
an analogous reason (Ref. 57). 

2.5.1 Velocity and Aerodynamic Angle Effects 

As has been observed previously (Section 2.4.1), 
velocity changes the dynamic pressure, which affects the 
control effectiveness. Table 2.5-1 illustrates some typical 
values of the transfer function gain, Kp, at different speeds, 
and the variation is as expected. Altitude also affects 
dynamic pressure in that increasing altitude decreases 
-atmospheric density; hence, dynamic pressure decreases. 


45 



TABLE 2.5-1 

VELOCITY EFFECTS ON TRANSFER FUNCTION GAIN, K p 


Velocity 

A w/AS h 

Ap/A6 a 

Ar/A6 a 

70 m/s 

-3.88 

-0.50 

o 

« 

1 

94 m/s 

-5.17 

-0.66 

-5.39 

117 m/s 

-6.47 

-0.83 

-6.74 


Aerodynamic angle variations can cause large changes 
in the system eigenvalues and can be expected to have sig- 
nificant effects on the numerator of the transfer function as 
well. Table 2.5-2 illustrates variations in Kj and Kp as 
olq and 0 q vary. (See Section A. 4. 3 for the definition of Kj . ) 
The invariability of Kj with sideslip indicates that the 
control effectiveness does no; depend on sideslip. This is 
a function of the aerodynamic data used here (Appendix B) , 
as the data does not model the effects of on control 
effectiveness. The steady-state gain, Kp, shows a large 
dependence on sideslip because this gain depends on the pole 
and zero locations, which themselves vary with Sq. 

The sign changes in Kp (at small c*q) as 8 q varies 
are due to changes in the number of unstable poles and non- 
minimum-phase zeros. From Eq. (2.5-6), it can be seen that 
such a change results in a Kp sign change if the sign of Kj 
remains the same. 

One of the major effects of angle-of-attack varia- 
tions (best seen in the X gain of theAp/A'S. transfer 

I «- 
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TABLE 2.5-2 

AERODYNAMIC ANGLE EFFECTS ON TRANSFER FUNCTION GAINS (V Q =94 m/s) 




k f 


ci (deg) 

A (deg) 


& V /A6 a 



Av/A£ r 

Ar/i6 r 

tp/u T 


6P/A« a 

Ar/A6 r 

5.0 

0.0 

- 1.04 

2.46 

18.39 

-4.91 

- 0.95 

-7.11 

1.90 

- 0.10 

2.98 

-0.53 

5.0 

10.0 

0.33 

- 3.86 

•32.02 

6.64 

0 . 3 / 

3.04 

-0.82 

- 0.10 

2.98 

-0.53 

15.0 

0.0 I 

- 1.58 

0.41 

2.47 

-0.66 

- 0.89 

-5.39 

1.44 

- 0.14 

1.71 

-0.57 

15.0 

10.0 

- 2.37 

1.43 

9.88 

-2.64 

- 1.71 

-11.85 

3.17 

- 0.14 

1.71 

-0.57 

25.0 

0.0 I 

- 0.84 

0.03 

0.49 

-0 13 

- 0.27 

-4.29 

| 1 .15 

- 0.11 

-0.001 

-0.44 

25.0 

10.0 

- 1.50 

0.01 

0.42 

-0.11 

- 0 . C 5 

-3.83 

[ 1.02 

- 0.11 

-0.001 

-0.44 


function) is the loss of aileron roll control at high a Q . This 
loss of aileron roll effectiveness, combined with the continued 
effectiveness of the rudder for roll and yaw control, leads 
to the conclusion that this aircraft is rolled more 
effectively with the rudder at high angles of attack. 

An examination of the transfer function zeros 
(Table 2.5-3) indicates that nonminimum-phase zeros are 
quite prevalent, although often accompanied by right -half- 
plane poles, i.e., they often occur in unstable systems. 
Right-half-plane zeros are important in control system 
design because closed-loop poles of a system with a simple 
loop closure migrate from the open-loop poles to the zeros 
as the loop gain is increased. Therefore, in a system with 
right-half-plane zeros, too high a gain may result in an 
unstable closed-loop system. 

When is not zero, there is mode coupling, and 
a control input excites all modes. This is indicated by 
Table 2.5-4, which presents the poles and zeros of three 
transfer functions at a flight condition where 8 0 is non- 
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TABLE 2.5-3 

EFFECTS OF ANGLE OF ATTACK ON TRANSFER FUNCTION ZEROS 


Transfer 

Function 

« 0 (d#g) 



Zeros 


Aw/A« h 

5 

0.0084 


0.0824 

-25.89* 

U 

IS 

•0.0025 

0.1167 

-21.07 


25 

-0.0281 


0.1500 

-25.56 

Av/AA. 

5 

•0.168 


-1.465 

86.43 


15 

•0.137 


•0.512 

1325.0 


25 

•0.164 


0.178 

133.7* 

Ar/A5 

5 

0.350 


2.914 

•1.06 

ft 

15 

2.06 


9.63 

-0.385 


25 

0.036 

* j 

1.725 

-0.237* 

AP/A6. 

5 

•0.141 


1.399 

0.0270 

ft 

15 

-0.152 

ij 

1.295 

0.0278 


25 

•13.67 


22.49 

0 . 0300* 

Av/A«„ 

5 

0.0198 


•0.840 

-40.97 

r 

15 

0.0956 


-0.393 

-45.74 


25 

0.0795 


•0.188 

-41.70* 

Ar/A6 

5 

0.030 

»J 

0.686 

-0.910 

* 

15 

•0.025 

ij 

1.720 

-0.391 


25 

0.023 


1.694 

-0.238* 


5 

0.027 


-1.80 

1.79 

r 

15 

0.028 


•2.98 

3.54 


25 

0.030 


-2.66 

3.05* 


•Accompanied by right-half-plane poles 


TABLE 2.5-4 

POLE-ZERO COMPARISON AT ct Q » 15 DEG, 8 Q - 10 DEG 


Poles 

Short Period 

Dutch Roll 

Roll 

Spiral 

Phufold 

-0.461 tj 0.915 

-0.0695 tj 1.953 

-0.390 

-0.026 

-0.0227 tj 0.153 


Zeros of 

-30.75 

-0.0938 tj j.033 

-0.435 

-0.033 

*0.0009 tj 0.114 

| Zeros of Ar/A* r | 

-0.593 tj 1.001 

0.153 *J 1.301 

-0.391 


-0.0313 tj 0.147 

Zeros of ip 

-0.449 tj 1.161 

-0.131 tj 1.373 


0.025 

-0.0204 tj 0.143 
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zero. Note that the lateral mode poles are not canceled in 
Aw/Afi^ and that the longitudinal mode poles are not canceled 
in the Ar/A6 r and Ap/A6 & transfer functions. 

The effects of velocity and aerodynamic angles on 
control of the example aircraft can be summarized as follows: 

• Lower velocities lead to decreased control 
effectiveness, as demonstrated by transfer 
function gains. 

• Non-zero Bo does not affect Kj , but does 
change the poles and zeros so that all 
modes are excited. 

• Mean angle of attack leads to significant 
changes in control effectiveness, so much 
so that the rudder is more efficient than 
the aileron for producing roll at high 
angles of attack. 

• Nonminimum-phase zeros are prevalent in 

the aircraft transfer functions at increased 
angle of attack. 

2.5.2 Angular Rate Effects 

Although no explicit effects of nominal angular 
rates on the control effectiveness are included in the 
specific aerodynamic data used here, angular rates cause 
significant changes in the transfer functions due to pole 
and zero shifts. This is apparent in Table 2.5-5, which 
shows changes in transfer function gains due to nominal 
pitch rate. The initial value of the transfer function, Kj , 
does not vary with because Kj depends only on the con- 
trol effectiveness. The steady-state gain, Kp, does vary 
with qg because of the pole and zero variations. As above, 
sign variations in indicate the appearance of unequal 
numbers of nonminimum-phase zeros and right half-plane poles. 
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TABLE 2.5-5 

EFFECTS OF PITCH RATE ON TRANSFER FUNCTION GAINS (a o =15deg) 



k f 



K i 


Aw/ A«5 h 

A v/A« a 

ar/A4 a 

Ap/A4 a 

Av/A5 r 

Ar/A<5 r 

Ap/A4 r 

A«f/A4 h 

tr/L6 r 

Ap/A6 a 

0 

-I.3S 

0.41 

2.47 

-0.66 

-0.89 

-5.39 

1.44 

*3* 

r-* 

o 

1 

-0.58 

1.71 

12 

-5.67 

-0.02 

-2.61 

2.44 

0,04 

5.68 

-5.32 

-0.14 

-0.57 

1.71 

24 

-2.51 

1 

0.17 

-0.39 

0.72 

-0.37 

0.84 

1.56 

» 

o 

h* 

*• 

-0.57 

1.71 


Mean wind-axis roll rate, like has no effect on 
control power but does change the mode shapes significantly. 

The steady-state transfer function gain varies with p^ as 
shown in Table 2.5-6. This variation is fairly smooth, com- 
pared to the effects of q QI and the only sign changes are 
in the roll-rate transfer functions. Because non-zero roll 
rate creates lateral-longitudinal coupling, any control dis- 
placement excites all of the normal modes. 

TABLE 2.5-6 

EFFECTS OF ROLL RATE ON TRANSFER FUNCTION GAIN, K p (a *15 deg) 


T <0244 


P w (deg/ sec) 
o 


Av/A<$ a 

&r/A5 a 

Ap/A6 a 

A v/ A <5 

r 

A r/A6 r 

Ap/A6 r 

0 

-1.58 

0.41 

2.47 

-0.66 

-0.89 

-5.39 

1.44 

15 

-1.64 

0.24 

0.69 

-0.11 

-0.52 

-1.49 

0.23 

30 

-1.40 

0.22 

0.32 

0.06 

-0.51 

-0.73 

-0.13 

45 

-1.25 

0.22 

0.18 

0.10 

-0.62 

-0.51 

-0.29 

60 

-1.23 

0.22 

0.10 

i 

0.04 

1 

-1.34 

-0.62 

-0.27 


Conclusions about angular rate effects on the con- 
trollability of the subject aircraft are as follows: 

• Non-zero nominal angular rates do not 
change control effectiveness, but do 
change mode shapes and/or speeds, as 
well as zero locations. 
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• Nominal qo has a large effect on Kp, 
primarily due to the creation of non- 
minimum-phase zeros and unstable modes. 

• Nominal po primarily changes mode 
shapes rather than pole locations, 
and it causes any control de- 
flection to excite all response 
modes . 


2.6 DYNAMIC VARIATIONS DURING EXTREME MANEUVERING 

Aircraft may be especially prone to departure from 
controlled flight during air combat maneuvering because such 
maneuvers are executed using the highest possible aircraft 
performance, and pilot workload during maneuvering flight is 
high. Although it is possible to fly most combat maneuvers 
in a smooth, coordinated manner, even small errors can cause 
difficulty due to instability, unfamiliar coupled mode shapes, 
or changes in control effectiveness. 

Many air combat maneuvers include periods of high 
angle-of-attack flight, in order to produce a large normal 
force for climbing or turning. High angular rates also are 
typical of many air combat maneuvers. High normal accelera- 
tion may be accompanied by large q Q , and large p Q may be 
generated to rapidly orient the lift force in a desired 
direction. 

Referring to the earlier sections in this chapter, 
the difficulties involved in extreme maneuvering become 
clear. High angles of attack and pitch rate destabilize the 
normal modes of motion and reduce the available control 
power, while high roll rate causes lateral-longitudinal 
coupling and produces mode shapes unfamiliar to the pilot. 


i 


i 

} 

v 


51 



1 


i 


t 


The first two of the following sections examine the 
changes in aircraft stability and control along two typical 
air combat trajectorie The third section approaches the 
same problem from a different viewpoint, examining the effects 
of an elementary target-tracking pilot model on aircraft 
stability . 

2.6.1 Wind-Up Turn 

In a wind-up turn, the aircraft is rolled and high 
load factor is commanded, resulting in a high pitch rate. As 
airspeed bleeds r.ff (which may occur even at maximum thrust), 
angle of attack is increased, and the aircraft stability 
decreases. Five points taken from a typical wind-up turn 
time history are described in Table 2.6-1, and the corre- 
sponding eigenvalues are given in Table 2.6-2. Of special 
interest is the Dutch roll damping, which decreases so that 
the Dutch roll mode becomes unstable as the wind-up turn 
progresses. These are not symmetric flight conditions, so 
it is expected that the Dutch roll eigenvector also con- 
tains angle-of-attack perturbations. 

In addition to lateral-longitudinal coupling and the 
general reduction in damping, the control effectiveness also 
decreases, as illustrated by the initial value of the trans- 
fer function, shown in Table 2.6-3. It is necessary for the 
pilot to use rudder as the roll control at high angles of 
attack, and this can cause sideslip perturbations, which can 
lead to further problems. 

These representative points from a wind-up turn 
demonstrate the deterioration of the stability and control 
of the example aircraft as it executes one form of air 
combat maneuver. 
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Working 

Point 


TABLE 2.6-1 

WIND-UP TURN WORKING POINTS 


Description 


Roll and Turn 
t * 0 sec 


Rapid Turn 
t * 13 sec 


Turning 

t * 30 sec 


Turning 
t * 52 sec 


Turning 

t * 75 sec 


Flight Condition 


TABLE 2.6-2 

WIND-UP TURN EIGENVALUES 


Working 

Point 


Short Period 

Dutch Roll 

Roll 

Spiral 

Phugoid 

-0 . 935? jl 76 

”0 . 383i j 3.43 

-1 61 

0.064 

-0.051* JO. 126 

-0. 952i J2 . 19 

-0 . 3u2 * j 4 .37 

-1 .04 

-0.043 

0.002* JO. 191 

- 0 . 777; j 1 .74 

-0.181tj4 .04 

-0.62C 

- 0 . 074 

-0 012 i jO . 209 

-0 . 456i j 2 . 66 

0 . 0766* j 3 26 

-0 7 58 

- 0 . 096 

-0,035ij0.327 

-0.438iJ2.16 

0.213+J2.55 

-0 543 

-0.154 

-0.057ij0.214 


















TABLE 2.6-3 


TRANSFER FUNCTION GAIN, Kj , ALONG THE WIND-UP TURN 


Working 

Point 

4w/AS h 

Ap/A« a 

Ap/a6 r 

Ar/A5 a 

Ar/A« r 

Av/A« a 

Av/A5 r 

1 

-0.42 

15.80 

3.80 

0.321 

-2.87 

- 0 . 01 S 

0.127 

2 

- 0.40 

12.12 

3.45 

0.095 

-2.91 

- 0.011 

0.131 

3 

-0.49 

5.91 

2.14 

-0.035 

-2.03 

- 0.002 

0 090 

) 

- 0.31 

1.32 

1.08 

0.037 

- 1.06 

0*002 

0 . 04 S 

5 

- 0.17 

-0.001 

0.24 

0.067 

-0. '6 

0.002 

C .029 


2.6.2 Rolling Reversal 

A rolling reversal combines a rapid pull-up with a 
rapid rolling maneuver, resulting in a "corkscrew-like" path 
through space. The combination of a high-acceleration pull- 
up and rapid rolling is expected to produce unstable modes 
with considerable lateral-longitudinal coupling. Table 

2.6- 4 describes the rolling reversal working points examined 
here. The corresponding eigenvalues, shown in Table 2.6-5, 
illustrate the changes in aircraft stability as the rolling 
reversal progresses. Due to the high involved in this 
maneuver, the Dutch roll mode is unstable throughout most of 
the maneuver. 

The initial and final working points of the rolling 
reversal are symmetric flight conditions so there is no 
lateral-longitudinal coupling during these phases of the 
flight. This is demonstrated by the eigenvectors of the fast 
modes at the first working point, which are shown in Fig. 

2.6- 1. The intermediate working points all occur during the 
aircraft's roll and involve significant 1 a t era l - longi tudi nal 
coupling. The eigenvectors of the fast modes at Working Point 
3, shown in Fig. 2.6-1, demonstrate this There is significant 




TABLE 2.6-4 

ROLLING REVERSAL WORKING POINTS 


Working 

Point 


Description 

Flight Condition 

Hi-G Pull-up 

V Q ■ 217 n./s u Q » 25' deg 

t * 0 sec 

* 15 deg/sec 


e Q ■ 30 deg 

Roll 

Vq * 163 m/s c*Q • 26 deg 

t * 4 sec 

p Q ■ -10 deg/sec q^ * 15 deg/sec 


$ 0 - -90 deg e Q » 50 deg 

Roll 

V « 101 m/s Oq “ 26 deg 

t * 10 sec 

p Q * -15 deg/sec * 15 deg/sec 


r Q * -5 deg/sec ■ -180 deg 


9 q * -30 deg 

Roll and Pull-up 

V Q - 117 m/s a Q - 23 deg 

t ■ 15 sec 

p Q * 10 deg/sec qg • 15 deg/sec 


r Q * 5 deg/sec * 45 deg 


00 “ -55 deg 

Final Pull-up 

V ■ 152 m/s a Q ■ 20 deg 

t * 22 sec 

p Q - 0 deg/sec ■ 10 deg/sec 


<J» 0 - 0 deg 0 O - -10 deg 


Working 

Point 


TABLE 2.6-5 

ROLLING REVERSAL EIGENVALUES 


Short Period 

Dutch Roll 

Roll 

Spiral 

Phugoid 

-0.873± J4 . 07 

0.191SJ4.54 

-0.910 

0.139 

-0.178 -0.056 

1 

-0.691sj3.15 

0 . 180± J3 . 55 

-0.716 

-C.169 

0. 0124t J0.406 

-0.403tjl.Sl 

0.1071J2.10 

-0.409 

-0.224 

0.122 ± JO. 124 

-0.394tJ2.34 

0 . 131± J2 . 70 

-0.711 

-0.093 

-0.153 i JO. 124 

-0.593±J1.78 

-0.010tJ3.83 

-0.818 

-0.082 

-0.073 i JO. 095 
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Figure 2.6-1 Eigenvectors of Rolling Reversal 


angie-of-attack (or Aw) motion in the Dutch roll mode and large 
roll-sideslip perturbation in the short period mode. Even the 
roll mode contains a significant angle-of-attack excursion. 


The aircraft control effectiveness follows trends 
similar to the aircraft stability, i.e., the control effec- 
tiveness is degraded throughout the first half of the maneu- 
ver, but it improves during the second half, as illustrated 
in Table 2.6-6. The pilot must use the rudder as a roll con- 
trol during the middle portion of this maneuver. This diffi- 
culty is complicated by high angular rates, lateral-longi- 
tudinal coupling, and extreme attitudes. 


An aircraft executing a rolling reversal exhibits 
unstable modes, lateral-longitudinal coupling, and reduced 
control effectiveness as the maneuver progresses, all of 
which make the pilot's task more difficult. 
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TABLE 2.6-6 


TRANSFER FUNCTION GAIN, Kj , ALONG THE ROLLING REVERSAL 


Working 

Point 

Aw/A6 h 

AP /&« a 

ap / a« r 

Ar/Afi a 

tr/&6 r 

Av/AA^ 

Av/A6 r 

1 

- 0.62 

1.87 

1.53 

0.204 

^ 2 ! 25 " 

0.004 

0.109 

2 

- 0.36 

0.597 

0.781 

0.142 

-1.29 

0.003 

0.063 

3 

- 0.012 

-0.001 

0.282 

0.053 

-0.465 

0.001 

0.023 

4 

- 0.23 

0.683 

0.722 

0.044 

-0.771 

0.001 

0.036 

5 

- 0.42 

2.91 

1.73 

-0.031 

-1.513 

0.002 

0.006 


2.6.3 Effects of Proportional Tracking 

Rudimentary piloting effects can be examined by 
assuming that the pilot attempts to control the aircraft's 
attitude. This can be modeled by a proportional feedback 
of angular deviation (pitch or roll angle) to the appropriate 
control surface (elevator or aileron). Pilot lags or time 
delays are neglected. As an example, pitch attitude control 
is chosen, and the feedback gains are set so that the effec- 
tive pitch moment due to pitch angle, Mg, is a multiple of 
the pitch moment due to angle of attack, M . Since this 
is achieved by elevator feedback, there are changes to the 
coefficients Xg and Zg as wel] . The multiplying factor is 
denoted by "i" in the following tables, and the feedback 
gain that produces equal Mg and M a (i = 1) for the reference 
flight condition is 0.64 deg elevator per deg of pitch 
angle. This simple model also disregards pitch-rate feed- 
back, which the pilot also might normally provide. 

For symmetric flight conditions, this loop closure 
does not affect the lateral eigenvalues; the longitudinal 
eigenvalue variations are shown in Table 2.6-7. The short- 
period mode is both increased in frequency and decreased in 



TABLE 2.6-7 


EIGENVALUE CHANGES DUE TO PROPORTIONAL TRACKING — 
SYMMETRIC FLIGHT CONDITIONS 


V Q ■ 94 m/s a Q * 15 

deg q Q * 0 deg/sec 

i 

Short Period 

Phugoid (Pitch 
Angle/Speed) 

0 

-0.407± jl . 099 

-0 . 017± jO . 137 

1 

-0. 342± jl . 552 

-0 . 082± jO . 068 

2 

-0. 319±jl . 906 

-0.157 -0.0522 

4 

-0 . 302± j2 . 47 

-0.218 -0.0260 


V Q * 94 m/s a Q = 15 

deg q Q = 12 deg/sec 

i 

Short Period 

Phugoid (Pitch 
Angle/Speed ) 

0 

-0 . 417 ± jl . 13 

0.045 -0.069 

1 

-0. 351tjl . 57 

-0.096 -0.058 

2 

-0.327±jl . 92 

-0.140 -0.063 

4 

-0 , 307 ±j2 . 48 

-0.178 -0.064 


damping by the addition of a pitch attitude-to-elevatcr 
feedback. For straight-and-level flight, attitud control 
increases the phugoid damping while decreasing th j.tural 
frequency, resulting in the conversion of the phugoid mode 
into two real modes — a pitching mode and a speed mode. 

In steady pitching motion, pitch-attitude control results in 
increased stability for the pitch angle mode with relatively 
little effect on the speed mode. 

Coupled flight conditions lead to significant 
effects on the lateral modes due to the longitudinal loop 
closure. In Ref. 9, a pitch attitude-to-elevator loop 



closure resulted in an unstable lateral mode when the sub- 
ject aircraft was in a steady sideslip. As shown in Table 
2.6-8 for the example aircraft used in this report, pitch 
angle-to-elevator feedback generally has a beneficial 
influence on the lateral eigenvalues when this aircraft is 
rolling at zero sideslip. However, the presence of a non- 
zero nominal sideslip angle results in a mildly diverging 
speed mode for moderate-to-large feedback gains. There is 
an indication that large attitude feedback gains destabilize 
the Dutch roll mode (for a non-rolling aircraft) or the short 
period mode (when the aircraft is rolling and slipping). 

Regarding the pitch attitude-to-elevator feedback 
as a simple pilot model, this examination confirms the earlier 
result that pilot control of the longitudinal motion of an 
aircraft could result in the destabilization of the aircraft 
when the vehicle is in a steady sideslip. This would be due 
to the pilot disregarding the lateral-longitudinal coupling 
present in asymmetric flight conditions. 

This simple attitude feedback underlines the 
necessity for considering lateral-longitudinal crossfeeds 
when designing a stability augmentation system for a high- 
performance aicraft . To achieve acceptable performance, 
it may be necessary to design a system that recognizes the 
aircraft flight condition and adjusts its gains to suit the 
situation. 


2 . 7 CHAPTER SUMMARY 

This chapter has presented a study of the dynamic 
charac eristics of a high-performance aircraft, with special 
emphasis on the effects of extreme flight conditions on 
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TABLE 2.6-8 

EIGENVALUE CHANGES DUE TO PROPORTIONAL TRACKING - 
SIDESLIP AND ROLL EFFECTS 


V Q = 94 m/s ct Q = 15 deg 8 Q = 0 deg p WQ = 


Short Period 

Dutch Roll 

Roll 

-0. 407± jl . 099 

0.0738±j2.25 

-0.442 

-0. 342± jl . 552 

0.07381 j2.25 

-0.442 

-0.319±j 1.906 

0. 07381 j2 . 25 

-0.442 

-0.302±j2.47 

0 . 07381 j2 . 25 

-0.442 


0 deg/sec 


Phugoid 


-0.0545 -0 . 017± jO . 137 

-0.0545 -0 . 082± jO . 068 

-0.0545 -0.157 -0.0522 

-0.0545 -0.218 -0.0260 


V Q = 94 m/s 


Short Period 


a. =15 deg g = 10 deg p w = 0 deg/ sec 


Dutch Roll Roll Spiral 


Phugoid 


-0.024±j0.146 
-0 . 032± jO . 124 
-0.287 |o . 133 
-0. 0247±j0.108 


= -39 deg/sec 


Phugoid 


0.303 


2.11 

2.11 

2.03 

2.08 



-0.0315 

0.081 

- 0.021 

-0.152 


V Q = 94 m/s a Q = 15 deg B Q = 10 deg p WQ 


Short Period Dutch Roll Roll Spiral 


-0 . 511± j 1 . 32 
-0 . 349± j 1 . 58 
-0.201±jl .77 
-0 . 038± jl . 93 


-0 . 057± j2 . 28 
-0 . 116± j 2 . 34 
-0 . 202± j2 . 44 
-0 . 306± j 2 . 80 


-0 . 246± jO . 199 
-0 . 157 ± jO . 413 
-0 . 215± jO . 480 
-0 . 270±j0 . 518 


0 . 071± jO . 207 
-0.193 -0.050 

-0.249 -0.0035 

-0.275 10.0155 


V n = 94 m/s a = 15 deg 


= 0 deg p w = 


Short Period Dutch Roll Roll Spiral 


-39 deg/sec 


Phugoid 


0 . 464± j 1 . 29 
0 . 400± j 1 . 56 
0.334ijl .82 
0.169^12.10 


-0 . 032 ±,j 2 . 34 
-0.044 t j2 . 37 
-0.073' J2.42 
-0.201 *-.) 2 . 67 


-0 . 293± jO . 144 
-0. 167± jO. 500 
-0 . 208+ jO . 600 

-0 . 246^.1 0 . 666 


0.043±j0.189 
-0.207 -0.062 

-0.229 -0.035 

-0.238 -0.0228 



































aircraft stability and control. The chapter first examines 
previous studies of the dynamics, aerodynamics, and control 
in this flight environment. The difficulty of measuring 
angular rate and translational acceleration effects leads 
to limited availability of this aerodynamic data, which 
has a significant impact on the simulation and analysis of 
dynamic departures. The survey of stability and control 
indicates a need for additional developments in these areas. 
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EFFECTS OF CONFIGURATION VARIATIONS 
ON AIRCRAFT DYNAMICS 


3.1 OVERVIEW 

Variations of aircraft configuration lead to changes 
in the aircraft's eigenvalues, eigenvectors, and control 
effectiveness. Section 3.2 presents the effects of changes 
in the most important longitudinal stability derivatives on 
the mode shapes and speeds. Similar effects caused by changes 
in lateral stability derivatives are detailed in Section 3.3 . 
The effects of aircraft mass and rotational inertia variations 
are given in Section 3,4 . Section 3.5 presents a general dis- 
cussion of possible departure modes and illustrates some of 
the possible departure time histories. Section 3.6 is a 
summary of the chapter. 


3.2 VARIATIONS DUE TO LONGITUDINAL STABILITY DERIVATIVES 

The longitudinal stability derivatives determine the 
aerodynamic force and moment contributions to the longitudinal 
perturbation equations, and these stability derivatives can 
vary considerably from aircraft to aircraft. This section 
surveys the changes in normal mode shapes and speeds for 
different ranges of the most important longitudinal stability 
derivatives . 

Three aerodynamic derivatives dominate the short 
period motion of the- aircraf t : C mc} , C ma , and C^. The sig- 

nificance of these terms can be seen in reduced-order approxi- 
mations to the damping ratio, ;, and natural frequency, u of 
this mode: 
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Over the range of likely values of the coefficients, the CmqCz^ 
product usually is considerably smaller than C ma , so it can 
be expected that C m and C% primarily affect short period 

4 

damping, while C ma changes only the short period natural fre- 
quency . 


The specific effects of varying Cz Q and C m q are 
illustrated in Table 3.2-1. The primary short period eigenvalue 
changes occur in the damping, as expected. The phugoid mode 
does change somewhat, as an increase in the lift-curve slope 
(Cz a more negative) increases the phugoid natural frequency 
at essentially constant damping ratio. An increase in the 
magnitude of pitch damping (Cm ) decreases phugoid frequency 
and damping significantly. 

TABLE 3.2-1 

Cz AND C m EFFECTS ON EIGENVALUES 
ns m q 


ii 

O 

> 

94 m/s , a Q = 

15 deg) 

c z a ( rad 1 ) 

Short Period 

Phugoid 

-2.0 

-0.345±jl .096 

- 0. 0146t jO. 128 

-3.2 

-0. 396+ jl . 099 

-0 . 0162±,j 0 . 135 

-5.0 

-0.471+jl .099 

-0 . 0181+ jO. 145 


Cm q (rad X ) 

Short Period 

Phugoid 

-34.4 

-17.2 

-5.73 

-0.6421 jl .058 
-0.399± jl .099 
-0.2361 jl. 098 

-0 . 013 + j 0 . x30 
-0 . 017i j 0 . 137 
-0 . 0201 jO . 143 


I 




Ac 




Variations in the longitudinal stability deriva- 
tives affect the lateral modes only when the aircraft is in 
asymmetric flight. In that case, mode coupling occurs — 
Table 3.2-2 illustrates the changes in the Dutch roll and 
short period modes for asymmetric flight as C mq varies. The 
transfer of damping from the short period mode to the Dutch 
roll mode for non-zero 6 q has been observed in Section 2 . 4 . 
In the case of reduced C mq , this transfer is essentially 
unchanged, indicating that short period damping due to Cz a 
(see Eq. (3.2-1)) is transferred whereas damping due to C mq 
is not. This is supported by the observation that Dutch 
roll eigenvectors for non-zero Bq (Fig. 2.4-4) include much 
more angle-of-attack motion than pitch rate. 


TABLE 3.2-2 

EFFECTS OF C mq ON EIGENVALUES IN ASYMMETRIC FLIGHT 


ct Q = 15 deg Bq 35 o deg p WQ = 0 deg/sec 


Short Period 

Dutch Roll 

C mq = -17.8 (rad -1 ) 
C mq = -5.73 (rad -1 ) 

-0 . 407 ± jl . 099 
-0. 236± jl . 098 

-0.074tj2.251 
-0 . 074± j 2 . 251 


a Q = 15 deg B Q = 5 deg p W() = 0 deg/sec 


Short Period 

Dutch Roll 

Cm q = -17.8 (rad -1 ) 
C mq = -5.73 (rad -1 ) 

-0 . 382+ j 1 . 201 
-0.205±j 1.205 

-0 0946tj2.024 
-0. 0956+ j2. 024 


a Q = 15 deg B Q = 

0 deg p WQ = -39 6 deg/sec 


Short Period 

Dutch Roll 

C m = -17.8 (rad" 1 ) 

q 

-0.464tjl.294 

-0 . 0319t j 2 . 342 

C m * -5.73 (rad -1 ) 

-0.291±j 1.295 

-0. 0343± j2 . 337 
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None of the stability derivatives discussed has any 
effect on the control effectiveness, which determines the 
transfer function gain, Kj ; therefore, the initial response 
to control inputs does not vary with Cz a and C m . However, 
stability derivative variations do affect the transient res- 
ponse through changes in poles and zeros. 

Variations in the aircraft's center of gravity cause 
variations in aerodynamic moment coefficients. The center 
of gravity (c.g.) is the rotational center of the aircraft. 

For fixed aerodynamic center of pressure, c.g. variation 
leads to static margin variation; hence, the moment rela- 
tionships are altered. 

The static margin is the distance between the c.g. 
location and the aerodynamic center, and it is usually expressed 
as a fraction of the mean aerodynamic chord, c. Figure 3.2-1 
details the changes in the longitudinal eigenvalues as the 
static margin is varied from 0.33 through its usual reference 
location of 0.17 to -0.15. This results in a C m variation 
from -1.17 to 0.554 rad , as well as changes in C,,^ from 
-21.1 rad -1 to -11.2 rad -1 . 

The eigenvectors (Fig. 3.2-2) change considerably 
as the c.g. moves aft. While the short period and phugoid 
modes still are recognizable at a static margin of 0.06 
(Cm a = -0.17 rad -1 ), a transition region is entered as the 
c.g. moves further aft. At a static margin of 0.01 
^C ma = -0.02 rad -1 ), a new ("third”) oscillatory mode which 
displays significant perturbations in all longitudinal states 
is evident. Two real convergences comprise the other longi- 
tudinal modes — one fast attitude mode occurring at constant 
velocity and flight path angle and one slow velocity mode 
that involves significant flight path angle variations. 
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Control effectiveness changes only slightly with 
c.g. location, but the transfer function gain, Kp, does 
vary as shown in Table 3.2-3. The zeros of the Aw/A$ h trans- 
fer function are given in Table 3.2-4. Both of these tables 
indicate the well-behaved nature of the numerator of the 
elevator-to-angle of attack transfer function. 

TABLE 3.2-3 

C.G. LOCATION EFFECTS ON TRANSFER 
FUNCTION GAIN, K F 

(V Q = 94 m/s, cxq = 15 deg) 


Static Margin 

Aw/ A 6^ 

Ap/A5 a 

Ar/A6 r 

0.17 

-5.17 

-0.66 

-5.39 

0.05 

-5.76 

-0.57 

-5.51 

CO 

o 

o 

1 

-6.15 

-0.51 

i 

-5.59 


TABLE 3.2-4 

COMPARISON OF ZEROS OF Aw/ A 6. AT 
THRFF C.G. LOCATIONS 


Static Margin 

Z 1 ,2 

Z 3 

0. 17 

-0. 0025+ j0. 1167 

-21 . 07 

0.05 

-0. 0038±.j0. 1238 

-20.998 

-0.03 

-0. 0047i jO. 1283 

-20.950 


This section examines the effects of longitudinal 
stability derivatives, and the following conclusions are 

made : 


67 




I 


I 


J. 


• Reductions in pitch damping (C mq ) and 

lift-curve slope (Cg^) reduce the short 

period damping without changing the 
frequency . 

• Lateral-longitudinal coupling produced 
by sideslip transfers damping due to 

to the Dutch roll mode from the shoit a 
period mode, but damping due to C mq re- 
mains in the longitudinal plane. 4 

• The c.g. location affects Cn^ directly, 
and a rearward c.g. location results in 
the creation of a new unstable oscillatory 
mode that exhibits significant perturba- 
tions in all longitudinal variables. In 
addition, two stable real mod^s are 
created — a fast attitude mode and a 
slow velocity mode. 

• The numerator of the elevator-to-anfle of 
attack transfer function is well behaved 
at rearward c g. locations. 


3.3 VARIATIONS DUE TO LATERAL-DIRECTIONAL STAEILITY 
DERIVATIVES 


In this section, variations in 


8’ Cn r 


c l. 


md 


C n are studied to determine their effects on the lateral- 
l *r 

directional modes. Experimental information on the value 

of C n £ is limited, so an investigation of possible effects 
of C n . on eigenvalues is included in this section. 

D 


The effects of C n£ , C^ g , and C lp variations are 
shown in Fig. 3.3-1. Increases in magnitude cl C n g and C^g 
cause an increase in the frequency uf the Dutch roll mode. 
Additionally, both parameters cause some change in damping, 
«?ith larger C n g increasing the damping ratio and larger 
Ci magnitude decreasing the Dutch roll damping ratio. 

A D 

Larger C t magnitude causes the piral mode to be more stable. 


) 


! 
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•••• C, VARIATION FROM 0.0 TO -0.5rod 1 
*P 

— C,^ VARIATION FROM - 0.095 TO -0.34 rod*’ 

- - C„ VARIATION FROM 0.0 TO 0.42 rod*’ 

P 

• REFERENCE VALUE 


j 0» t- » O H 9 

-L0 0 19 

Figure 3.3-1 Effects of Cig, C n g, and C]_ Variations 

on Lateral-Direction? 1 Eigenvalues 

Figure 3.3-1 indicates that at clq of 15 deg the 
effects of Ci appear in the Dutch roll damping ratio and 
in the roll convergence mode. The presence of these effects 
is predicted by the approximate lateral-directional equa- 
tions, Eqs. (2.4-1) to (2.4-5), but the approximate equa- 
tions give an inaccurate indication of tbe size of these 
effects. Equations (2.4-1) and (2.4-4) predict much larger 
variation in the roll inode eigenvalue due to Ci variation 
than is indicated in Fig. 3.4-1. Conversely, Eqs. (2.4-2) 
and (2.4-5) predict a much smaller variation in Dutch roll 
damping ratio than occurs in the complete model of the 
subject aircraft. 


jw !*•<*') 


ROU 



DUTCH 

ROLL 


o 0 «IS* 


SPIRAL 




) 

\ 


69 


v 


I 


i 


The lack of accuracy of the approximate equa- 
tions results from the fact that the subject aircraft does 
not conform to the assumptions upon which the approximate 
equations are based. The subject aircraft exhibits much 
more roll motion in Dutch roll mode (as shown by the 
eigenvectors in Fig. 2.4-3) because it rolls easily (due 
to low I /I ratio) and because its roll-yaw and roll- 
sideslip cross derivatives (Ci r> Cp^, and Cig) are not 
small when properly compared to the roll and yaw-sideslip 
derivatives (C x , C nr> and C ng ). 

Table 3.3-1 compares the lateral-directional eigen- 
values for C 0 g and C j_g variations in the presence of non- 
zero nominal pitch rate. A comparison of the eigenvalues 
indicates that the effects of' steady pitching are independent 
of the changes in aerodynamic coefficients. In all cases , 
positive pitch rate destabilizes the Dutch roll and spiral 
modes and speeds up the roll mode . 

TABLE 3.3-1 

EFFECTS OF C ng AND C lg VARIATIONS IN THE 
PRESENCE OF STEADY PITCH RATE 




Dutch Roll 

Roll 

Spiral 

Ref ereno ~ 
Values 

q 0 » 0 

q Q = 12 deg/sec 

-0 . 0738± j 2 . 251 
0 . 114 i j 2 . 320 

-0.443 

-0.839 

-0.0545 

0.0220 

C n .> = 0.0 

0 -1 
(rad ) 

q 0 = 0 

q Q = 12 deg/sec 

-0.0043±jl .855 
0 . 214 ±jl . 996 

-0.499 

-1.016 

-0.127 

0.0063 

Ci = -0 . 338 
^ -1 

<1 0 - 0 

-0.0505±jl .996 

-0.459 

-0.0752 

( rad ) 

q^ = 12 deg/sec 

0 . 150 i j 2 . 638 

-0.904 

0.0135 


( 

/ 

V. 
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The approximations of Eqs. (2.4-2) and ( 2 . 4 - 5 ) 
indicate that the change in the Dutch roll mode's real 
part should be approximately proportional to the change in 
C nr - Table 3.3-2 indicates that this is not true in this 
example, for as Cn r increases by a factor of 20, the Dutch 
roll mode's real part cnly doubles. The other significant 
effect caused by a veduction in yaw damping magnitude is 
to destabilize the spiral mode. 


TABLE 3.3-2 

EFFECTS OF C n VARIATIONS 
11 r 


Cn r _ 

Dutch Roll 

Roll 

Spiral 

-1.200 

-0 . 0983± j2 . 2468 

-0.4612 

-0.1396 

-0.646 

-0 . 0738± j2 . 2508 

-0.4425 

-0.0545 

-0.0602 

-0.0477±j2. 2523 

-0.4304 

0.0422 


The derivative C n s is an acceleration derivative 

P 

(analogous to C m ^) that arises because the aerodynamic flow 
field exhibits some lag in rearranging itself following a 
change in aerodynamic angle. Acceleration terms are approxi- 
mations to these flow field dynamics. Experimentally, 
acceleration derivatives are difficult to measure, and 
they usually are combined with the rotary derivatives. 
Analytical studies often make the assumptions that 


c ng 

= 0 




A 

% 

= Cn P 

+ sin 

a 0 

c ng 

^nr 

— n 

~n r 

- cos 

a o 

C n& 


( 3 . 3 - 7 ) 

( 3 . 3 - 8 ) 

( 3 . 3 - 9 ) 
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Reference 46 examines some of the possible effects 
of this approximation on aircraft time histories and param- 
eter identification. It is concluded that significant param- 
eter identification "errors can occur if the acceleration 
derivatives are not included in the model when necessary. 

Two types of C n g effects have been examined. In 
the first type, non-zero C n s is added to the model with no 

p 

other change; this effectively decreases the total damping. 

In the second type, C 0p and Cn r are adjusted by subtract- 
ing the C n £ terms given in Eqs. (3.3-8) and (3.3-9) to 
maintain nearly constant damping. In the variable damping 
case (in which C nr and C np are constant), the only change is 
in the Dutch roll eigenvalues. There is significant change 
in damping and a slight change in frequency. This indicates 
that Cns primarily affects the Dutch roll mode. The constant 
damping cases are quite different, in that the changes in 
C n g, C np , and Cn r effectively cancel, as far as the Dutch 
roll mode is concerned. In this case, however, there are 
significant changes in the roll and spiral modes. These 
changes are not due to C n s (which does not affect the roll 

p 

and spiral modes) but rather to the corresponding changes in 
Cn r and Cn p . Therefoic, combining C n ^ with C 0p and C nr leads 
to erroneous eigenvalue calculations in modes not directly 
affected by C n g at all. 

Table 3.3-3 illustrates the constant darning results 
with and without steady rolling motion. Non-zero C n j affects 

P 

the frequencies of the oscillations only slightly in the 
presence of rolling motion, but there are significant changes 
in damping of the various modes. Dutch roll damping ratio 
changes for non-zero Cnj when the vehicle is rolling, and 



\ 
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TABLE 3.3-3 

EFFECTS OF C n g IN THE PRESENCE OF STEADY ROLLING 

(C n and C n Adjusted to Maintain Constant Damping) 
P I* 


T-0251 



Short Period 

Dutch Roll 

Roll Spiral 

Phugoid 

c ng 

= 0 (rad -1 ) 

P\V 0 * 0 deg/ sec 




-0.407ijl.099 

-0.073Sij2.251 

-0.443 -0.0545 

-0 . 0165i jO . 137 

c ng 

* 0 (rad -1 ) 

PWo = ~ s 9.6 deg/se 

c 



-0.464ijl.294 

-0 . 0319i j2 . 342 

-0 . 293i j0 . 144 

0. 04271 jO. 189 

c ng 

= 0.55 (rad^ 1 

) Pw 0 * 0 deg/ sec 




-0.407ljl.099 

-0 . 0757i j2 . 256 

-0.525 0.0321 

-0.01651 JO. 137 

c n§ 

= 0.55 (rad 1 

) Pwo = “39.6 deg 

/sec 



-0.48H jl.290 

-0 . 04111 j 2 . 345 

-0.2S7ij0.124 

-0.06301 jO.183 


i 

\ 

V 


this is not observed for the non-rolling, constant damping 
case. The conclusion is that an improper C n jj identification 
may not appear in the Dutch roll mode for non-rolling flight, 
but it San cause a significant change in the Dutch roll mode 
during a rolling maneuver . 


The following points summarize the findings of 
this report concerning the effects of lateral stability 
derivative variations: 






The large amount of rolling in the 

subject aircraft means that Ci 

••-Tj.dyn 

and C_ ?.re poor indicators of 

n r , dyn 

the roll mode's eigenvalue and the 
Dutch roll mode's real part. 


C n g and C^ variations primarily affect 

the Dutch roll frequency, as indicated 

bv C . Variations due to steady 

ng.dyn 
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pitching motion are independent of the 
effects caused by changes in C n g and C^g. 

• Due to the large amount of rolling motion 

in the Dutch roll mode, has a larger 

Ir 

than expected effect on Dutch roll damp- 
ing,- while the effect of C n is smaller 
than might be expected. 

• C n £ by itself primarily affects the Dutch 

roll mode, but subtracting its effect 
from C np and C Dr (holding the damping 

constant ) leads to large roll mode and 
spiral mode changes and small Dutch roll 
mode variations. 

• In the presence of steady rolling, there is 
a significant C n * effect on the Dutch roll 

P 

mode even for the constant damping case. 
Coupled flight conditions should be 
investigrted when the identification 
of C n g is desired. 


3.4 VARIATIONS DUE TO MASS AND INERTIA EFFECTS 

The variations in the aircraft modes due to mass and 
rotational inertia changes are examined in this chapter. The 
rotational inertias considered span the range of aircraft 
types (from wing-heavy to fuselage-heavy), and the mass 
variations range from light to heavy wing loadings. 

Table 3.4-1 details the aircraft eigenvalue trends 
as aircraft mass varies. The rotational inertia matrix 
is held constant (as if a point mass was added or subtracted 
at the vehicle’s center of gravity), so the ratios between 
the mass and inertia also vary. The mass change represents 
a change in relative density, u, which is defined as 

y = m/pSc (3.4-1) 
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TABLE 3.4-1 

EFFECTS OF AIRCRAFT MASS ON EIGENVALUE LOCATION 
(V Q = 94 m/s, a Q = 15 deg) 


Mass 

Short Period 

Phugoid 

Dutch Roll 

Roll 

Spiral 

3060. kg 

-0.572ijl.07 

-0 . 028± jO . 186 

-0.120±j2.25 

-0.442 

-0.054 

5055. kg 

-0.442ijl.09 

-0 . 019i jO . 150 

-0.084ij2.25 

-0.442 

-0.055 

8175 kg 

-0 . 367 i j 1 . 10 

-0.013±j0.120 

-0. 062i j2 . 25 

-0.443 

-0.055 


whore c is 3 reference length. The relative density relates 
the aircraft mass to the air density, and therefore indicates 
the relative magnitude of aerodynamic and inertial effects. 

Variations of the aircraft mass have a large effect 
on the damping of the rotational oscillations — Dutch roll 
mode and short period mode — but only a negligible effect 
on their frequencies. The spiral and roll modes also are 
essentially unchanged. Mass variations change the phugoid 
mode's natural frequency and damping ratio, because this 
mode involves the interchange of kinetic and potential 
energy and is highly mass dependent. 

The rotational inertia of the aircraft describes the 
distribution of the aircraft mass about the center of gravity 
Most of the inertia is due to the fuselage and the wing, and 
the relation between them leads to the designation of a 
specific configuration as "wing-heavy' or "fuselage-heavy." 

Th<- yaw inertia is approximately 10 to 15 percent 
larger than t he pitch inertia, and the roll inertia can be 
from 4 to 12 times smaller than the yaw inertia. High per- 
formance fighters emphasize high rolling performance, and 
tend to be f u.selage-heavv . Transport aircraft are built for 
cruising efficiency and tend to be wing-heavy. 
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Table 3.4-2 details the results of two rotational 

inertia investigations, one with constant rolling inertia 

(I ) and varying fuselage inertia (I and I), and one with 
x y z 

constant fuselage inertia and varying rolling inertia. In 
neither of these cas* is the phugoid mode significantly 
affected, underlining the conclusion that the phugoid is a 
translational mode rather than a rotational mode. For the 
uncoupled reference flight condition, the short period mode is 
affected only by a change in pitch inertia; an increase in 
pitch inertia causes a major decrease in short period fre- 
quency and a small decrease in damping ratio. 

TABLE 3.4-2 

EFFECT OF ROTATIONAL INERTIA ON EIGENVALUES 


Fuselage Inertia Varied — I x Held Constant (Vq * 81 n/s , a Q - 25 deg) 


V X x 

y 1 * 

Short Period 

Dutch Roll 

Roll 

Spiral 

— 

Phugoid 

4.4 

4.0 

-0. 534* J 2. 01 

0.102t jl.59 

-0.483 

-0.071 

-0.018tj0.156 

8.8 

8.0 

-0.329±jl.43 

0. 085*J1 . 54 

-0.366 

-0.050 

-0. 018ij0.157 

13.2 

12.0 

-0. 261*J1 .17 

0. 0788* j 1.52 

-0.327 

-0.039 

-0. 017 =j 0.157 


Roll Inertia Varied - I y 

and I z Held Constant (V Q « 94 m/s, a Q * 15 deg) j 

V 1 * 

y 1 * 

Short Period 

Dutch Roll 

Roll 

Spiral 

Phugoid 

4.4 

4.0 

-0. 407* jl . 099 

-0. 0428*jl . 826 

-0.313 

-0.0537 

-Q.0165*j0.137 

8.8 

8.0 

-0.407±jl . 099 

-0 . 0738±J2 . 251 

-0.443 

-0.0545 

-0. 0165s JO. 137 


The lateral modes are affected by an increase in 
rolling inertia, in that both Dutch roll frequency and roll 
mode response are slowed significantly. There is a signi- 
ficant decrease in Dutch roll damping and very little change 
in the spiral mode. Larger values of yaw inertia lead to 
somewhat different effects. Both roll and spiral modes are 
significantly slower, and there is some decrease in Dutch 
roll damping ratio and frequency. 
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The effects of a large rolling inertia on the 
rolling/slipping stability are shown in Fig. 3.4-1. The 
stability boundaries are generally similar to those shown 
in Fig. 2.4-10, although they differ in detail. The Dutch 
roll mode is stable for the higher rolling inertia. The 
phugoid instability combines with an unstable spiral mode 
at high sideslip angles and roll rates to form a fast, 
highly unstable oscillation. 


SlOESUP. B 0 ld««) 



m pnuoao 

V /\ STAKE 


F.gure 3.4-1 Effects of Large Rolling Inertia 

on Aircraft Stability 
(I x = 14,370 kg-m2) 


Conclusions concerning the effects of mass and 
inertia variations are summarized as follows: 


Mass increases reduce Dutch roll and 
short period damping, but mass varia- 
tions do net have large effect on the 
short period or Datch roll mode fre- 
quencies. They do not cause signifi- 
cant changes in the roll or spiral 
modes . 
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• The phugoid mode is a translational mode 
and is greatly affected by mass varia- 
tions. The phugoid eigenvalue does not 
depend strongly on rotational inertia. 

• Increases in pitch inertia reduce short 
pei iod frequency and damping. 

• Increases in roll inertia increase the 
roll mode time constant, decrease the 
Dutch roll frequency significantly, and 
modify the effects of mode coupling due 
to asymmetric flight. 

• Increases in yaw inertia primarily slow the 
roll and spiral modes, and there is some 
effect on Dutch roll damping and frequency. 


/ 
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3.5 CLASSIFICATION OF DEPARTURES 

Departure from controlled flight can occur in two 
ways. Unforced departures are due to instabilities in the 
basic aircraft. Even if the pilot does not move the controls, 
small perturbations in the aircraft states build up until the 
aircraft can no longer be controlled. In a forced departure , 
the basic aircraft may or may not be unstable, but the addi- 
tion of a pilot loop closure creates an unstable vehicle- 
pilot system. The two following sections discuss these de- 
parture classes. 

3.5.1 Unforced Departure Modes 

Unforced departures occur when the pilot cannot or 
does not stabilize an unstable vehicle. The vehicle eigen- 
values directly indicate the open-loop system stability in 
this case, so that many of the stability boundaries that have 
been shown in this report can be classed as u nforced departure 
boundaries . 


k 


i 


! 

I 

S. 


78 



The speed of the fast modes (Dutch roll mode, short 
period mode and roll mode) is such that the stability of 
these modes is critical, and the least stable of these usually 
is the Dutch roll mode. This mode can become unstable in two 
ways, either resulting from negative "spring terms" or result- 
ing from negative damping, both of which can be influenced by 
aerodynamics and coupling effects. 

The approximate equations for the Dutch roll mode (Eqs. 
(2.4-1) and (2.4-2)) indicate a pure static instability for large, 
negative C n g, but the exact result is somewhat more complex. As 
in the case of longitudinal mode coupling due to positive Cm a 
(Section 3.2), directional instability can cause the Dutch roll 
mode to couple with the classical roll and spiral modes, and it 
can lead to a new oscillatory mode, analogous ro the so-called 
"roll-spiral" or "lateral-phugoid" mode. Figure 3.5-1 illus- 
trates a case in which negative C D g causes an oscillatory 
mode that has low natural frequency and is highly unstable. 

By comparison to the conventional Dutch roll mode (Fig. 

2.4-3), there is a significant change in mode shape. There 


SPIRAL ROLL-SPIRAL ROLL 


Chp * * 0 . 85 ? 
C. A ^*-0J44 
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Figure 3.5-1 An Example of Lateral-Directional Eigenvalues 

for Negative Directional Stability 
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is a 180-deg phase change in the yaw-rate component, as well 
as a substantial roll angle change. In addition, the spiral 
mode has gained significant Av, Ar, and Ap components. 

The departure caused by this type of instability is 
shown in Fig. 3.5-2. Although the linear model indicates 
that this motion is an oscillation, it. is so unstable that 
only part of a period appears on the . ae history plot. The 
first few seconds of the motion exhibit a rapid roll-yaw 
angular motion. The pilot would sense a rapid rotational 
divergence about this axis and might refer to it as a 



time. t|it«) Time, t (••«) 


Figure 3.5-2 An Unforced Departure Due 

to Negative C n ^ 
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Negative yaw damping leads to a more conventional 
destabilizing of the putch roll mode, which retains its 
characteristic mode shape. The time history of a departure 
due *o dynamic Dutch roll instability is shorn in Fig. 3.5-3, 
and the difference in shape from the departure due to static 
instability is apparent. 



Figure 3.5-3 An Unforced Departure Due to 

Negative Dutch Roll Da-nping 

The large amount of rolling motion in the Dutch roll 
mode indicates that this may be what pilots refer to as 
"wing rock." This is uncertain, however, since "wing rock" 
also could be a roll-spiral oscillation or a limit cycle caused 
by an aerodynamic non-linearity. In any case, a pilot sensing 
such an oscillation probably would unload t'e aircraft by 
reducing the angle of attack, removing the aircraft from the 
region of instability. 
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A third type of unforced departure can occur at high 
sideslip angles. The modes of motion ai«. shown in Fig. 3.5-4, 
a». ' the unstable roll mode is seen to exhibit a mixed rolling- 
yawing departure characteristic. This divergence is rapid 
and, as shown in the stability boundari 3 of Fig. 2.4-2, can 
appear with only small aerodynamic angle changes from a much 
more benign flight condition. 
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Figure 3.5-4 Eigenvalues and Eigenvectors for a Flight 
Condition with Large Sideslip Angle 


This unstable roll mode is essentially a pure roll 
about the stability x-axis, but, because of the large nominal 
aerodynamic angles, it appears as a roll-yaw motion in body axes. 
It is possible that a pilot flying an aircraft at these large 
angles would interpret a stability- axis rolling depai ture as 
a "nose slice." 


3.5.2 Forced Departure Modes 

Control inputs fiom a pilot or control system can 
force an aircraft to depart from controlled flight in two 
ways. In the first way, .Jverse response to pilot inputs 
moves the nominal flight condition into an unstable region 
where an unforced departure can or;cu r . The second possible 
cause of a forced departure is an improper loop closure that 
creates an unstable closed-loop system. Departure preven- 
tion procedures are quite different for the two cases; in 


82 



the former case, positive control action is necessary for 
recovery, whi'e a neutralization of control inputs might 
allow a recovery from the latter departure. 

An example of an improper loop closure was pre- 
sented in Section 2.6. In that case, the target-tracking 
pitc: attitude-to-elevator loop closure caused lateral mode 
stability problems at coupled flight conditions. This un- 
derlines the necessity of including control cross-couplings 
in situations where the system itself is coupled . 

The target -tracking example can be considered as 
a situation in which the pilot learns to control the air- 
craft at one flight condition but does not change his con- 
trol strategy as the flight condition changes. This is 
emphasized by the observation that coupled flight conditions 
often exhibit drastic changes in the shapes of the normal 
modes, so that a pilot might apply the wrong control action. 

The investigation of Section 2.5 demonstrates that 
control effectiveness problems may leave the pilot no alter- 
native but to apply a poor control combination. For example, 
a. 25-deg angle of attack, the roll moment due to aileron is 
essentially zero. The rudder would nave to be used for roll 
control but this brings an unavoidable sideslip response 
with it. This sideslip could drive the vehicle into the 
roll divergence region cited in the last section. 

As an example of unexpected control response, Fig. 
3.5-5 shows a departure caused by an aileron input. Nor- 
mally, the result would be a significant negative roll rate, 
but the sideslip and yaw rate build up so rapidly that the 
instability of the basic aircra r t results in a rapid rolling 
departure with positive roll rate. 
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Figure 3.5-5 Aileron Input for Negative C n ^ 

Situations in which improper pilot inputs are likely 
are discussed in Section 2.5. These situations are character- 
ized by reduced stability of the open-loop system (due to high 
angle of attack or pitch rate) and highly coupled modes caused 
by an asymmetric flight condition (such as non-zero 3 q orp^). 
Nonminimum-phase zeros often appear and can cause great dif- 
ficulty if "tight" control is attempted. 


The conclusions regarding departure modes are summar- 


ized as follows: 


I 





An unforced departure (one due to an ; 

unstable open-loop system) is most 

likely to appear in the lateral oscilla- / 

tory mode. Static instability results 

in a rapid rolling-yawing departure, 

while dynamic instability causes an 

unstable "wing rock" motion. 

An unstable wind-axis roll divergence 
can appear at extreme aerodynamic angles. 

Forced departures can occur when de- 
graded control response causes the pilot 
to fly the aircraft into a flight con- 
dition where unforced departures are 
likely. 

Mode coupling or unexpected nonminimum- 
phase zeros can change the control re- 
sponse so that a "normal" control loop 
closure leads to an unstable closed-loop 
system. 




3 . 6 CHAPTER SUMMARY 



This chapter has presented effects of configurational 
variations on aircraft dynamics. Relationships between mode 
approximations and exact results are discussed for longi- 
tudinal, lateral, and coupled motions, and examples of 
various departure types are presented. It is shown that 
the effects of aerodynamic parameter variations are modified 
by the coupling which results in asymmetric flight, par- 
ticularly in regard to the transfer of damping (due to 
rotary derivatives) from longitudinal to lateral-direc- 
tional modes (and vice-versa). Time histories of 1 inear ized- 
model response illustrate departure characteristics similar 
to those experienced in flight. 
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4. 


PREVENTION OF DEPARTURE FROM CONTROLLED FLIGHT 


4.1 OVERVIEW 

As indicated by earlier developments in this report 
and the summary of prior work in Section 2.2, there is ample 
reason to consider designing stability augmentation systems 
for the specific purpose of preventing departure. Aircraft 
design is dominated by performance requirements, and even 
unconstrained configuration modifications may not provide 
adequate stability or control response (especially during 
extreme maneuvering). Appendix A and Chapters 2 and 3 
demonstrate how linear-time-invariant models of aircraft 
dynamics can be derived for studying stability and control 
response during difficult maneuvers. These models are used 
to illustrate stability augmentation system concepts in the 
present chapter. 

Unlike earlier studies of departure prevention, 
the powerful tools of linear-optimal control theory are 
applied to the problem in this chapter. Since new ground is 
broken and methods which are unfamiliar (in the departure 
prevention context) are presented, the objective is to pro- 
vide preliminary guidelines for Departure-Prevention Sta- 
bility Augmentation System (DPSAS) development. Therefore, 
a simple optimal controller — the continuous-time linear - 
optimal regulator — is applied to departure prevention. A 
linear-optimal regulator is a feedback control law of the 
form. 

Au( t ) = -KAx(t) (4.1-1) 
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where Au(t) is the vector of control command perturbations, 
Ax(t) represents the vector of the aircraft’s dynamic states, 
and K is the gain matrix which scales the state measurements 
for proper stabilization ar.d compensation of -the aircraft's 
motion. (An equivalent discrete- time linear-optimal regu- 
lator, for which the state is measured and control is com- 
manded at a fixed sampling interval, can be derived for a 
digital flight control system.) This control law has several 
qualities which are desirable for the present study, in which 
i 1 is assumed that system dynamics are known exactly and that 
all states are measured precisely: 



• The control gains guarantee stability 
of the closed-loop system. 

• Complete longitudinal/lateral-directional 
coupling is assumed and is accounted for 
in the design process. 

o 

• The control design technique identifies 
all significant crossfeeds and intercon- 
nects, as well as feedback gains. 

• Tradeoffs between the amplitudes of state 
perturbations and of control motions are 
specified in the design process. 

In addition, a gain-scheduling algorithm which accounts for 
varying maneuver conditions is developed. 

The control design techniques applied to the DPSAS 
can be generalized to full command augmentation systems for 
a high-performance aircraft. Reference 58 shows how 
practical command-response control laws can be developed for 
a highly coupled aircraft, a tandem-rotor helicopter. These 
control laws satisfy classical step-response criteria, adapt 
to flight condition, honor rate- and displacement-limits on 
control actuators, and use incomplete (possibly noisy) 
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feedback measurements. A command-response system for a high 
performance fighter is described in Ref. 59. It adapts to 
flight condition to provide uniform handling qualities 
throughout the flight regime. These control laws are de- 
veloped for direct implementation in a digital computer 
and use low sampling rates. This extension of the DPSAS to, 
a complete flight control system, while promising, is a 
subject for future study. 

o 

The remainder of this chapter is directed to a 
brief explanation of linear -optimal regulator design and 
extensive application of this control design approach to 
DPSAS examples. Section 4.2 presents the linear-optimal reg- 
ulator and a discussion of the parameters used in computing 
control gains. Control designs are developed for a refer- 
ence aircraft over a wide range of angles of attack, pitch 
rates, sideslip angles, and roll rates at a single altitude- 
velocity point — the central flight condition of 6100 m, 

94 m/s — in Section 4.3 . The symmetric and asymmetric 
variations in flight condition are considered separately, 
in order to make the differentiation between pullup and 
sideslip-rolling effects more apparent. Control gains are 
computed at 32 maneuvering conditions to obtain the results 
of Section 4.3; with eight states fed back to four control 
effectors, over 1000 gains are generated. In Section 4.4 , 
these gains are correlated with each other and with maneuver 
conditions to identify candidate interconnects and gain- 
scheduling relationships. Negligible and constant gains also 
are identified in the process. The chapter is summarized in 
Section 4.5. 
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4.2 


THE LINEAR-OPTIMAL REGULATOR 


Optimal control theory provides a useful and 
practical multi-input, multi-output control system design 
tool. Linear-optimal control methods are based on the 
differential equations that describe the vehicle in the 
time domain (Eq. (A. 3-3)), and they produce feedback con- 
trollers which exhibit desirable properties. 

The problem is to find a controller for the system 
described by Eq. (A. 3-3), which exhibits a linear feed- 
back structure (Eq. (4,1-1)) and minimizes a scalar-valued 
cost functional of the state and the control : 

00 

J - \ (Ax T QAx + Au T RAu) dt (4.2-1) 

•'o 

This controller is called a linear-optimal regulator , and 
it is derived in Refs. 60 to 62. 

The designer's freedom rests in his choice of the 
weighting matrices, Q and R. The design procedure consists 
of the choice of Q and R, the computation of the Riccati 
matrix, an evaluation of closed-loop performance, and the 
adjustment of Q and R as discussed in Section A. 4. 4. 

The linear-optimal regulator is a tool for design- 
ing a Departure-Prevention Stability Augmentation System 
(DPSAS). It is not a limiter, because no limits are placed 
on the pilot's control authority, and it is not an auto- 
matic spin-recovery system, because open-loop anti-spin con- 
trol settings are not implemented. The DPSAS is intended to 
augment stability and to minimize the gyrations which pre- 
cede loss of pilot control. The DPSAS makes full use of 
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available control power, and, in this respect, could com- 
pete with the pilot's control commands; however, the de- 
signer can specify the amount of control-surface displace- 
ment which normally is available to the DPSAS. Basing the 
system on the linear-optimal regulator, the OPSAS can be 
designed to use less than full control authority for expected 
magnitudes of aircraft maneuvers, leaving a percentage of 
control authority free for manual commands. 


The primary objective of this chapter is to iden- 
tify the basic effects of varying flight condition on the 
structure of a DPSAS. To keep the number of varying param- 
eters to a minimum in this demonstration, the state and 
control weighting factors (Eq, (A. 4-21) and (A. 4-22)) are 
chosen at a single symmetric flight condition and held con- 
stant throughout the sweep of 32 maneuvering conditions. Q 
and R elements which provide satisfactory eigenvalues, 
acceptable time response, and reasonable control gains are 
chosen at the central flight condition of this sweep 
(cIq » 15 deg, V * 94 m/s fps, H * 6100 m) . Thus, it is expected 
that eigenvalues and control gains will vary with flight con- 
dition, but the rms-values of state and control perturbations 
should remain relatively constant. 


The Q and'R elements are used as design parameters 
which can be interpreted as the following maximum allowable 
rms perturbations: 


Throttle setting: 
Elevator deflection: 
Aileron deflection: 
Rudder deflection: 
Euler angle: 

Body angular rate: 
Body velocity: 


100% of full scale 

20 deg 

60 deg 

30 deg 

30 deg 

25 deg/sec 

9 m/s 


90 



’W *' « '- J -t *. 'r* - -*•*- *gs^ fi y 'i 


I 




These values indicate that throttle setting, ele- 
vator, aileron, and rudder are allowed to vary between their 
limits and that angles of attack and sideslip must be held 
within 5.6 deg (the 9 m/s body velocities correspond to 
aerodynamic angles of this number). Table 4.2-1 indicates 
that the primary effects of the loop closures at the cen- 
tral flight condition are to increase short period, Dutch 
roll, and phugoid damping and to quicken the roll and spiral 
modes . 


TABLE 4.2 -1 

EFFECTS OF DPSAS AT THE CENTRAL FLIGHT CONDITION 



Open-Loop Characteristic* 

Closed-Loop Characteristic* 

Dynamic 

Mode 

Natural 
Frequency , 
rad/*ec 

Damping 

Ratio. 

Time 

Constant , 
sec 

Natural 
Frequency , 
rad/scc 

Damping 

Ratio, 

Time 

Constant , 
sec 

Short Period 

1.17 

0.35 

- 

2.67 

0.72 

- 

Dutch Roll 

2.25 

0.03 

- 

2.35 

0.69 

- 

Roll 

- 

- 

2.26 

- 

- 

0.36 

Spin) 

- 


18.34 

- 

- 

1.07 

Phuyoid 

0.14 

0.12 

• 

0.16 

0.99 

- 



i 


This illustrates implicitly that the linear-optimal 
regulator design can produce stricter tracking than indicated i 

by the choice of Q and R elements. The Aa and AS require- j 

ments can be met only by increasing damping and decreasing j 

time constants. This infers that Euler angles and body 
angular rates also are closely regulated, even though the 
weighting of the corresponding elements in Q is light. Table 
4.2-1 also indicates that the selection of equal weights on 
Av and Aw (and, therefore, on Aa and AS) drives the natural 
frequencies and damping ratios of the short period and 
Dutch roll to similar values. 

) 
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The DPSAS gain matrix for this flight condition is 
listed in Table 4.2-2. The gain matrix illustrates why damp- 
ing is increased in the closed-loop system; rate feedbacks 
are large. The classical longitudinal/lateral-directional 
partition can be observed in the gains; the control algo- 
rithm actually computes coupling gains on the order of 10~7 
due to the use of single-precision arithmetic. These gains 
can be ignored. The elevator is seen to be the primary 
longitudinal controller, as throttle feedback gains are small 
(the principal effect of throttle control is to damp the 
phugoid mode). Lateral-directional control largely parti- 
tions along the roll and yaw axes. Although the gains shown 
in Table 4.2-2 have reasonable magnitudes, they could be 
reduced by reducing the values of (Eq. (A.4-21)). Tran- 
sient response would be altered, but the system would remain 
stable . 


TABLE 4.2-2 

DPSAS GAIN MATRIX AT THE CENTRAL 
FLIGHT CONDITION 


Control Output 

Pitch 

Anglo, 

dog 

Body x-axli 
Velocity, 
D/a 

Pitch 
Rato , 
dog/aoc 

Body z-axla 
Velocity , 

D/9 

Body y-axis 
Velocity, 
m/a 

Taw 
Rato, 
dog /aec 

Roll 

Rato, 

dog/aoc 

Roll 

Anglo, 

dog 

Throttle Sotting, 
Fraction of Full 
Seal* 

-0.019 

0.069 

-0.005 

0.020 

0 

0 

0 

0 

Elevator Anglo, 
dog 

-0.700 

1.332 

-1.540 

-0.988 

0 

0 

0 

0 

Aileron Anglo, 

dog 

0 

0 

0 

0 

-0.509 

0.138 

2.540 

1.894 

Bud dor Anglo, 
dog 

0 

0 

0 

0 

1.942 

-3.096 

0.225 

-0.100 


The performance of the linear-optimal regulator is 
demonstrated by comparing open- and closed-loop response to 
perturbations in angle of attack, sideslip angle, and roll rate. 
Figure 4.2-1 illustrates that a 1.1-deg Aa perturbation is 
moderately damped without the regulator ai.J well-damped with 
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the regulator. Figure 4.2-2 shows that the lightly damped J 

natural motion resulting from a 1-deg A 6 initial condition 

creates a substantial amount of roll as well a$ yaw. The ') 

regulator damps the oscillation and limits the roll angle 

excursion to 20 percent of its open-loop value, providing 

significant decoupling of lateral and directional motions. 

This decoupling effect is confirmed by the aircraft's closed- 
loop response to a roll-rate disturbance of 1 deg/sec (Fig- 
ure 4.2-3). This initial condition creates a small side- 
slip oscillation and triggers the spiral mode (indicated by 
the underlying exponential response trend in roll angle). 

The regulator damps the oscillation, reduces the sideslip 
response by 70 percent, and stabilizes the roll angle. 

Having obtained a representative design point for 
the DPSAS at the central flight condition, the effects of 
maneuvering on control gains, aircraft stability, and time 
response are examined in the next section. 

4.3 DPSAS CONTROL LAWS 

The control gains obtained at the central flight 

condition would stabilize the aircraft for some range of * 

v 

nominal angles and angular rates; however, changes in the 
aircraft's dynamics (reflected by variations in F and G) 
would lead to less-than-optimal regulation. It is neces- 
sary, therefore, to redesign the control gain matrix at 
each maneuvering condition in order to assess the full 
possibilities for preventing departure with the linear- 
optimal control law. 

Two separate maneuvering condition sweeps have been 
conducted, using the reference aircraft flying at 6100 m 
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Figure 4.2-2 Directional Response at the Central Flight Condition 
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and 94 m/s in both cases. The first is a longitudinal 
sweep , in which a range of angles of attack and pitch rate 
are considered. As indicated in Chapters 3 and 4, there is 
a significant change in lateral-directional dynamics during 
pullup maneuvers, although the longitudinal and lateral- 
directional axes remain uncoupled. The lateral-directional 
sweep varies sideslip angle and stability-axis roll rate, 
introducing full coupling about all three axes. In the 
first sweep, control gains and closed-loop characteristics 
change, but the DPSAS structure is conventional, i.e., gains 
are partitioned along usual lines. The second sweep generates 
unconventional DPSAS structures as well as gain variations. 

(In both cases, the control law is described by Eq. (4.1-1). 

K contains zero sub-matrices in the first sweep but not in 
the second. ) 

4.3.1 Longitudinal Sweep 

This section presents the effects of angle of attack 
and pitch rate on closed-loop eigenvalues, DPSAS control gains, 
and aircraft response. Section 2.4 showed that the reference 
aircraft has an unstable Dutch roll mode at high a Q and 
unstable Dutch roll, roll, and spiral modes at high qQ. These 
conditions are stabilized by the DPSAS. Using the state 
and control weighting factors discussed in the previous sec- 
tion, linear-optimal regulators are designed for 15 maneuver 
conditions (oiq varies from 5 to 25 deg, in 5-deg increments , 
and qQ is 0, 12, and 24 deg/sec). This sweep represents 
relatively low load factors (n =0.4 to 1.4 "g's"), it covers 
the normal cxq range, and it exceeds the normal qQ range. Con- 
sequently, these flight conditions do not literally repre- 
sent coordinated pullup maneuvers, although they introduce 
the same symmetric coupling terms in the F matrix (Eq. (A. 3-4)) 
that occur in the pullup. 



Closed-loop stability at the 15 conditions is sum- 
marized by Table 4.3-1, where it can be seen that all modes 
are stable and at least moderately damped. Control power 
does not change with q Q , but it does change with a Q ; con- 
sequently, the closed-loop stability at a given a Q is rela- 
tively independent of q Q . There is a gradual decrease in 
Dutch roll damping as a Q increases, and roll response 
becomes more sluggish. This happens because rudder and 
aileron are less effective at the higher angles, while the 
elements of R which weight the cost of using these surfaces 
remain unchanged. A heavily damped roll-spiral oscillatory 
mode occurs at a Q of 20 and 25 deg. A coupled roll-spiral 
mode can degrade handling qualities, so adjustment of 

TABLE .4.3-1 

CLOSED-LOOP STABILITY IN THE LONGITUDINAL SWEEP 


Maneuver 

Condition 

Short Period 

Dutch 

Roll 

Roll 

Spiral 

Phugoiri 

“o' 

V 

u • 
n 

C 

V 

C 

T . 

T , 

w n* 

C 

deg 

deg/soc 

rad /sec 

- 

rad /bee 

- 

sec 

ROC 

rad/scc 

— 

5 

0 

2.43 

0.71 

1.90 

0.52 

0.14 

1.18 

5.54* 

8.06* 


12 

2.43 

0.69 

1.95 

0.53 

0.14 

1.23 

.43* 

16.82* 


24 

2.44 

0.67 

2.02 

0.55 

0.14 

1.32 

.65* 

18.31* 

10 

0 

2.60 

0.72 

1.94 

0.57 

0.17 

1.44 

0.17 

0.95 


12 

2.60 

0.71 

1.99 

0.58 

0.17 

1.20 

2.36* 

14.95* 


24 

2.60 

0.69 

2.08 

0.59 

0.17 

1.30 

1.60* 

17.25* 

15 

0 

2.67 

0.72 

2.35 

0.69 

0.36 

1 .07 

0.16 

• 

0.99 

* 


12 

2.67 

0.71 

2.42 

0.67 

0.37 

1.17 

2.31 

18.23 


24 

2.68 

0.70 

2.51 

0.66 

0.34 

1.30 

1.57* 

21.15* 

20 

0 

2.72 

0.72 

2.95 

0.57 

1.0« + 

0.97* 

0.16 

0.98 


12 

2.72 

0.72 

2.97 

0.56 

0.82 

1.06 

2.22* 

18.46* 


24 

2.72 

0.71 

3.00 

0.55 

0.66 

1.34 

1.50* 

20.21* 

25 

0 

2.75 

0.62 

2.22 

j 

0.40 

0.80 + 

0.90* 

6.09* 

9.58* 


12 

2.78 

0.61 

2.23 

0.39 

l 

0.79 + 

0.93* 

2.57* 

67.57* 


24 

2.81 

0.60 

2.25 

0.39 

0.78* 

0.99* 

1.79* 

58.41* 


•2 Real Roots 

tRoll-Splral 
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Q and R could be necessary to eliminate this characteristic. 

The normally oscillatory phugoid mode degenerates I 

into two real modes at most maneuvering conditions consid- 
ered here. The over-damped phugoid mode may result from the 
"cost" associated with Au perturbations, which could be 
relaxed in future DPSAS designs. 


There are 16 non-trivial DPSAS gains generated for 
the pullup maneuver. Scheduling of these gains is dis- 
eased in Section 4.4, and 12 of the gains (four each for 
elevator, aileron, and rudder) are presented here. Table 
4.3-2 lists these gains for a Q of 5, 15, and 25 deg and q Q 
of 0, 12, and 24 deg/sec. The first subscript of k indi- 
cates the control effector (in the order used in Table 4.2-2) 
and the second subscript indicates the feedback variable 
(also ordered in Table 4.2-2). 

Longitudinal Gains (k ^ to k »^) - The gains main- 
tain an orderly progression with both ag and q^; none change 
sign and most follow a single increasing or decreasing trend 
with the two flight variables. Table 4.3-2 shows that the 
pitch-rate gain ^23) is dominant at all maneuver conditions 
and has a maximum variation of less than 25 percent , which 
is representative of the variations of most gains at most 
conditions. 



Directional Gains (k^ to k 4 S ) - Substantial varia- 
tions in rudder gains can be expected with increasing cxq. 

The fuselage blocks the flow over the vertical tail at high 
a 0 , and the rudder side force transforms into stability 
axis roll and yaw moments differently at different angles of 
attack. Unlike the longitudinal gains, there is a dramatic 
change in the directional gains as a Q increases from 15 to 
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25 deg (Table 4.3-2). The yaw gains (k^g and k 4 g) have 
significant changes with both a Q and q Q . The roll-angle 
gain (k^g) is noted to change sign as otg progresses from 
5 to 15 deg, while the increased roll-rate gain (k^) 
attempts to provide stability-axis yaw damping. 

Lateral Gains (k ^ g to k ^g) - Trends in the aileron 
gains also have large variation with due to geometric 
transformation, loss of rudder effectiveness, and aileron 
yaw effects. More 1 gains change sign, and pitch rate has a 
greater effect on gain magnitude. There is an abrupt reduc- 
tion in the use of aileron for roll control (kgy and kgg) 
at an ctg of 25 deg, which is accompanied by increased aileron 
use for yaw control (k^g and kgg). 

It was noted earlier that pitch rate destabilizes 
the Dutch roll, roll, and spiral modes. The coupled nature 
of this phenomenon has an interesting effect on the secondary 
lateral-directional control paths, i.e., the yaw feedback to 
the roll moment controller (and the converse), such as k Qe .» 
kgg, k^y, and k^g at the lower angles of attack. These gains 
have as great or greater variation with q Q as with the 
change from 5- to 15-deg , which is not the case for the 
primary control paths (yaw-to-rudder and roll-to-aileron ) . 

As in the previous section, the performance of the 
DPSAS in the pullup flight condition is assessed by compar- 
ing open- and closed-loop time responses. Figure 4.3-1 
illustrates the aircraft’s open- and closed-loop responses 
to an initial sideslip perturbation when otg is 15 deg and 
q^ is 12 deg/sec. The oscillation grows at a moderate rate 
without stability augmentation but is damped in one cycle 
with the control loops closed. At higher angle of attack 
(25 deg) and the same pitch rate, the open-loop oscillation 
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grows at a faster rate, and the closed-loop oscillation 
takes two cycles to disappear (in keeping with the reduced 
damping ratio of the Dutch roll mode). 

4.3.2 Lateral-Directional Sweep 

Nominal values of sideslip angle and stability-axis 
roll rate are varied in this section, and their effects on 
closed-loop eigenvalues, DPSAS control gains, and aircraft 
response are presented. The development of this section 
follows the previous section, although the results presented 
for asymmetric flight are somewhat different from those of 
the longitudinal sweep. The Q and R matrices are the same 
as before, and linear-optimal regulators are designed at 
18 points. Sideslip angles of 0, 5, and 10 deg are con- 
sidered in combination with stability-axis roll rates of 
0, ±13, ±26, and ±39 deg/sec. (For a given sideslip angle, 
roll rates of opposite sign have different dynamic effects.) 
Angle of attack, velocity, and altitude are fixed at 15 deg, 
94 m/s, and 6100 m, respectively. 

© 

Table 4.3-3 presents the natural frequencies, 
damping ratios, and time constants of the aircraft, with 
the linear-optimal regulator loops closed. The most strik- 
ing result, in comparison with Table 4.3-1, is that the 
lateral-directional closed-loop roots evidence relatively 
little variation with maneuver condition. There are no 
roll-spiral or phugoid degeneracies, and all parameters 
stay within 40 percent of their mean values. Short period, 
Dutch roll, and phugoid natural frequencies decrease with 
increasing 6 Q magnitude and increase with increasing p Wo 
magnitude. Roll time constant and damping of the short 
period and phugoid modes are largely independent of 8 q 
magnitude but decrease with p^ magnitude. Dutch roll 
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TABLE 4.3-3 

CLOSED-LOOP STABILITY IN THE LATERAL-DIRECTIONAL SWEEP 


Man 

Con 

euver 

dition 

Short Period 

Dutch 

Roll 

Roll 

Spiral 

Phugoid 

9 o- 

dee 

p *0' 

deg/sec 

V 

r ad/ sec 

C. 

n* 

rad/sec 

C. 

T . 

sec 

T , 

sec 

V 

rad/aec 

c, 

0 

39 

2.95 

0.56 

2.52 

0.91 

0.51 

1.13 

0.34 

0.77 


36 

2.79 

0.60 

2.35 

0.84 

0.40 

1.09 

0.26 

0.80 


13 

2.64 

0.66 

2.40 

0.75 

0.37 

1.07 

0.19 

0.89 


0 

2.67 

0.72 

2.35 

0.69 

0.36 

1.07 

0.16 

0.99 

5 

-39 

2.95 

0.36 

2.48 

0.93 

0.52 

1.07 

0.35 

0.77 


-26 

2.79 

0.61 

2.28 

0.85 

0.38 

1.06 

0.27 

0.78 


-13 

2.65 

0.66 

2.35 

0.76 

0.38 

1.06 

0.20 

0.87 


0 

2. ,4 

0.72 

2.35 

0.69 

0.35 

1.07 

0.16 

0.99 


13 

2.59 

0.66 

2.45 

0.75 

0.36 

1.09 

0.19 

0.91 


36 

2.76 

0.60 

•> 

0.83 

0.38 

1.14 

0.26 

0.81 


39 

2.92 

0.56 

2.51 

C. 90 

0.46 

1.22 

0.34 

0.78 

10 

-39 

2.87 

0.56 

1.94 

1 0.92 

0.32 

1.06 

0.35 

0.76 


-26 

2.71 

0.62 

2.06 

0.83 

0.31 

1.05 

| 0.28 

0.78 


-13 

2.60 

0.69 

2.16 

0.73 

0.30 

1.06 

0.20 

0.86 


0 

2.63 

0.73 

2.16 

0.67 

0.30 

1.08 

0.16 

0.99 


13 

2.50 

0.71 

2.31 

0.70 

0.31 

1.11 

0.19 

0.92 


26 

2.65 

0.61 

2.25 

0.81 

0.32 

1.18 

0.27 

0.82 


39 

2.82 

0.56 

2.21 

0.89 

0.34 

1.30 

0.35 

0.79 


damping increases with Py; Q magnitude and is little affected 
by 8 q. The spiral mode time constant increases with P\y Q 
magnitude, although its minimum value occurs at more nega- 
tive P^ as increases. 

Examples of the DPSAS gain variations with sideslip 
angle and roll rate are plotted in Fig. 4.3-2 and 4.3-3. 

The most apparent trend is that primary gains, i.e., those 
which would be non-zero in symmetric flight, change very 
little with 8 q and p w ^ while crossfeed gains have sub- 
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stantial variation with maneuver condition. The standard 
deviation of each gain, computed over the 18 lateral - 
directional sweep conditions, is an indication of its 
variation from a constant value. The average standard 
deviation for the primary gains is 16 percent, and for the 
crossfeed gains it is 422 percent. As discussed in Sec- 
tion 4.6, this is a first indication of gain-scheduling 
requirements, suggesting that many primary gains are nearly 
constant and that most secondary gains must be scheduled 
(unless they are negligible). 

Gain variations are seen to depend on whether the 
vehicle is sideslipped "into" or "out of" the roll. (The 
vehicle is sideslipped into the roll when Bq and p^ have 
opposite sign, e.g., when the nose is left and the left 
wing is moving down; it is sideslipped out of the roll when 
the 3igns are equal.) Figures 4.3-2 and 4.3-3 illustrate 
gain variations for positive Bq only; for negative Bq, the 
variations with p WQ are changed. The graphs of primary 
gains for negative Bq are mirror images of those for posi- 
tive Bq (Fig. 4.3-2). The graphs of crossfeed gains for 
negative B Q shift up or down, in opposition to the Bq trend 
shown in Fig. 4.2-3 Primary gains can be monotonic or 
convex functions of p ; crossfeed gains are monotonic in 
P W q and always pass through zero when both Bq and p^ Q 
are zero (Gains for symmetric flight are inidcated by 
in Fig. 4.3-3). 

The crossfeed gains are shown to be non-trivial for 
even moderate values of Bq and P^q, and those shown in Fig. 
4.3-3 can be interpreted as nonlinear control elements . 

Note that each gain could be approximated by a function of 
the form 

Gain = c 1 B () + c 2 P Wq (4.3-1) 
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where and Cg are appropriate constants. Then the con- 
trol signals represented by these four graphs would be 

(4.3-2) 
(4.3-3) 
(4.3-4) 

(4.3-5) 

where the constants are derived by regression analysis 
(Section 4.6). The P^qAw and p^AG terms can be recognized 
as analogous to so-called "pseudo-8" or "pa" crossfeeds, which 
have been incorporated in the SAS of modern high-performance 
aircraft. (An additional "pa"-tyt)e primary gain is indi- 
cated in Table 4.3-2. The roll rate-to-rudder gain could be 
approximated by ca^; therefore, the associated rudder com- 
mand would be cdgAp.) Nonlinearities in the curves of Fig. 
4.3-3 suggest that higher-order fits than Eq. (4.3-1) to 
(4.3-5) are required if design performance is to be obtained 
over a wide range of 8 q and P\v Q - 

Examples of open- and closed-loop response at two 
asymmetric flight conditions are shown in the next two fig- 
ures. Figures 4.3-4 and 4.3-5 show that roll rate intro- 
duces substantial longitudinal response to a directional 
input and that the addition of sideslip angle leads to 
qualitative changes in response shapes. Roll rate alone 
introduces regular oscillations in the aircraft's open-loop 
response (Fig. 4.3-4). The DPSAS damps the oscillation 
within 1$ cycles, although excitation of the phugoid mode 
leads to a slow decay in Aa (The effective time constant 
(-Cu> n ) of the phugoid is 6 sec at this flight condition). 


A5 h * c l 6 0 Ar + c 2 p w 0 Ar 

AS a * °3 6 0 4q + °4 p w 0 aq 
4«r - c 5 B 0 A9 f 

“a ■ c 7 6 0 Aw + c 8 p W 0 iw 
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When the aircraft has developed a large mean side- 
slip angle as well as roll rate, the open-loop perturba- 
tion motions tend to meander, as several modes are involved 
in each motion (Fig. 4.3-5). For example, the initial Aa 
appears to be damping out, but after 5 sec, it begins to 
wander. Pitch and roll angle develop offsets which are 
continuing to increase at the end of the time period shown. 
The DPSAS restricts the maximum initial excursions of A0 , 

A<t>, and AB to less than half their open-loop values and 
eliminates the meandering characteristic. 

Plotting B rather than AB in Fig. 4.3-5 is a 
reminder that the DPSAS provides stab .lity about a refer- 
ence flight condition, in this case, 10-deg sideslip angle 
and -39 .6-deg/sec roll rate. With the assumption that 
these values are commanded by the pilot, it can be seen that 
the DPSAS does not limit aircraft maneuverability — in fact, 
it expands the flight envelope by stabilizing the aircraft 
in conditions which could not be controlled by the unaided 
pilot. Although non-zero Bq is not normally desired in 
maneuvering current high-performance aircraft, future air- 
craft, particularly those with direct side force control, 
could use this capability to tactical advantage. 

This section has presented linear-optimal DPSAS 
designs for the reference aircraft and for a variety of 
maneuvering conditions. The next section of this chapter 
demonstrates how control gains can be adapted to flight con- 
dition. 


Ill 



4.4 


CONTROL-LAW ADAPTATION FOR VARYING FLIGHT CONDITIONS 


This section presents results for a procedure which 
adapts the control gains of a high-performance aircraft to 
varying flight conditions, including gain correlations for 
the reference aircraft model. The gai^s are scheduled by 
finding functional l'elationships between aircraft flight 
variables and the control gains at the corresponding flight 
conditions. 

Previous methods for scheduling control gains have 
been successful and Indicate that gain scheduling is a 
sound approach. The methodology typically is based on 
single input/single output concepts, e.g., maintaining 
constant loop gain. These previous methods, however, pro- 
vide inadequate insight for scheduling a multivariable 
system. 


The method is a logical extension of previous woi k 
to multivariable systems. It involves three steps, and 
it places minimum reliance on past experience and intuition. 
The three steps are: 


• The determination of means and stand- 
ard deviations of the control system 
gains . 

• The determination of correlation 
coefficients between gains and flight 
variables . 

• The determination of functional re- 
lationships (or curve fits) between 
the chosen flight variables and the 
gains . 


This new gain scheduling procedure, discussed in Section A. 4. 5, 
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is simple to use, the results are easy to implement on a 
digital computer , and the procedure can have broad appli- 
cation. / 

The longitudinal and lateral-directional sweep 
control gains discussed in Section 4.3 have been correlated 
with a number of flight variables. In order to identify the 
individual effects of longitudinal and lateral-directional 
mean motions, the correlations for each sweep are done 
separately. In a flight system, the gains should be cor- 
related jointly, and additional factors — such as weight, 
altitude, and velocity — must be considered. 

For each sweep, a list of candidate independent 
variables is established, and various functions of their 
variables are correlated with the 32 gains associated with 
each flight condition. Functions considered included poly- 
nomials of order one and two , 

• 

Gain * b Q + b^ (4.4-1) 

Gain * bQ + b^m + bgm 2 (4.4-2) 

and linear regressions in two variables, 

\ 

Gain * b o + + b 2 m 2 ( 4 . 4 . 3 ) 

Given a flight variable, y, independent variables, m, of 
the form y, y , 1 /y, 1 /y , and y|y| are considered in the 
polynomial regress ions . Equation (4.4-3) is used with 
mj^yj^ and m 2 = y 2 . The objective of the computations is to 
find the functional approximation to each gain which has the 
greatest correlation with the linear-optimal gain at all 
conditions in the particular sweep. In many cases, alter- 
nate functions have similar correlation coefficients, so 
more than one schedule could be considered in implementation. 



113 




The correlation between gains also is of interest, 
as it suggests which gains can be scheduled as functions of 
other gains, and it helps to identify control interconnects , 
o This correlation can be computed, using Eq . (A. 4-37), by 
defining the first gain as k and the second as £. The 
following results indicate that an aileron-rudder inter- 
connect could be considered for the stability augmentation 
system of the reference aircraft. 



The procedure followed in establishing gain-sched- 
uling requirements is to compute means and standard devia- 
tions (as percentages of the means) of tne gains, to corre- 
late gains with flight variable functions, and to corre- 
late gains with other gains. Sixteen crossfeed gains are 
identically zero, leaving sixteen gains for scheduling. 


Table 4.4-1 summarizes the findings for DPSAS 
gains in the longitudinal sweep, presenting the mean and 
standard deviation of each gain. The independent variables 
which provide the best gain schedule are listed, along with 
the correlation between the actual and scheduled gain values. 
Also listed is the gain which exhibits the highest cross- 
correlation, (calculated by applying Eq. (A. 4-37) to all 
pairs of gains) and the value of that cross-correlation. 

For example, the gain A6^,/A0 exhibits a mean of -0.016 and 
a standard deviation of 32% over the chosen set of longi- 
tudinal flight conditions. A gain schedule using normal 
load factor (u Z q) and pitch rate (qg) produces a scheduled 
gain whose correlation factor with the actual gain is 0.89. 
Finally, A6 T /A9 exhibits strongest cross-correlation with 
A 6 t / u, and the correlation factor is 0.98. The flight 
variables considered as possible scheduling variables are 
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TABLE 4.4-1 

GAIN CORRELATIONS FOR THE LONGITUDINAL SWEEP 


Longitudinal Gains 


Gain 

A6 t /A6 

A6 t /Au 

A6 T /Aq 

A<$ t /A\v 

A6 h /A0 

A V 6 u 

A6 h /Aq 

A6^/Aw 

Mean of Gain 

-0.016 

0.062 

-0.006 

0.013 

-0.763 

1.640 

-1.621 

-0.928 

Standard Deviation 
of Gain, % of Mean 

32 

16 

18 

54 

17 

20 

6 

26 

Best Scheduling 
Variables 

n 7, > Q 

n z» q 

n 2 , q 

n z , q 

n z , q 

n 2 , q 

cos a 

cos a 

Scheduled/Actual 
Gain Correlation 

0.89 

0.93 

0.95 

0.95 

0.85 

0.94 

0.88 

0.89 

Gain of Highest 
Cross Correlation 

A6 t /A u 

A6 t /A6 

A6 h /Au 

A6 h /Aq 

A6 h /Aq 

A6 T /Aq 

A<S a /Ar 

A6 & / A<J> 

Gain Cross 
Correlation 

0.98 

0.98 

0.91 

0.89 

0.84 

| 

0.91 

0.92 

0.94 


o 

Lateral Gains 


Gain 

A6 /Av 
a' 

A6 /Ar 
a 

A6 a /Ap 

A6 /A0 

Cl 

A6 r /Av 

A6 / Ar 

r 

A6 r /Ap 

A6 r /A(j> 

Mean of Gain 

-0.446 

0.076 

2 . 054 

1.764 

2.480 

-3.567 

0.404 

-0.206 

Standard Deviation 
of Gain , % of Mean 

251 

1713 

50 

40 

50 

32 

60 

111 

Best Scheduling 
Variables 

cos a 

cos a 

cos a 

cos a 

cos a 

cos a 

"2 

n z ,q 

Scheduled/Actual 
Gain Correlation 

0.93 

0.83 

0.82 

0.90 

0.93 

0.91 

0.94 

0.98 

Gain of Highest 
Cross Correlation 

A6 r /A$ 

A6 a /A<(> 

A6 r /Ar 

A6 r /Ar 

A6 r /Ar 

A<$ r /Av 

A6^/Aw 

i 

A6 t /Au 

Gain Cross 
Correlation 

0.94 

0.98 

0.98 

0.99 

1.00 

1.00 

0.89 

0.90 


♦All independent variables evaluated at nominal flight condition, 
M 0" subscript omitted. 
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angle of attack (a 0 ). pitch rate (Qg). normal load factor 

fn Zo ), sinoiQ, and cosa^ This set of independent variables 
includes those actually varied (a Q and q Q ) and some likely 
functions cf them, and it serves to illustrate the DPSAS gain 
scheduling procedure. 

The only gain which is reasonably constant in 
Table 4.4-1 is A6 h /Aq, the pitch rate-to-elevator gam, as 
its standard deviation is just 6 percent of the mean value. 

As is the case for all gains which are most highly corre- 
lated with cos cxq, the best functional fit is given by 

2 4 

Gain = bQ + b^/cos + bg/cos (4.4-4) 

The roll rate-to-rudder gain, A6 r /Ap, is best approximated 
by a polynomial in load factor, 

Gain = b n + b.n^ + b 9 nj (4.4-5) 

u 1 z o z z 0 

and the remaining gains a-re "best" fit by linear functions 
of n ZQ and q Q . This is not to say that alternate functions, 
e.g., n|g and q Q , would not be better. These correlations 
are "best" for the independent variables and functions con- 
sidered, and the average correlation with actual gain values 
is over 0.9. 

All gains have a correlation greater than 0.89 with 
at least one other gain except AS^/AS, which has a corre- 
lation above 0.8 with several gains. The suitability of a 
SAS aileron-rudder interconnect is evaluated by noting the 
correlation between the aileron and rudder gains for each 
feedback variable (Av, Ar, Ap , and A<J>). The correlations 
(not necessarily the maximums, and not necessarily in the 
table) are 0.93(Av), 0.95(Ar), 0.83(Ap), and 0.70(A<J>), indi- 
cating a good possibility for combining Av and Ar feedbacks 
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with little performance degradation, and a moderate possi- 
bility for combining Ap and A<f> feedbacks as well. 


4.4.2 Lateral-Directional Sweep 

Correlation results for the lateral-directional 
sweep are shown in Table 4.4-2 which indicates that no gain 
means are identically zero, although several appear neglig- 
ible. The flight variables considered for scheduling are 
sideslip angle (Bq) and stability-axis roll rate (p^), which 
are chosen to illustrate the DPSAS gain scheduling method. 

In an actual application, additional independent variables 
could be included in the search. 


Unlike the longitudinal sweep, it appears that 13 

o 

gains could be considered constant , with standard deviations 
of less than 8 percent of the mean value. Five gains are in- 
adequately scheduled by the chosen independent variables and 
functions, as their correlations are below 0.75. Two of these 
are the rudder gains shown in Fig. 4.3-2, which can be seen to 
be more complex than the polynomials and linear combination 
considered here. Higher-order curves would fit these gains, 
although they are candidates for the constant-value approxi- 

O 

mation because their standard deviations are low. 


Seventeen gains are most closely correlated with 

Fty and are fitted best by second-order polynomials in pv^. 

The eleven gains which are most correlated with are fitted 

u 2 

almost as well by second-order polynomials in B Q , B Q , 1 /Bq, 

2 


1/6J, or B 


0 


BqI 


Three of the four linear Bq-P^q fits are 

adequate, with A6 r /Ar requiring an improved fit (along with 

A6 /Aw, A<5 /Av, A6,/Ar, and A6 /Av) . Most pairs have strong 
el.. a a r 

correlation with at Least one other gain. The correlations 
associated with SAS aileron-ruuder interconnect are 0.95(Av), 


( \ 
v 

i 
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TABLE 4.4-2 

GAIN CORRELATIONS FOR THE LATERAL-DIRECTIONAL SWEEP 



All independent variables evaluated at nominal flight condition, "O’* subscript omitted 
































0.50(Ar), 0.13(Ap), and 0.04(A4>), indicating that either 
rudder or aileron would require additional Ar and Ap feed- 
backs in parallel with the interconnected control path. 


/ 


4.4.3 Additional Considerations 

The longitudinal and lateral-directional sweeps were 
conducted to illustrate the separate effects of oiq, qg, Bq, 
and p^ 0 on DPSAS gains. Furthermore, a limited set of in- 
dependent variables and scheduling functions were examined. 

At a minimum, these sweeps should be combined in a single 
correlation procedure to obtain a single multi-variable 
schedule for each gain. Altitude, velocity, and weight 
effects should be added, principally through indicated air- 
speed, Mach number, and the ratio of weight-to-dynamic 
pressure. Permutations of the independent variables, e.g. , 
body-axis rather than stability-axis rates, may provide 
better correlation or may be easier to implement in a par- 
ticular system. 


The present results suggest that primary gains 
schedule largely on longitudinal variables and that cross- 
feed gains schedule primarily on lateral-directional vari- 
ables . This observation derives from the fact that most 
primary DPSAS gains are nearly constant as £q and p^q 
change, while crossfeed gains are zero in symmetric flight. 
Any approximations made in gain scheduling must be validated 
by direct simulation, as this is tantamount to changing 
the gains from their linear-optimal values, thus altering 
closed-loop response. 


An entirely separate issue is the on-board deter- 
mination (either through measurement or estimation) of the 
independent variables to be used for gain scheduling — 
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particularly a Q and 3 Q , which are notably difficult to 
measure. Two potential problems are inaccurate steady- 
state measurement, which leads to inaccurate calculation of 
gains, and superposition of perturbations on the mean values, 
which could cause longitudinal motions to drive lateral- 
directional motions (and vice versa) through oscillatory 
gain changes. The solution to both problems, should they 
occur, is found through s tate estimation , which allows all 
available measurements to be blended in a unified estimate 
of nominal and perturbation motion variables (Ref. 63). As 
an example, measurements of a, q, n , airspeed, and 5. 

Z D 

could be used to estimate a Q , q Q , Aa, and Aq. If the DPSAS 
is incorporated in a full command augmentation system , pilot 
commands could be direct indicators of the desired (or nominal) 
state; therefore, they could be used for gain scheduling 
(Ref. 58). This is a topic for further study. 

4 . 5 CHAPTER SUMMARY 

This chapter has presented design principles for 
stability augmentation systems (DPSAS) which prevent depar- 
ture from controlled flight. Linear-optimal control theory 
has been used to develop control structures for departure 
prevention, and the effects of maneuvering condition on 
optimal feedback and crossfeed gains have been explored. 
Examples of aircraft response to longitudinal, lateral, and 
directional initial conditions illustrate the well-controlled 
behavior which the DPSAS provides, and closed-loop eigenvalues 
show that variation in aircraft dynamic characteristics is 
minimized for a wide range of maneuvering conditions. 

In many respects, stabilizing the reference aircraft 
in a pullup maneuver is a more challenging task than accounting 
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for the coupling which results from sideslip and roll rate; 
however, lateral-directional maneuvering results in signifi- 
cant linear-optimal gains which improve aircraft response. 

In combination with gain-scheduling functions which depend 
on mean values of angles and angular rates, the DPSAS con- 
trol algorithms are seen to produce nonlinear crossfeeds 
which are analogous to control structures being employed 
in modern high-performance aircraft . 

Linear-optimal control theory solves many aircraft 
control problems which have been difficult to overcome with 
past design techniques. It is easy to use, it guarantees 

e 

system stability, and it accomnodates aircraft with limited 
control authority. 


! 


1 - . T 


5. 


CONCLUSIONS AND RECOMMENDATIONS 


This report has illustrated how linear systems • 

analysis can be used to characterize the stability of air- 
craft during maneuvering flight. 1+ also presents a design 
procedure for stability augmentation systems which prevent 
departure fiom controlled flight. The key to linearizing 
the dynamics of the aircraft is that an accelerated flight 
condition can oe used as a reference path . A linear model 
can provide a good description of the aircraft's perturbation 
response (to initial conditions, control inputs, and dis- 
turbances) even when the aircraft has large aerodynamic angles 
and ftngular rates. Control systems designed for fully 
coupled linear models and adapted to changing flight con- 
ditions can provide protection against inadvertent depar- 
ture from controlled flight. 


5.1 CONCLUSIONS 


A detailed examination oi the dynamics of the ref- ( 

erence aircraft has led to generalizations concerning air- 
craft stability and control. These include the following: 


• The aircraft's stability (as shown by its 
eigenvalues) is most affected by changes 
in the nominal longitudinal variables 
(V Q , a Q , and Qq) , while the mode shapes 

(as described by the aircraft eigenvectors) 
are most affected by non-zero nominal 
values of the lateral variables (8 q and 
p wq ) ■ Asymmetric flight leads to iongi- 

tudinal-var iable response in typically 
lateral-directional modes, and vice-versa. 
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Nonminimum-phase zeros in the aircraft 
transfer functions occur frequently in 
asymmetric flight, and the transfer 
function numerators are changed sub- 
stantially by non-zero q 0 . 

Extreme maneuvers are often characterized 
by rapid changes in both mode stapes and 
speeds due to large values of a Q , p^, 

and qQ. Highly coupled, unstable natural 
modes can result. 

Elementary loop closures wjiich are sta- 
bilizing in symmetric flight can lead to 
unstable system dynamics in asymmetric 
flight . 

The departure parameter, ^ng^yn’ has 

limited value in predicting aircraft 
departure. It provides no information 
regarding Dutch roll damping; hence, it 
does not predict departure due to nega- 
tive damping (as is the case for the 
subject aircraft). 

Unforced departures occur when one of 
the fast modes (short period, Dutch roll, 
or roll mode) is unstable. These depart- 
ures primarily take the form of a Dutch 
roll instability, with fast rolling-yawing 
motions or oscillatory divergence. The 
roll mode can become unstable at extreme 
angles, where it exhibits a fast roll- 
yaw divergence. 

Forced departures occur as a result of 
pilot action. This can happen when de- 
graded response to control inputs causes 
the pilot to fly the aircraft into a 
flight regime where unforced departures 
are likely or when pilot actions de- 
stabilize the aircraft directly. 

Guidelines for the design of Departure- 
Prevention Stability Augmentation Sys- 
tems (DPSAS) have been presented. An 
adaptive-control design procedure, using 
the linear-optimal regulator for 
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fixed-point design followed by gain 
scheduling, is shown to provide a non- 
linear control structure containing 
crossfeeds as well as feedback gains. 

The resulting DPSAS has similarities 
to the flight control systems of cur- 
rent high-performance aircraft. How- 
ever, the new design is based on "quad- 
ratic synthesis" techniques, which 
provide a unified set of control gains 
for all axes from a single set of 
vector-matrix design equations. 

• The linear-optimal DPSAS prevents departure 
not by limiting the maneuvering ability of 
the aircraft but by stabilizing the air- 
craft in all foreseeable maneuver conditions. 

• The linear-optimal control law can be readily 
extended to a full Departure-Prevention 
Command Augmentation System (DPCAS) which 
accounts for control actuator rate limits 
and allows essentially unlimited pilot con- 
trol authority (within che physical limi- 
tations of the aircraft). 

• The maneuverability envelope of the subject 
aircraft could be materially expanded through 
the incorporation of DPSAS/DPCAS concepts, 

as identified in this report. 


5 . 2 RECOMMENDATIONS 


The following recommendations are made as a result 
of this study: 


• Departure prevention studies for high- 
performance aircraft should be extended 
to transonic and supersonic flight re- 
gimes . 

• The high angle-of-attack/high angular 
rate problems of additional aircraft 
types, including transports, helicop- 
ters, and general aviation aircraft, 
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are amenable to coupled linear analy- 
sis and bear investigation. 

• Design requirements for a DPCAS should be 
investigated. Digital implementation and 
incorporation of active control concepts 
for improved maneuverability should be 
considered . 

• It is recommended that improvements to the 
subject aircraft's maneuverability envelope 
due to DPSAS/DPCAS implementation be explored 
in ground-based piloted simulation and flight 
test . 


Dynamic coupling is a significant factor in the 
maneuvers of high-performance aircraft, and a full under- 
standing of its effects is an important facet of preventing 
departure from controlled flight. This report has shown 
how linear models of the aircraft's motions can be used to 
investigate the stability and control of maneuvering flight, 
and it has demonstrated the flexibility and ease with which 
linear-optimal control theory can be used to design departure- 
preventing control systems. 



APPENDIX A 


ANALYTICAL APPROACH TO AIRCRAFT DYNAMICS 


A. 1 OVERVIEW 

This appendix develops the analytical approach taken 
in this report. The goal of this work is to analyze the air- 
craft stability and control problems that arise at high 
angles of attack and during rapid maneuvers. The nonlinear 
equations of motion of a vehicle in atmospheric flight are 
developed in Section A. 2 . A large body of theory and experi- 
ence relating to linear system analysis and control synthesis 
exists (Ref . 64 to 67), yet a full and complete linearization of 
the aircraft problem is not readily available; therefore, a 
rigorous application of linear analysis should provide new 
insights regarding departure. Section A. 3 presents the full 
linearized equations of motion in a general form, and it 
includes a discussion of methods for choosing the point of 
linearization. Section A. 4 presents an overview of linear 
system analysis and control methods. 


A. 2 NONLINEAR EQUATIONS OF MOTION 

The complete nonlinear rigid-body equations of 
motion are derived in this section. They are developed using 
"flat-earth” assumptions, i.e. the effects of earth curva- 
ture and rotation are assumed negligible. This means that 
earth-fixed and inertial reference frames are equivalent. 
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There are four coordinate systems of interest in 
the study of rigid-body motions of aerodynamic vehicles. 
They are described as follows: 


• Inertial-Axis System - This frame is 
fixed in inertial space, and is the 
frame from which the inertial velocity 
and angular rare of the vehicle are 
measured. 

• Body-Axis System - The forces and 
moments on the vehicle, and there- 
fore the dynamic equations, are best 
expressed in a body-fixed reference 
frame. The stability-axis system is a 
special body-fixed reference frame. 

• Velocity-Axis System - An especially 
useful reference frame for navigation 
and guidance, the velocity-axis system 
relates the vehicle velocity vector to 
inertial axes. 

• Wind-Axis System - Since the aerodynamic 
forces and moments depend largely on the 
body-velocity vector orientation, the 
wind-axis system, which provides straight- 
forward body-wind relations, is useful. 


For the moderate velocities of interest in this 
report (typically subsonic, i.e., below about 340 m/s), 
the equivalence of earth-fixed and inertial reference 
frames is a good assumption. The origin of the in- 
ertial reference frame used here is located on the sur- 
face of the earth, with the x- , y-, and z-axes in a north- 
east-down orientation. Since the simplest statement of 
Netwon's Second Law is given in an inertial reference 
frame, this frame plavs an important nart in the derivation 
of the dynamic equations. 


The various body-fixed axis systems have a common 
origin, located at the body center of mass, and are fixed 
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in orientation with respect to the vehicle. Generally, the 
body x-axis extends forward out the vehicle's nose, the 
y-axis extends out the right wing, and the z-axis extends out 
the bottom of the vehicle. The x-z plane is usually a plane 
of geometric symmetry, if the vehicle has one. There are a 
number of possible body-fixed reference frames, and the one 
fixed by the builder is simply referred to in this report 
as the body-axis system . For any nominal flight condition, 
body-fixed axes can be chosen so that the x-axis is aligned 
with the velocity vector, and the z-axis is in the body- 
axis x-z plane. This set of body-fixed axes is referred to 
as the stability-axis system . 

Since body axes are the only axes in which the 
vehicle rotational inertia matrix is constant, the rotational 
dynamics equations are usually (though not exclusively) 
expressed in this frame. The body frame also is the one in 
which the pilot, and all sensors and control surfaces are 
located; for this reason, the body frame is considered the 
basic frame of reference in this report. Figure A. 2-la 
illustrates the body frame orientation with respect to the 
inertial-axis system. 

The velocity - and wind - axis systems have a common 
origin (the vehicle center of mass) and a common x-axis, 
which is oriented along the vehicle's inertial velocity vec- 
tor. The velocity reference frame y-axis is parallel to the 
inertial x-y plane. This results in simple relations be- 
tween inertial and velocity axes, so that these axes are 
useful for navigation and point-mass trajectory calcula- 
tions. Figure A. 2-lb illustrates the relationship between 
inertial and velocity axes. 

The wind reference frame's z-axis is located in 
the x-z plane of the body frame; this reference frame is 
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b) Inertial-Velocity Orientation 



c) Wind-Body Orientation 
Figure A. 2-1 Reference Frame Relations 
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very useful in dynamic calculations because the orientation 
angles between the wind frame and body frame have large 
influences on the aerodynamic forces and moments. Figure 
A.2-lc illustrates the orientation between body and wind 
axes. Figure A. 2-2 summarizes the transformations between 
reference frames. Any necessary transformation can be 
identfied from this figure, noting that the Euler angles 
are given in the order of yaw, pitch, and roll, as speci- 
fied by the arrows. For example, a transformation from 
inertial to body axes is composed of a right-handed yaw 
through an angle, <J> , then a right-handed pitch through 
an angle, 0, and then a right-handed roll through an angle, 
<J>. These three single-angle transformations can be com- 
bined to form an inertial-body transformation as follows: 

, 8 , <J0 =H®(<*>) H^( 0 ) Hj (ijj ) 



_ 1 

0 

0 " 

~COS0 

0 

-sin8" 

~cos\)> 

sinijj 

0 " 


0 

cos<J> 

sin<J» 

0 

1 

0 

-simp 

COSlj> 

0 


_0 

-sin4> 

cos<J> 

sin8 

0 

COS0 

0 

0 

1 _ 


(A. 2-1) 

For orthogonal matrices such as these, the matrix inverse, 

-1 T 

( ) , is equal to the transpose, ( ) . ^ 

V 

In the remainder of Section A. 2, the vehicle's 
equation of motion is derived as a single state-vector 
equation of the form, 

x = f(x,u) (A. 2-2) 

where x is the state vector, u is the control vector, f is 
the vector system dynamics equation, and disturbances are 
neglected. The state vector is a 12-element vector, and 
the nonlinear state equations are readily derived as four 
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Figure A. 2-2 Reference Axis Transformations 
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sets of three equations representing 

• Translational Kinematics 

• Rotational Kinematics 

• Translational Dynamics 

• Rotational Dynamics 











The kinematic equations relate the vehicle's translational 
and rotational velocities to its position in inertial space, 
and they involve body-axis/inertial-axis relationships. The 
dynamic equations describe the changes of the vehicle veloc- 
ities caused by the applied forces and moments; they are best 
derived in a body-fixed frame of reference. 

A. 2.1 Kinematics 

Kinematics is the study of the motion of a body 
without regard to the forces which cause that motion. In 
this section, the relations between the vehicle’s position 
and velocity are examined. The translational and angular 
position of the vehicle are given relative to inertial space 
by the inerti*al position vector, xj , and by the inertial- 
body Euler angle vector, 

(A. 2-3) 


(A. 2-4) 

It is important to note that the Euler angle "vector" is not 
a x rue vector in physical space; it is an ordered triple cf 
right-handed rotations which occur about different axes of 
different reference frames. 

The translational and angular rate vectors are 
often expressed in body axes, as in the following: 
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The body-axis translational rate vector, v Q , is an expression, 
in body axes, of the derivative of the inertial position vec- 
tor. This relationship supplies the first part of the non- 
linear state equations of notion: 


-I * H B — B 


(A. 2-7) 


where Hg is the inverse of the inertial-body transformation 
derived in Eq. (A. 2-1). 

The body angular rate vector also can be related to 
the derivative of the Euler angle vector by noting that the 
Euler angle derivatives occur in three different reference 
frames. The resulting transformation is constructed in 
Eq. (A. 2-8), where the individual transformations are the 
same as those of Fig. A. 2-2 and Eq. (A. 2-1): 


4>] TO] ("0 

0 + H® 6 + H®H* 0 

oo tii 


* S^B 

The non-orthogonal transformation, Lg, is 


(A. 2-8) 


1 0 < -sin0 

0 cos$ sin4>cos8 
0 -sin$ cos4»cos0 


(A. 2-9) 





The ordering of the transformations in Eq. (A. 2-8) arises 
from the ordering of the Euler angles. As can be seen from 
Fig. A. 2-la, the angular rate, 0, occurs about the Xg axis, 
while the rate, 0, occurs about the y ^ axis, and iji occurs 
about the axis. The inverse of Eq. (A. 2-8) supplies the 
rotational kinematic part of the vehicle nonlinear state 
equations, and is given by: 



(A. 2-10) 


The relations between translational position and velocity, 
Eq. (A. 2-7), and between angular orientation and velocity, 
Eq. (A. 2-10), comprise the kinematic portions of the non- 
linear state equations. 


A. 2. 2 Dynamics 

The dynamics of the vehicle involve the interaction 
between the vehicle motion and the forces that produce that 
motion. This involves the application of Newton's Second 
Law, which equates the applied force to the time derivative 
of inertial translational momentum of a body. For rota- 
tional motion, this equivalence becomes one between torque 
and the derivative of angular momentum, measured in an inertial 
reference frame. 


An expression for the inertial translational accel- 
eration, expressed in body-axis variables, is needed. 

This can be derived from Eq. (A. 2-7) by taking the derivative 
of both sides. It is important to note that the transforma- 
tion matrix is time-varying. This result is 

*1 - H B^B * (A. 2-11) 
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(A. 2-12) 


and Ug is the cross-product equivalent matrix for u>g given by 


~I 

“b 


' 0 -r q 
r 0 -p 
-q p 0 


(A. 2-13) 


This leads to the body-axis equation.* 


-B * H I-I “ “B^B 


(A. 2-14) 


The applied specific forces consist of gravita- 
tional forces and aerodynamic forces. The gravity force is 
especially simple in inertial axes, as it is confined to the 
vertical axis: 



(A. 2-15) 


The specific contact force can be broken into two components, 
one of which is due to aerodynamic forces and one of which is 
due to thrust : 


i B / m 


X/m" 
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Z/m 


(A. 2-16) 
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(Capital letters are conventionally used in aerodynamics to 
denote the force components.) 

The translational dynamic equation is formed by 
equating the sum of the aerodynamic and gravitational spe- 
cific forces to the inertial translational acceleration of 
the vehicle. When all vectors are expressed in body axes 
and the derivatives of the body- axis velocities are iso- 
lated on the left-hand side, the vector equation is 

IB * + + ~ 5 b2b (A. 2-18) 


To construct the rotational dynamic equation, an 
expression for the time derivative of angular momentum mea- 
sured in inertial axes is necessary. The angular momentum, 
hg, is most easily expressed in body axes; neglecting rotat- 
ing machinery, it is the product of the moment-of-inertia 
matrix (constant in body axes) and the angular rate vector, 

h B - I B h b (A. 2-19) 


where the inertia matrix contains all products and moments 
of inertia: 



-I 

-i 


xy 

xz 


I 

-I 

xy 

y 

yz 


-i 

i 

x z 

yz 

z 


(A. 2-20) 


136 



I 


The time derivative of the a.ngular momentum, expressed in 
inertial axes, is easily derived by noting that the trans- 
formation Hg is time-varying: 

*1 * H B T b2b + ‘‘b’bHb 

“ h b^b-b + h b w b^ B— B (A. 2-2: 


The contact moments consist of aerodynamic and 
thrust components. These are defined as 


Gp 

°bL 


(Capital letters are conventionally used for the moment 
components.) The rotational dynamic equation can be formed 
by equating the applied torques to the derivative of the 
angular momentum. All vectors are expressed in body axes, 
and the derivative of body-axis angular rate is isolated 
on the left-hand side of the equation, resulting in 


J B + 2b ) “ *B “B^-B 


A. 2. 3 Summary of Equations 


Section A. 2 has presented the various reference 
frames of interest and has dei ived the equations of motion 
of an atmospheric vehicle. The 12-element state vector con- 
sists of three positions, three angular orientations, three 
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translational rates and three angular rates. The state 
equations were found by examining the translational and 
rotational kinematics and dynamics, and are repeated here: 


*j * H l v B (A. 2-25) 

v B « L^ 1 Wg (A. 2-26) 

“ <lB + !B)/m + H ?£i (A. 2-27) 
2B " -'b^b'bSb (A. 2-28) 


These equations fall into the general state equa- 
tion form , 


X = f (x ,u) 


(A. 2-29) 


by defining the state vector as 

T T T T T t t 
* * I Si — B 2b — B 


(A. 2-30) 


and noting that the aerodynamic forces and moments are func- 
tions of the states, controls, disturbances and, to some 
extent, the state time history. These nonlinear state equa- 
tions are useful because they are general enough to allow a 
thorough analysis of the departure prevention problem. They 
are expressed in state-space form, which is notationally 
efficient and which makes the subsequent linearization an 
easily followed process. 


A. 3 T.INEAR EQUATIONS OF MOTION 

While the nonlinear equations derived in the pre- 
vious section can be solved on a digital computer, they are 
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not easily analyzed by general techniques, and general 
closed-form solutions cannot be obtained. Many of the im- 
portant dynamic attributes of the aircraft can be pre- 
served and the analysis facilitated by developing cor- 
responding linearized equations of motion, as is done in 
Refs. 64 to 67. 


A . 3 . 1 Derivation from Nonlinear Equations 


The linearization procedure begins with the con 
struction of a Taylor series expansion representing the 
nonlinear equations about some nominal trajectory: 


x = x Q + ax 






3f 

Ax + — 
- 3u 


Au + Higher Order Terms 


s-*o 

“ = B 0 


(A. 3-1) 


where the subscript ”0” indicates tne nominal value and the 
prefix "A" denotes a small perturbation. All except first- 
order terms are then neglected by arguing that the higher- 
order terms are small compared to linear terms. The results 
of this procedure are separated into a nonlinear equation 
describing the nominal trajectory (Eq. (A. 3-2)) and a linear 
equation defining the dynamics of the pertu r bations about 
the nominal trajectory (Eq. (A. 3-3)): 


-0 “ (A. 3-2) 

Ax = F Ax + G Au ( a o on 
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and 
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***0 

U-Uq 


G = t- 


3f 

3u 




H“Mo 


(A. 3-4) 


(A. 3-5) 


The linearization is straightforward because the 
nonlinear state equations (Eq. (A. 2-25) to A. 2-28)) are spe- 
cified in a general state-space format. Equations for the 
perturbations of the axis transformations are easily derived 
by taking the partial derivatives of each transformation with 
respect to the Euler angles of that transformation and mul- 
tiplying by the Euler angle perturbations. For the inertial- 
body transformation and its inverse, these transformation 
perturbations can be stated as follows: 


. „B . .. „B 

AHj « -Ap B 


(A. 3-6) 


where the cross-product operator (”) is employed, 


1 Mb - L Bq Av b 

and 

“4 * H B 0 4D b 


(A. 3-7) 
(A. 3-8) 


These relationships for the transformation perturba- 
tions are used to linearize the translational kinematic 
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equation, Eq. (A. 2-25), to give 

AXj = -(T^) H* Lg^ Av b + hJ Av b (A. 3-9) 

This equation clearly shows that the perturbation inertial 
velocity depends both on the perturbation body-axis velocity 
and the inertial-body Euler angle perturbations. 


The rotational kinematic equation, Eq. (A. 2-26), 
results in the following linear perturbation equation: 


B * ' L bJ L B q + 4 “ B 


(A. 3-10) 


where 


9 < L B V 




B 


(A. 3-11) 
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so that 



0 o tan0 Q ^ Q sece 0 0 

-^ 0 cos8q 0 0 

§Qsec0 o ^ 0 tan8 0 0 


(A. 3-12) 


The definition of LBq takes the form it does because the 
linear rotational kinematic equation was derived from the 
original form of the rotational kinematics equation, given 
in Eq. (A. 2-8). 
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Linearization of the dynamic equations requires 
consideration of the aerodynamic force and moment relation- 
ships. These are functions of the states, controls, and 
the past history of these variables. This dependence on 
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past values is caused by aerodynamic flow field effects 
and their propagation delays; unsteady aerodynamic effects 
can be represented as functions of the state rates. The 
formal linearization of the aerodynamic forces and moments 
is a lengthy but straightforward process which amounts to 
taking the partial derivatives of every aerodynamic force 
and moment vector with respect to the states, state rates, 
and controls. These partial derivatives are called 
stability derivatives . 

The difficulty revolves around the actual values to 
be used for each of these coefficients. This data is pro- 
duced primarily by wind tunnel testing, as described, for 
example, in Refs. 68 and 69 . There is a large amount of effort 
and expense involved in generating this data, so only the most 
important functional relationships can be examined. This 
often results in different data sets for each aircraft. 

For this reason, only general terms for the perturbation 
forces and moments are included in the following discussion, 
appendix B contains a discussion and example of the con- 
struction of force and moment stability derivative matrices 
from real data. Many stability derivative matrices are 
either known to be zero or are so small as to be neglected 
in all cases of interest. Assuming that altitude and orien- 
tation variations have negligible effect on contact foice 
and moment variations, and assuming insignificant state 
derivative and angular rate effects on thrust forces and 
moments, the perturbation aerodynamic forces and moments are 
as follows: 

A *B = 

(A. 3-13) 
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The linear translational dynamic equation is derived 
from Eq. (A. 2-27) and incorporates the perturbation aero- 
dynamic force expressions presented above. Note that the 
state-rate stability derivative matrices enter the equation 
in a different manner than the other terms; they must be 
moved to the left-hand side of the dynamic equations. The 
linear translational dynamic equation becomes 



(A. 3-17) 


143 

i 






The linear rotational dynamic equation is derived 
similarly, and the state-rate stability derivatives appear 
in the same way. The linear rotational dynamic equation is 



(A. 3-18) 


The four state equations (translational and rota- 
tional kinematic and dynamic equations) can be written in 
standard linear equation form (Eq. (A. 3-3)) by using the 
following state vector: 

Ax^ * £ax^ AVg AVg AWg J (A. 3-19) 

c 

The state equations then fall into the form: 

JgAXg « Fg AXg + Gg Au (A. 3-20) 

where the state-rate transformation matrix is 


■'B 


10 0 0 

0 10 0 

0 0 J 33 J 34 

0 0 J 43 J 44 


(A. 3-21) 
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and the three-by-three sub-matrices are 


J 33 * 1 


i \!h\ 

m L 3 -bJ 


'34 


J 43 * 
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.1 PMbI 
B L 3 — bJ 


J 44 “ 1 " 1 
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B .-I 
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The primed state dynamics matrix is 
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F 12 

F 13 

0 

0 

*1 

to 

to 

0 

F 24 

0 

F 32 

F 33 

F 34 

0 

0 

F 43 

F 44 


where the three-by-three sub-matrices are 
F 12 

,1 


-L H l L 
*0 B 0 B 0 


F 13 “ «B, 


22 


‘ **0 ^0 


F 24 * L 1, 


£ 


(A. 3-22) 


(A. 3-23) 


(A. 3-24) 


(A. 3-25) 


(A. 3-26) 


(A. 3-27) 
(A. 3-28) 

(A. 3-29) 
(A. 3-30) 

145 




T 


T 


r 


'32 ’ H ! 0 8i H b 0 L B 0 


33 


1 

pvi 

+ 1 

" 9 1b" 

~I 

m 

i 

CQ 

>1 

CO 

m 

_ 3 — B_ 

B o 


34 


1 



m 

- 9 4- 

+ 

B o 


F = I 
r 43 B 


3 — B 


3G B 

3 ^B 


F = I 1 
f 44 B 


L 3 ^B-< 


+ t b^b. 


\ IB ‘ 


and 


G B • 


0 

0 

G3 

L g 4. 


f 

(A. 3-31) r 


(A. 3-32) 


(A. 3-33) 
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(A. 3-35) 


(A. 3-36) 


The sub-matrices (of three rows and as many columns as con- 
trols) are 
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(A. 3-37) 



(A. 3-38) 


The complete state equation is produced by pre- 
multiplying Eq. (A. 3-20) by the inverse of the state-rate 


146 

f 





transformation matrix, giving the following result: 


4 b = J b 1f b 4 *b + J i lG B 

■ F b 4 *B + °B 


. (A. 3-39) 


This resulting linear system is analyzed through- 
out this report. It is important to note that this sys- 
tem specifies the state and control perturbations about their 
nominal values. Methods of properly choosing these nominal 
values are examined next . 

A. 3. 2 Generalized Trim Conditions 

An aircraft is in the trimmed condition when its 
controls are set to produce equilibrium in the equations of 
motion. Steady trim occurs when the aircraft is under no 
inertial acceleration, i.e., when translational velocities 
are constant and rotational rates are zero. The trim con- 
cept can be extended to dynamic flight conditions by defin- 
ing generalized trim as the condition in which control set- 
tings produce constant velocities and angular rates. In 
this case, the vehicle is not necessarily in steady equilib- 
rium, due to changing roll and pitch angles. 

The importance of these trim classifications lies 
in the use of trimmed flight conditions as nominal trajec- 
tories for linearization, since trim implies that the nomi- 
nal velocities, angular rates, and controls are constant, 
or, at most, slowly varying. Thus, the use of the trimmed 
condition as a nominal flight condition causes the linear 
equations to represent almost all of the system dynamics 
at that flight condition. This can be seen by writing the 
general equations for the total state rate: 



(A. 3-40) 
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% + a - b 
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(A. 3-41) 


If a set of nominal states and controls can be found so that 
the generalized trim condition • 


• -I 

X B = =0 

B 0 B 0 


(A. 3-42) 


is satisfied, then the perturbation state rates are equal 
to the total state rates. This can be seen by inserting 
Eq. (A. 3-42) into Eqs . (A. 3-40) and (A. 3-41). 


Another desirable characteristic of using the 
trimmed flight condition as a nominal for linearization is 
that the total state and control trajectories over a signi- 
ficant interval of time are the sums of the constant nominal 
values and the linear perturbation time histories. This 
relation is given as: 


X B (t > 

Ug(t) 

u(t) 
where vg^ , 





+ Av R (t) 
0 B 
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0 ” B 

+ Au( t ) 
and Uq are 


constant . 


(A. 3-43) 

(A. 3-44) 
(A. 3-45) 
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A solution method for the generalized trim problem 
can be derived by examining genera*! forms oi the nonlinear 
translational and rotational dynamic equations, which are 
derived in Sect, n A. 2. 2 as 


— B * — B ' Hb- 

4 * *2<V 4' H) 


(A. 3-46) 
» (A. 3-47) 


The aerodynamic forces and moments are assumed to be func- 
tions of Vg, u)g, and u. One element of the Euler angle vec- 
tor, , does not appear in the equations, and altitude 
is a parameter. 


Two different generalized trim problems then become 
possible, the first of which is the following: 


• Find the values of velocity and angular 
rate (vg and w^) that produce general- 
ized trim (vg = i£g® 0 ) for given controls 
and Euler angles (u, v B ). 


This problem consists of six equations (.the dynamic equa- 
tions) in six unknowns (the velocity and angular rates), 
and therefore it car. be expected to have a solution. The 
controls and Euler angles are the set points that determine 
the generalized trim solution, and it should be noted that 
altitude also has an effect on the solution. 


The second generalized trim problem can provide 
the trim controls that produce specific state values: 




Find the values of the controls and 
Euler angles (u and vg) that satisfy 
the generalized trim conditions 
(vg=w^ = 0) for riven val les of the 


velocity and 
(vg and <4). 


angular rate vectors 
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The existence of a solution to the second problem 
depends on the degrees of freedom and power of the controls. 
Many aerodynamic vehicles have a four-element control vec- 
tor; along with the Euler angles, <t> and 8, this results 
in a problem of six equations (the dynamic equations) in six 
unknowns (throttle, elevator, aileron, rudder, pitch angle, 
and rol 1 angle ) . 


One approach for solving either of these generalized 
trim problems is to use functional minimization . This 
approach, also called parametric optimization, requires all 
elements of u and Vg (or of v^ and ) to be specified as given 
or desired, but it ' s not tied to any particular flight path. 

In this case, Eqs . (A. 3-46) and (A. 3-47) are solved directly 
using an iterative process, e.g., a steepest-descent , 
accelerated gradient algorithm, or direct numerical search 
(Eef. 70). A scalar cost function, J, measures trim error; 
a quadratic form is appropriate for computing a norm of 
the vector error: 


•T 




In both generalized trim problems, the trim value is 
determined when J reaches a mini""’",. 


A . 3 . 3 Body-Axis Equations 

It is advantageous to express the vehicle state 
equations in oody-fixed axes. These are the axes in which 
the pilot, the sensors, and the control surfaces are located. 
Body axes are the only axes in which the moment -of - inert i a 
matrix is constant. Aerodynamic data collected from sting- 
mounted wind tunnel models or from flight test usually is 
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expressed in body axes. Consequently, body axes are con- 
sidered to be the basic axes in this report, and all of the 
equations given so far have been in body axes. 

A . 3 . 4 Stability-Axis Equations 

Stability axes also are body-fixed axes, so that 
they retain the convenient characteristics mentioned above. 

The axes are fixed in the body so that the x-axis is aligned 
with the nominal velocity vector. Although the axes are 
fixed in the body, they have different orientations at 
different nominal flight conditions. Further, the pertur- 
bation velocity vector is not expressed as three orthogonal 
velocity perturbations but as a velocity magnitu< pertur- 
bation and two body-velocity orientation angles: pertur- 

bation angles of attack (Aa) and sideslip (A$). This ”ec- 

tor is referred to here as Av^. it is important to note that 
stability axes are also the same as wind axes for a specific 
nominal flight condition: uo'^o^O. This is the justification 

for using a "W" subscript throughout this report for stability 
axis variables. 

Stability axes simplify certain aspects of the 
linear equations. In level nominal flight, stability x- 
and y-axes are horizontal and the z-axis is vertical. The 
lift and drag forces act along stability axes. Finally, it 
often is true that stability axes are close to the normal 
mode axes. This means that the basic modes of motion 
appear as motion about or along a single axis of the sta- 
bility reference frame. The frequency and damping of the 
basic modes then become simple functions of the aerodynamic 
stability derivatives expressed along stability axes; there- 
fore aerodynamic stability derivatives expressed in the 
stability-axis frame can be used as approximate indications 
of stability or instability. For example, assuming that 
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body axes coincide wi^h principal axes, the stability-axis 
rotational dynamic equations appear in part as: 


% = § pV 0 Sb T~ cosot O C l6 + r S sina o c ng AC + 


(A. 3-49) 


Af W - i pV o Sb T~ z cosa 0 Gn 8 " if sina 0 C l6 AS+ - • • 


(A. 3-50) 


where Appendix B details the definition of the individual 
symbols. The departure parameter C D g ^^(Section 2.2) is 
recognized as an element in the yaw equation. This is sig- 
nificant because these equations are derived by a simplification 
of the complete equations, and it is possible that other use- 
ful departure parameters can be derived from the same approach. 

Since stability axes are body-fixed axes with a par- 
ticular nominal orientation, it is convenient to derive the 
linear equ ti~>ns in body axes and simply rotate them to 
obtain a stability-axis set. This can be done by applying 
the transformation matrix, Kg , to the body-axis system 
matrices, Fg and Gg, as follows: 


% ■ K ™ f b k 0» 


(A. 2-51 ) 
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and, from Fig. 2.2-2, the wind-body transformation is 



COSCXq COSBq 

sinB 0 

sina Q cosBq 


-cosa Q sinB Q -sina Q 
cosBq 0 

-sina Q sinB 0 cosa Q 


while 



sint *o 

0 

-COSOq 


0 

1 

0 


(A. 3-54) 


(A. 3-55) 


and Jy ? 0 is a diagonal matrix whose elements are Cl . Vq, VqCosBq]* 
These relations (which assume aQ=S 0 =0) are used to transform 
body-axis equations to stability-axis equations for further 
analysis . 


A. 3. 5 State Ordering and Dimension 

The order of the states given by Eq. (A. 2-30) is 
the one that proceeds from considering translation before 
rotation and kinematics before dynamics. This order does 
not, however, group related states together. For example, 
many aircraft demonstrate a natural mode that is primarily 
composed of Au and A0 oscillations (the phugoid mode), and 
it is logical to regroup the states so that au and A6 fall 
next to each other. 

For some reference flight conditions, a further 
major division between longitudinal and lateral-directional 
variables can be made. The former variables describe motion 
withii. the vehicle plane of symmetry, while the lateral- 
directional variables describe motion out of the plane of 
symmetry. This division is useful because it allows a 
quick appraisal oi the extent and nature of cross couplings 
that arise in maneuvering flight. 



With this condition in mind, the ordering of states 
given in Table A. 3-1 is suggested for aircraft and similar 
vehicles. The six longitudinal states are first, followed 
by the six lateral-directional states. 

TABLE A. 3-1 
STATE ORDERING 


Body Axes 

Stability Axes 

Ax j 

Ax J 

Azj 

A 2 j 

A9 

A9 

Au 

AV 

Aq 

Aqw 

Aw 

Aa 

Av 

A6 

Ar 

Ar w 

Ap 

Ap w 

A4> 

A$ 

A4> 

Aty 

Ay x 

Ay : 


In any reduced-order approximation, this state order- 
ing also is useful. The four "outemost" states, for example, 
do not affect the inner eight, but are merely integral functions 
of them. Thus, Ax^, Az^ , Ay j and Aif» can be removed without 
changing the basic modes of motion of the vehicle. The 
first two of the remaining eight states are the primary 
states involved in the phug^id mode, while the next two 
represent the primary short period longitudinal mode. The 
four lateral-directional states often exhibit a Dutch roll 
oscillatory mode and roll and spiral convergence modes. 

Except in special cases, each of these modes involves most 
of the lateral-directional states. 
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The formation of models of order less than eight 
depends on an examination of the individual problem. For 
example, when there is no coupling between lateral-direc- 
tional and longitudinal modes, the eighth-order model can 
be split into two independent fourth-order models with no 
loss of accuracy. If the time-span of interest is short, 
it may be possible to neglect the slower modes (phugoid 
mode and spiral convergence) without decreasing the accu- 
racy of the results; however, as the eigenvectors of earlier 
chapters show, one runs the risk of missing significant 
coupling effects when "inner eight" states are eliminated. 

A different order reduction suggests itself when 
the generalized trim problem is examined. As discussed in 
Section A. 3. 2, the goal of the generalized trim procedure 
is to find a nominal flight condition with constant velocity 
and angular rates. Because they do not affect the velocity 
and angular rate states, Xj , y-j , Zj and ^ are dropped 
immediately. To completely control the six desired states, 
six controls are necessary, but most atmospheric flight 
vehicles have less than six control effectors. Noting that 
the two Euler angles, 6 and <j) , are involved primarily in the 
slow modes, these two states may be regarded as parameters. 
This results in a problem of four control throttle, ele- 
vator, aileron, rudder), two parameters ( t* , . ) , six states 
(u,v,w,p,q,r) , and six state equations (u , v ,w , p ,q ,r ) to 
define trim. 
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A. 4 TOOLS FOR LINEAR ANALYSIS OF AIRCRAFT STABILITY 
AND CONTROL 

The previous sections of this chapter have developed 
the linear aircraft model and demonstrated its validity along 
lighly dynamic trajectories. The use of a linear model is 
desirable because the large body of theory and experience 
relating to the analysis and control of linear systems then 
can be applied to the departure prevention problem. These 
tools are discussed in this section. 

A. 4.1 Eigenvalues, Eigenvectors, and Normal Modes 

The initial-condition response of a linear-time- 
invariant system is composed of a linear combination of a 
limited number of natural, or normal modes . Each normal 
mode is characterized by its time scale, giver, by the eigen- 
value of that mode, and the relative involvement of each 
state in that mode, indicated by the eigenvector of that 
mode. In physically realizable systems, the modes are 
described either by individual real eigenvalues or by pairs 
of complex eigenvalues. A first-order mode exhibits either 
an exponentially increasing response (positive eigenvalue) 
or an exponentially decreasing response (negative eigen- I 

value). A complex (second-order) mode oscillates at a 
frequency determined by the imaginary part of the eigen- 
value within an exponential envelope determined by the real 
part of the eigenvalue. Therefore, the oscillation can 
diverge, converge, or maintain constant amplitude. 

As an example of the eigenvalues involved in the 
normal modes of fighter aircraft, Fig. A. 4-1 illustrates 
th~ areas in the complex pLane which contain the five modes 
anc eight basic eigenvalues of a small high-performance 
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Figure A. 4-1 Approximate Root Locations of a 

High-Performance Aircraft 

aircraft. Because complex eigenvalues occur in pairs of 
complex conjugates, the lower half of the complex plane is 
symmetric with the upper half and is not shown. 

The eigenvectors indicate the relative involvement 
of the aircraft states in a given mode. Each of the modes 
shown in Fig. A. 4-1 may involve the motion of every state, 
but the following generalizations can be made for straight- 
and-level *f light . The longitudinal states (A8, Au, Aq , Aw) 
exhibit two second-order modes. The phugoid mode is a slow, 
lightly damped interchange of kinetic energy (speed) and 
potential energy (altitude) and primarily involves AO and 
Au. The short period mode is the rapid, well-damped angular 
oscillation, and is exhibited primarily by Aq and Aw. 
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The lateral-directional modes usually consist of a 
second-order and two first-order modes. The former is called 
the Dutch roll mode and is a fast, poorly damped yaw oscilla- 
tion about the stability z-axis. The roll convergence mode 
is a fast, stable mode, and it represents the aircraft response 
(generally about the stability x-axis) to a roll moment. 

Due to the angle between body and stability axes, the Dutch 
roll mode and the roll convergence modes appear in the three 
states Av, Ar, and Ap. The spiral mode is slow and frequently 
unstable. An unstable spiral mode is important in a piloted 
aircraft only if the time constant is so short that the 
pilot has difficulty keeping the wings level. 

These mode shapes change considerably as the flight 
condition varies from srtraight-and-level flight. Asymmetric 
flight conditions result in coupled longitudinal/lateral- 
directional modes, and angular rates cause large changes in 
the eigenvalues. In certain cases, modes combine: the 

roll and spiral modes can form a roll-spiral oscillation , 
for example. 

The eigenvalues of the outer four states (AXj, Ayj, 
and At|0 are zero, i.e., these states have no effect on the 
other variables and are pure integrations of the other 
state variables. As a consequence, control-loop closures 
have no direct effect on these modes unless the outer vari- 
ables are fed back directly. 

The eigenvalues are the roots of the characteristic 
equation of the system dynamics matrix, F, 

| XI - F| =0 (A. 4-1) 

where I is the identity matrix of the same order as the 
system matrix, and A is a scalar which must equal an eigen- 
value for Eq. (A. 4-1) to be satisfied. 
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As shown in Ref. 71, there is a set of vectors 
associated with the eigenvalues which have special. proper- 
ties. The eigenvectors , z^, are linear combinations of the 
elements of the state vector and are solutions to the equa- 
tions 


[\ ± l - F]z ± * 0 (A. 4-2) 

As in the case of eigenvalues, the eigenvectors of a second- 
order mode appear as complex conjugate pairs. The eigen- 
vectors contain the same information about the normal modes 
that is given by the time vectors of classical aircraft sta- 
bility analysis (Ref. 65). 

The modal matrix is the matrix of eigenvectors 
arranged columnwise. The inverse of this matrix transforms 
the state vector into normal mode space, in which each ele- 
ment of the vector, A^, represents a normal mode of the sys- 
tem: 

A£ = M _1 Ax (A. 4-3) 

(Instead of the two complex- conjugate eigenvectors of a 
second-order mode, it may be useful to use two real vectors, 
one composed of the eigenvector real part and one the imag- 
inary part . ) 

The linear equation of motion, Eq. (A. 3-3) can be 
transformed as well: 

a£ * M -1 FM A£ + M" 1 G Au (A. 4-4) 

The normal-mode system matrix, M -i FM , is composed of first- 
cr second-order diagonal blocks containing the system eigen- 
values, and the normal mode input matrix, M -1 G, indicates 
which inputs affect which normal modes. This alternate form 
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of the state equation is useful because it demonstrates a 
method of analyzing normal mode excitation. The excitation 
due to the state initial condition can be calculated from 
Eq. (A. 4-3), while the excitation due to control inputs is 
given by the normal mode input matrix, M _1 G. 

This section has discussed eigenvalue and eigen- 
vector concepts. Examples of the application of these 
analytical tools to the analysis of linearized aircraft 
models are contained in Sections 2.4 and 2.5. 

A. 4. 2 Controllability 

In a multi-input, multi-output system, certain 
normal modes may be unaffected by the system controls with- 
out this being apparent from the system dynamics and input 
matrices. This can not occur in an nth-order linear-time- 
invariant system whose controllability test matrix, T, has 
full rank: 

r » jV; FGj . ..j F n “ 1 Gj (A. 4-5) 

The presence of controllability is necessary for 
the construction of a complete system controller, and this 
property almost always exists in physical systems of inter- 
est. (Controllability tests show that the high-performance 
fighter investigated in this report is controllable through- 
out the range of flight conditions.) 

Of more interest is the investigation of control 
effectiveness throughout the flight regime. The difficulty 
is in devising a simple measure of control effectiveness, 
but this can be overcome, to some extent, by using the 
normal mode control input matrix, M”*G, which was intro- 
duced above. The rows of this matrix indicate the relative 
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importance of the aircraft's controls in affecting each 
of the normal modes. 

A. 4. 3 Transfer Functions 


f 

/ 


Specific input-output relationships in linear-time- 
invariant dynamic systems can be described by transfer func- 
tions, which are typically given as ratios of polynomials in 
the Laplace operator, s. The Laplace transform of the or- 
dinary differential equation of motion (Eq. (A. 3-3)) is 

(Sl-F) Ax(s) « G Au(s) (A. 4-6) 

where I, F, and G have been defined, and Ax(s) and Au(s) 
are Laplace transforms of the state and control vectors, 

AxCt) and Au(t). The input, Au(s), and the output, Ax(s), 
are related by a transfer function matrix, H(s), which is 
obtained when Eq. (A. 4-6) is pre-mult iplied by the inverse 
of (sI-F): 

Ax(s) » (sI-F)” 1 G Au( s) 

(A. 4-7) 

* H(s) Au(s) 

where 

H(s) - (sI-F)” 1 G (A. 4-8) 

Any scalar transfer function of intertst (for example, the 
effect of the i*k control on the motion variable), can 
be obtained from Eq. (A. 4-8) by evaluating two determinants 
derived from the matrices of Eq. (A. 4-6) (Ref. 72), 


AXj(s) |sI-F+g iJ j - | sI-F | 
Au i< 8) “ |sI-F| 


(A. 4-9) 
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where is an n by n matrix whose elements are zero, except 
for the j th column, which contains the i tl1 column of the G 
matrix. This transfer function is a ratio of polynomials in 
s, and the numerator and denominator can be factored to iden- 
tify the poles, p, and zeros, z, which describe the relation- 
ship between Au^(s) and AXj(s): 

. Ki(s- 2 1 )(s-z 2 )....(s-z m ) (A.4-10a) 

Au i (s) ~ (s-Pj^Xs-pg) (s-p n ) 


Alternatively, dividing by the individual poles and zeros, 
Eq. (A.4-10a) becomes 


Ax i (s) 

Auj(s) 



(A.4-10b) 


The poles of the transfer function are the roots of the sys- 
tem's characteristic equation, i.e., they are the system's 
eigenvalues, and they are identical for all transfer func- 
tions of the system described by F. The zeros depend on G 
as well as F; therefore, they vary from one transfer function 
to the next. The transfer function gain, Kp, is the steady- 
state value of the transfer function after all transients 
damp out, assuming that all transients are stable. The trans- 
fer function gain, , is (for most aircraft transfer func- 
tions) the initial state rate response to a control step. 


The initial value gain , Kj, is important because it 
determines the initial slope of a given state variable's 
step response. This can be seen by appTying the initial 
value theorem (Ref, 11) to the state variable's transform, 
which states that 
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lim x(t) = lim s x(s) 
t -*• 0 s -*■ 00 


(A. 4-11) 


if the limit exists. The Laplace transform for the pertur- 
bation state, Ax^, given a unity step input in control, 

Au j , is given as 


Ax i (s) 


Ax^s) 

AUj(s) 



(A. 4-12) 


Since there are more poles than zeros in the transfer func- 
tions of interest, the initial value of Ax i (t) is zero. 
Physically, this is an indication that the vehicle states 
do not change instantaneously in response to a control in- 
put. The lowest order non-zero state derivative is equal 
to the excess of poles over zeros. For most aircraft trans- 
fer functions of interest, the excess is one, and 


AXj(s) 


Ax i (s) s Ax i (s) 
s Ax j (s) = Au,(s) s = Au.(s) 


(A. 4-13) 


This leads to the calculation of the initial value of the 
state derivative response to a unity step input as 


Ax i (t*0) = Kj = lim Ax i (s) 

S-*-°° Au.(s) 

J 


(A. 4-14) 


Coupled aircraft transfer functions typ ' jally have 
seven zeros and eight non-zero poles, so Kj and Kp are re- 
lated as follows: 


Kp( -p i ) ( -p 2 ) ••• <-F 8 ) 

K I * (- Zl )(-z 2 ) . . . (-z 7 ) 


(A. 4-15) 


163 


A comparison of the signs of Kj and Kp, which are 
related by the signs of the poles and zeros, as indicated by 
Eq. (A. 4-15), is important because these signs give an indi- 
cation of the expected transient response. If the transfer 
function is stable and minimum phase, the signs of Kj and 
Kp are the same. The resulting response is sim 'lar to the 
solid line in Fig. 2.5-1. A nonminimum-phase zero causes 
Kj and Kp to have opposite signs, and the resulting response 
resembles one of the dashed lines in Fig. 2.5-1. In these 
responses, the initial response direction is away from the 
desired final value. 

The transfer function has been a fundamental tool 
of control system design in the past, and, although linear- 
optimal control theory serves that purpose in this report 
(Chapter 4), transfer functions can be valuable for under- 
standing details of the aircraft's dynamics. For example, 
nonminimum-phase zeros and sign changes in Kj can degrade 
handling qualities (Section 2.2). When poles and zeros 
are nearly equal, there is a cancelling effect which tends 
to remove the associated normal mode from the output vari- 
able’s response to the given input. Conversely, feedback 
paths between the transfer function’s output and input have 
negligible effect on that normal mode. In other words, the 
transfer function provides the information regarding the 
quality of controllability which was missing in Eq. (A. 4-5). 
These capabilities are put to use in Chapter 2. 

A. 4. 4 Optimal Control Theory 

A regulator is a feedback control law which is 
designed to maintain asymptotically stable output of a 
dynamic system, i.e., it bounds the fluctuations in the 
output, and it assures that the output goes to zero as time 


164 



I 


I 


I 

i- 


■\> 




increases. An optimal regulator minimizes a cos; (or 
penalty) functional of the output and control in stabilizing 
the dynamic system. A linear-optimal regulator minimizes 
a particular cost function — the time integral of quadratic 
functions of the output and control — for a linear dynamic 
system, and it takes the form of Eq. (4.1-1) (Ref. 50). A 
linear-optimal regulator can be designed for an aircraft 
near, at, or beyond its open-loop departure boundary. This 
design indicates the control loops which must be closed 
(either automatically or by the pilot) to prevent departure, 
providing asymptotic stability and minimizing a quadratic 
cost functional of the output and control. 


The basic design objective for the linear-optimal 
regulator is to define the feedback control law which mini- 
mizes a quadratic cost functional, J, of the perturbation 
output vector, A^(t), and the perturbation control vector, 
Au(t) : 


J “ f f*Ay T (t) QA^(t) + Au T (t) RAu(t)l dt 
■'0 L ~ -j 


(A. 4-16) 


The control vector contains all available aircraft control 
displacements — in this case, throttle setting ( A£ T ) , ele- 
vator (Afijj), aileron (A6 ), and rudder (A5 ). The output 
vector represents the measured aircraft variables and car be 
formulated as a linear combination of the aircraft's pertur- 
bation states, Ax ( t ) , state rates, Ax(t), and controls, Au(t): 
the present development uses the simplifying assumption, 

A^(t) » Ari ( t ) , where 


Ax T * £a0 Au Aq Aw Av Ar Ap A$J 


(A. 4-17) 


The state-weighting matrix is nonnegative-definite and 
symmetric . 
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(A. 4-18) 


and the control-weighting matrix is positive-definite and 
symmetric : 



r il 

0 • • • 

c n 


R = 

0 

• * 

T* • • • 

22 

0 

• • 

(A. 4-19) 
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-15) can be 

written as 

J = 

a 

q ll Ax l (t) + 


+ q Ax 2 (t) + r. - Au 2 (t ) + 
M nn n v ' 11 l v 1 



+ r Au^( t ) 
mm m 

]« 

(A. 4-20) 


and the cost functional is seen to be a weighted sum of the 
integrated-square values of the perturbation state and control. 
In the present case, minimizing the weighted sum of integrated- 
square values is equivalent to minimizing the weighted sum of 
root-mean-square (rms) values of the state and control . 

Equation (A. 4-20) provides a means of trading off 
the cost oi output errors against the cost of control, and 
it is simply this: choose each \ve:ghting coefficient in Q 

and R as the inverse of the maximum allowable mean-square 
value of the weighted variable, i.e., 
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(A. 4-21) 


Oil ’ 1/ax I 


max 


r . . = 1/Au , j ■ 1 to n 
JJ J max 


(A. 4-22) 


This normalizes each term in Eq. (A.4--20), s.- that us contri 
but ion to the integrand is unity when the variable equals its 
maximum value. The elements of R are specified by the con- 
trol authority which can be assigned to the DPSAS. For 
example, if 10 deg of elevator perturbation can be assigned 

to departure prevention, the corresponding element of R is 
2 

1/(10) = 0.01. If there is prior information regarding 

allowable state perturbations (as in a tracking task), the 
elements of Q are determined similarly. 


An alternate approach is to use the elements of Q 
and R as design parameters which can be varied until desirable 
transient response or eigenvalues are achieved. In such case, 
the equivalence of Q elements to allowable mean-square values 
is not lost, and it is possible to gain insight regarding the 
correspondence of rms-output errors and classical figures of 
merit in each particular case. 


The minimization of J must be accomplished subject 
to the dynamic constraint provided by the linear equation 
of motion, 

Ax ( t ) = FAx(t) + GAu(t) (A. 4-23) 

(It is assumed that F and G form a controllable pair.) The 
method of finding the control which minimizes J subject to 
Eq. (A. 4-23) is derived in numerous texts (e.g., Refs. 60 to 
62). In the special case of quadratic cost and linear sys- 
tem dynamics, the control solution is a linear feedback law 
(Eq. (4.1-1)). The gain matrix of this control law is 



-1 T 

K = R X G P 


(A. 4-24) 


where the symmetric matrix, P, is the steady-state solution 
of the matrix Riccati equation 

P = -PF - F T P + PGR -1 G T P - Q (A. 4-25) 

In other words, the DPSAS gain matrix is easily 
found by two mat. ix multiplications once P is found 

(Eq. (A. 4-24), but the solution for P (Eq. (A. 4-25)) appears 
formidable. There are, however, four recognized methods 
for steady-state solution of Eq. (A. 4-15)), all of which 
require digital computation (Ref. 60): these are direct 

integration of Eq. (A. 4-25), the Newton-Raphson method, the 
Kalman-Englar method, and the diagonalization/eigenvalue 
method. The choice between these methods must be based on 
grounds of numerical convenience and efficiency. 


The Kalman-Englar method (Ref. 73) has been used to 
generate the results which follow in later sections. In 
this technique, P is propagated to steady state using the 
recursive equation 


- ■ [ 0 
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(A. 4-26) 


The matrices 0^ , G^ 9 , 
(n x n) sub-matrices of 


and 0 22 


are 


tne appropriate 


©(At) 


-ZAt 

e 


where 


Z = 



-1 T 
-GR i G i 



(A. 4-27) 


(A. 4-2 i) 
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and the propagation interval, At, is small compared to the 
natural periods of the aircraft motion. 

The DPSAS design procedure is summarized and shown 
to be a straightforward technique once the appropriate 
general-purpose computer routines are programmed: 


• Define control authority available 
to the DPSAS, thus specifying R. 

Define allowable state perturba- 
tions, thereby specifying Q. 

• For the aircraft dynamics specified 

by the stability and control matrices, 

F and G, compute the feedback gain 
matrix, K, using Eq. (A.4-2&), (A. 4-27), 
(A. 4-26), and (A. 4-24). 

• The control law for the DPSAS is given 
by Eq. (4.1-1), using the gain matrix 
calculated in the previous step. 


The resulting DPSAS stabilizes the aircraft without using 
mroe control authority than that specified by R for state 
perturbations defined by Q. 


A . 4 . 5 Gain Scheduling Procedure 

Means and Standard Deviations - Two features which 
suggest that a gain be held constant are its standard devi- 
ation and mean value. Certain gains do not exhibit wide 
variations as the flight conditions change. This can be 
determined by constructing a table of means and standard 
deviations for the gains, as illustrated by Table A. 4-1. 

In the table, Gain 6 displays a low standard deviation and 
a large mean value. This indicates that the gain should 
probably not be scheduled, i.e., that its mean value can be 
used at all flight conditions. 
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Another feature of each gain is its relative mag- 
nitude compared to other gains of a similar class. Gain 4 
is small compared to Gain 6 and also exhibits a wide varia- 
tion in magnitud-; thus it may be desirable to schedule 
Gain 4, if its variation with flight condition is coherent, 
or set it to zero. Simulations should be done to fully 

TABLE A. 4-1 

EXAMPLE OF MEAN-STANDARD DEVIATION TABLE 

Standard S.D. Per Cent 
Mean Deviation of Mean 


Gain 

4 

-.00598 

.0269 

449.0 

Gain 

5 

-.1013 

.032 

31.7 

Gain 

6 

2.55 

.046 

1.82 

etc 






determine the zeroed gains' effects on controlling the air- 
craft. Gain 5 is a logical candidate for scheduling. The 
gain magnitude is not negligible, and it displays enough 
variation to warrant scheduling. 

Correlation Between Gains and Flight Variables - 
The aircraft dynamic model varies in a complex but deter- 
ministic way with flight conditions. If the closed-loop 
response of the aircraft is maintained essentially invariant 
by automatic control, it is reasonable to assume that the 
necessary control gains also vary in a complex but deter- 
ministic way with flight conditions; hence the gains and 
flight conditions should be correlated, 

The search for gain/flight variable dependencies 
begins by determining correlation coefficients between the 



I 


gains and all available flight variables. One method of 
determining the correlation coefficient between a set of 
gains (dependent variables), and a flight variable 

(independent variable), m, is given by the following: 



In Eq. (A. 4-29), I is the number of flight conditions for 
which the gains are known, ^ is the value of the gain 
observed at flight , oint i, and tr^ is the value of the 
flight variable at flight point i. The variable £ is the 
mean value of the gain, given by 

I 

k = ~ 2 k. (A. 4-30) 

i=l 

and m is the mean value of the flight variable. The closer 
the magnitude of p is to one, the better the correlation 
between the gain and the flight variable. Independent 
variables ^hich can be considered for gain scheduling 
include indicated airspeed (IAS), body-axis velocities 
(u,v,w) angles of attack and sideslip (a, 8), angular rates 
(p.q,r), and control trim positions. These variables can 
be inverted, squared, and so on, in the search for high corre- 
lation. An example of a correlation coefficient ta^le is 
shown in Table A. 4-2. The circled values are the high corre- 
lation coefficients betv^en gains and independent variable 
functions . 
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TABLE A. 4-2 

EXAMPLE OF CORRELATION COEFFICIENTS TABLE 




a 

a* 

1/a 

IAS 

IAS2 

P 

q 

r . . .etc . 


Gain 

1 

.838 

(939) 

.297 

(93$) 


.804 

.018 

.0359 


Gain 

2 

.879 

.743 

.168 

.745 

.541 

(890) 

.028 

.245 

i 

Gain 
etc , 

3 

.672 

.786 

.069 

.782 

(904) 

.724 

.447 

.339 

i 

\ 


Curve Fitting - The third step in the gain-scheduling 
procedure is to construct a smooth relationship between the 
gains which are to be scheduled and the most highly corre- 
lated flight variables. Multiple regression 

k » bg + b i n i + tj 2 m 2 + ••• + kn m n (A. 4-31) 

and polynomial regression, 

k = bQ + b^m^ + ^2 m l + • • • * t> Q m° (A. 4-32) 

can be used. Equation (A. 4-31), the multiple regression, 
uses n different independent variables, while Eq. (A. 4-32), 
the polynomial regression, uses powers of the highest cor- 
related independent flight variable, n^, to estimate the 
gain, k. A method for determining the regression coeffi- 
cients, b^, is shown next. 

A multiple regression analysis determines the 
regression coefficients bQ, b-,...,b in Eq . (A. 4-31) so 
that the sum of the squared error between the regression 
estimate, k, and the tiue value of k is minimized. For 
I flight conditions, the function to be minimized is 


i 
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V. 



(A. 4-33) 


J - £ (k.-k .) 2 

i=l 1 


To minimize J, set 


9J _ 9J_ = M_ - a (A. 4-34) 

9b n 9b- " • • ' 9b 
0 1 n 


Then the expression for the regression coefficients becomes 


A b = C 
or 

b = A -1 C 


(A. 4-35) 
(A. 4-36) 


where A and C are defined as in the least-squares formulas 
of Ref. 63. These are the values of b^ which minimize the 
squared error. 


To determine just how good the b^ values are, the 
correlation coefficient for the multiple regression fit can 
be found, as in Eq. (A.4-°9): 



(A. 4-37) 


The closer p is to one, the better the fit of the multiple 
regression model. When using & multiple regression, the 
more independent variables chosen, the higher the value of 
p will be, until n=I and p=l. 
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An alternate way of estimating gain values is to 
use only one of the flight variables in a polynomial regres- 
sion. An n^k-order polynomial regression analysis deter- 
mines the regression coefficients, b Q , b^,...,b n in Eq. 

(A. 4-32) so that the sum of the squared error between the 

/N 

regression estimate, k,of Eq. (A. 4-31) and the true value 

of k is minimum. In this case, there is only one kind of 

independent variable, m. , for each gain value, but it is 

1 . 2 

raised to various powers, i.e., m-^ * nK , mg = , until adequate 

correlation is achieved. 

The analysis for the polynomial regression pro- 
ceeds as in the multiple regression, starting at Eq. 

(A. 4-33). The polynomial regression can be considered a 
special case of the multiple regression. 

A . 4 . 6 Pr ogram ALPHA 

The construction of the complete linear equations 
of motion, their analysis, and the design of feedback con- 
trollers has been programmed in ALPHA — Analysis Program 
for High Angle-of-Attack Stability and Control. Figure 
4.4-2 illustrates the structure of this program. Input 
data consists of aircraft inertial and aerodynamic char- 
acteristics. The aerodynamic data can be entered as con- 
ventional stability derivatives, dimensionless derivatives, 
or full tables of nonlinear force and moment characteristics. 
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Figure A. 4-2 ALPHA - Analysis Program for High 

Angle-of-Attack Stability and Control 

The program executes a three-step procedure. The 
first step consists of steady or generalized trim calcula- 
tion, if desired, and the construction of the complete body- 
axis system dynamics and control input matrices. During 
the second step, the linear system is modified, if required, 
to include any axis transformation, order reduction, or fixed 
stability augmentation loop closure. The final step consists 
of the analysis of the resulting system. Eigenvalues, eigen- 
vectors, transfer functions, linear-optimal stability aug- 
mentation systems, and time histories can be calculated and 
plotted. 
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The executive structure of ALPHA includes logic 
to modify the dynamic model on succeeding passes through 
the program and to vary the analysis type as certain param- 
eters are varied over a given range of interest. Program 
ALPHA provides an efficient tool for the thorough analysis 
of aircraft high angle-of-attack stability and control. 
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APPENDIX B 


AIRCRAFT AERODYNAMIC MODEL 

The reference aircraft is a small, supersonic fighter 
type designed for air superior *y missions. Mass, dimensional, 
and inertial characteristics are listed in Table B-l. The 
aerodynamic data set is a composite of sub-scale wind tunnel 
measurements for two configurations of the reference aircraft; 
hence, the numerical results presented here do not represent 
a specific aircraft in detail. 


TABLE B-l 

CHARACTERISTICS OF THE REFERENCE AIRCRAFT 


Mass , kg 6 > 124 • 

l x , kg-m 2 7,186. 

l y , kg-m 2 57,350. 

l z , kg-m 2 63,450. 

I xz , kg-m 2 “149. 

Reference Area, m^ 17.30 

Mean Aerodynamic Chord, m(c> 2.46 

Wing Span, m(b) 

Length, m -73 

Reference Center of Gravity (c.g.) 0.25c 

Flight c.g. 0.17c 



The control variables are elevator (or horizontal 
tail, 6^), leading/trailing edge flaps (6 f ), ailerons (6 a ), 
rudder ( 6 r ) , speed brake (6 SB ),and ‘ irust setting ($ T ). 

The ranges of these variables are listed in Table B-2. 


TABLE B-2 

CONTROL VARIABLE RANGES 


6 

6 

6 

6 

6 

<S 


h 

f 

a 

r 

SB 

T 


-20 to +5 deg 

0 to 100 percent 
-60 to +60 deg 
-30 to +30 deg 
0 to 45 deg 
0 to 100 percent 


The aerodynamics of the aircraft are represented 
by 45 coe 'icients which are functions of angle of attack, 
sideslip angle, elevator deflection, and flap setting. 

Using simplified conventional notation, the 6 total coeffi- 
cients are described as follows: 



C 


C x (a,B) + AC^a • B ) 4 AC^(a , 6g S ) - AC (a , 6 f , 6 h ) 
+ C x (a) (B-l ) 

q 

Cy ( a , 6 j )+ACy( a , S ) + ACy(a , 6 , 6 f ) + ACy(oi,6 f ,6 a ) 

+ AC y (a,6 r ) +^ [c Yr (a)r + C yp (a)p + C Y £(a)eJ 


(B-2) 
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C„, = C z (a . 8) + AC z (a , B , ; + AC z (a , 6gg) + AC z (a , 6^ 6^) 


+ if c z <“> 

q 


(B-3 ) 


* C 1 (a , 6 f ) + C-^cx, 8) + AC-^a, 8, <5 f ) + AC-^a, 6 f , 6^) 


+AC 1 (a,6 r ) + 


2V [% 


(a)r + C 1 (a)p + C 


H Mi ] 


(B-4) 


%. " C m 0 +C m <«'< , ) +C Z„(^) + '' c n ,( = .S,8 f )*JC |r U,6 f ,5h) 


+/ V“' 6 SB )+ ff C m 

q 


( B-5 ) 


C - C n (o,6.) +C n (o,B) + AC„ (Ax) g + AC (a .8,6.) 
n^nQin 1 T b n f 


+AC (a ,6., 5 ' + AC (a, 6 ) 
n v f ’ a n r 


^ [ C n r ( ° 


(a)r + C («)p*C 
r p 8 




(B-6 ) 


The first step in finding the individual terms of 
the perturbat: forces and moments, Eqs. (A. 3-13) to (A. 3-16), 
is to evaluate the derivatives of the force and moment co- 
efficients with respect to the nondimensional states at the 
nominal flight condition. This results in nondimensiona 1 
stability derivatives, such as the following: 
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Many of these derivatives contain the partial deriva- 
tives of the nondimensional wind-axis translational velocities 
(V/Vq, 8, a) with respect to the nondimensional body-axis 
translational velocities (u/Vg, v/Vg, w/Vg). This matrix of 
derivatives, evaluated at the nominal flight condition, is 
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The dimensional stability derivatives are formed 
by taking the derivatives of the dimensional aerodynamic 
forces and moments with respect to the dimensional state 
variables. These dimensional derivatives contain the non- 

3Fy 3Fy 

dimensional derivatives; and are examples of these 

derivatives : 
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The complete dimensional stability derivative 
matrices are 
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These stability derivative matrices are used in Eqs. (A. 3-13) 
to (A. 3-16) to determine the perturbation forces and moments, 
which themselves determine the aerodynamic terms in the linear 
perturbation equations of motion. 
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